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Since the May meeting of the Administrative Com- 


mittee there has been considerable activity directed 


toward the selling of advertising and otherwise relieving 
the financial problems of PGAP. In particular, there was 


a meeting at Boulder, Colo. during the URSI General 


Assembly which is briefly reported below under news 
from that chapter. To date, over $6000.00 in firm com- 


_ mitments have been received for advertisements to ap- 
pear in this and the next three issues. This represents ap- 


a 


proximately 60 per cent of the total estimated amount 


‘needed in addition to our anticipated income from mem- 


bership dues of approximately $12,000.00. It seems 
hopeful that within the next few months we will solve 


our financial problems and the publication difficulties 


to which they give rise. 


Administrative Committee 


The following are the current members of the Ad- 
ministrative Committee and their assignments: 

Elected Members—Terms Expire 1958: J. 1. Bohnert, 
Chairman, H. G. Booker, Library Subscriptions, D. D. 
King, Awards, J. B. Smyth, PGAP Editor. 

Elected Members—Terms Expire 1959: A. Dorne, 
News and Views Editor, F. T. Haddock, Radio Astron- 


~ omy Issue of PRoceEpINGs, J. W. Herbstreit, Technical 


Symposia, R. L. Mattingly, Vice-Chairman and Mem- 


bership. 
Elected Members—Terms Expire 1960: J. W. Findlay, 


1958 National Convention Program, R. K. Moore, 
_ Paper Procurement, Wm. H. Radford, No official as- 


 signment, O. G. Villard, Chapter Organization. 


Ex Officio Members (Past Chairmen): P. S. Carter, 


~ term expires 1959; D. C. Ports, now assigned to Ad- 
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ADMINISTRATIVE NEWS vertising, term expires 1960; A. H. Waynick, now as- 
BE ances signed to Student Membership, term expires 1958. 


Honorary Member: L. C. Van Atta. 

Assisting Members—A ppointed: S. A. Bowhill, Prop- 
agation Papers Review, C. Goatley, Secretary-Treas- 
urer, S. M. King, 1958 NATIONAL CONVENTION REC- 
orp, V. H. Rumsey, Antenna Papers Review, R. C. 
Spencer, PROCEEDINGS-TRANSACTIONS Liaison. 

Comment: Questions or comments regarding any 
phase of the activity of the Administrative Committee 
are, of course, the prerogative of any PGAP member. 
One such letter recently has been received. It is pub- 
lished here, together with the reply that was sent the 
writer, because it concerns an issue which may trouble 
other members. Any further comments will be welcomed 
by this department. 


“Dear Dr. Mattingly: 


“This is to express profound concern over the recent 
‘election’ to fill four vacancies on the PGAP Adminis- 
trative Committee. Such ‘rubber-stamp’ procedures 
as this, in which only four candidates are nominated 
for an equal number of posts, make a mockery of the 
democratic process of voting, and the ballot becomes 
equal to a proxy. 

“Incidentally, it is a constant source of wonderment 
to students from abroad to witness, in many of our 
organizations, such hypocritical practices, which seem 
to conflict with the American ideals we assert, and that 
are very difficult to explain away. 

“It is respectfully suggested that the PGAP By-Laws 
be appropriately modified to allow for genuine elections 


in the future.” 
ROBERT COHEN 


857 18th Street 
Boulder, Colo. 
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“Dear Mr. Cohen: 

“Thank you for your letter of June 22, 1957, in which 
you indicate your concern over the procedures used in 
the election to fill vacancies on the PGAP Administra- 
tive Committee. 

“As I am sure that you realize, it is not the intent of 
the Administrative Committee to usurp the democratic 
birthright of the American citizen. However, we who 
are interested in the PGAP must do what we can to in- 
sure that the PGAP is administered by individuals 
with proven capabilities and work capacity and whose 
business affiliations place them in a position to carry 
out the uninspiring and often arduous business details 
of our organizations. In addition, as explained in the 
note accompanying the recent ballot, we must strive for 
a reasonable balance of representation between the an- 
tenna and propagation areas, between various geograph- 
ical regions and to introduce spokesmen for such special 
fields, as radio astronomy, which we feel fall properly 
within our cognizance. It is just this sort of thinking 
that dictates the apparent arbitrary manner in which 
elections are held in other scientific organizations. 

“As you imply, we could submit a list of, say, eight 
names from which the membership could make their 
choices for the four vacancies. As a matter of fact, we 
have done substantially this in the past. However, it 
has been our experience, and the experience of many 
other similar societies that this procedure does not 
serve the best interests of the group. 

“T hope that you do not imagine that the write-in 
votes go for naught. While it is true that a write-in 
election would be difficult to bring about, write-in 
nominations are carefully tabulated, including the 
single-vote cases, and where there is a definite expres- 
sion for certain candidates, these are considered by 
our nominating committee for subsequent elections. 
Furthermore, our local Chapters have a good deal to 
say in the choice of candidates. As an example, Professor 
Moore, recently elected to the Administrative Commit- 
tee, was suggested by our Albuquerque Chapter. 

“Let me assure you that there is no conspiracy on 
the part of the Administrative Committee to restrict 
admission to this body. The Administrative Committee 
sincerely tries, within the limitation of human frailties, 
to make the choice of the membership-at-large their 
choice, provided that certain necessary qualifications, 
established often by bitter experience, are met. 

“T shall take pleasure in submitting your letter for 
consideration by the Administrative Committee at its 
next meeting. Thank you again for your expression of 
feeling.” 

ROBERT L. MATTINGLY 
Vice-Chairman, PGAP. 


CHAPTER NEWS 
Chicago 
The Chicago Chapter meeting took place September 
13, 1957. Norman J. Foot, Assistant Chief Engineer in 
charge of the Antenna Section of the Hallicrafters Air- 
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borne Laboratory spoke on “Introduction to the Scimi- 
tar Antenna.” The Scimitar is a very broad-band and 
and highly efficient antenna, the impedance and patterns 
of which are nearly independent of frequency. It is 
related to the family of antennas whose shapes are 
specified entirely by angles. The meeting was well at- 
tended. 


Denver-Boulder 


During the two-week period—August 23 to Septem- 
ber 5, 1957—the International Scientific Radio Union 
(URSI) held its triennial meeting in Boulder, Colo. This 
Twelfth General Assembly, was the first such meeting 
held in the United States since 1927. Many of the ses- 
sions of the seven commissions of the Union contained 
material of great interest to the members of the PGAP. 


The occasion drew together many of the national officers 


of both the IRE and the PGAP. To take advantage of 
this, a joint business meeting of IRE executive and 
PGAP administrative committees was held at the 
Boulder Laboratories of the National Bureau of Stand- 
ards on August 29th. The officers of the Denver-Boulder 
section of I[IREand Denver-Boulder Chapter of the PGAP 
participated. 

The attendance was approximately twenty-five and 
included among others, J. T. Henderson, President, 
IRE; Dr. Yasujiro Niwa, Vice-President, IRE, and 
Chairman of the Tokyo Section; Dr. Isaac Koga, Vice- 
Chairman of the Tokyo Section and the newly-elected 
Vice-President of URSI; Dr. J. I. Bohnert, National 
Chairman, PGAP; R. L. Mattingly, Vice-Chairman, 
PGAP; C. Goatley, Secretary, PGAP; R: S. Kirby, 
S. B. Peterson, and W. C. Worcester, Chairman, Vice- 
Chairman, and Secretary, respectively, of the Denver- 
Boulder Section, IRE; and H. V. Cottony and A. P. . 
Barsis, Chairman and Secretary, respectively, of the 
Boulder-Denver Chapter of PGAP. 


Los Angeles 


The Los Angeles Chapter meeting was held on Sep- 
tember 12 at the IAS Building, 7600 Beverly Boulevard. 
The speakers were Dr. Joseph A. Barkson of Hughes 
Aircraft Corporation, who spoke on “Coupling of Rec- 
tangular Waveguides,” and Kenneth C. Kelly of Hughes 
Aircraft, who spoke on “Recent Annular Slot Array 
Experiments.” - 


Philadelphia 


The new officers elected for 1957-1958 are: Chairman, 
Edmond J. Forbes, Radio Corporation of America, 
Bldg. 1-5, Camden, N. J.; Vice-Chairman, Edwin A. 
Fink, Philco Corporation, 4700 Wissahickon Avenue, 
Philadelphia, Pa.; Secretary, Frank T. Barclay, ITE 
Circuit Breaker, 601 E. Erie Avenue, Philadelphia, Pa. 


Washington, D. C. 


The new officers elected for 1957-1958 are Harry 
Fine, Chairman; Herschel T. Ward, Vice-Chairman; 
and Robert J. Adams, Secretary. 
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contributions 


Radiation from a Radial Dipole Through a 
Thin Dielectric Spherical Shell* 
M. G. ANDREASEN} 


Summary—tThe radiation from and the reaction on a radial dipole 
placed inside a thin dielectric spherical shell are investigated. Nu- 
merical results have been obtained for the gain pattern and the 
_ impedance change of an elementary dipole. The nearer the center of 
the shell the dipole is placed, the less the pattern is found to be in- 
fluenced by the shell, the opposite being the case for the impedance 
of the dipole. 


INTRODUCTION 


HE use of a dielectric shell to protect an antenna 
ab against rain water or the action of the wind, as is 

well known, may strongly influence the electrical 
properties of the antenna. Part of the energy radiated 
by the antenna may be reflected back into the line feed- 
ing the antenna, thus changing the antenna impedance 
and the broad-band properties of the antenna. Further- 
more, because of internal multiple reflections from the 
shell, the radiation pattern may be considerably dis- 
torted. 

To avoid strong reflections back into the antenna line 
it seems reasonable from a point of view of geometrical 
optics that the antenna should not be placed near any 
focus of the shell. On the other hand, the radiation pat- 
tern of the antenna is influenced less by the shell if it is 
placed nearer a focus of the shell. Therefore, a com- 
promise usually must be made as to the position of the 
antenna inside the dielectric shell. The electrical thick- 
ness of the shell should be either as small as permissible, 


* Manuscript received by the PGAP, May 16, 1956. 
+ Siemens und Halske, Munich, Germany. 


due regard being given to the mechanical strength of 
the shell, or it should be chosen equal to half a wave- 
length. For a half-wavelength shell, the electrical per- 
formance of the surrounded antenna is influenced less if 
the dielectric constant of the shell is higher. 

It is the aim of this paper to present the results of an 
investigation of the influence of a thin dielectric spheri- 
cal shell on the radiating properties of a radially di- 
rected Hertz dipole arbitrarily placed inside the shell. 
An investigation of the special case of a dipole placed 
at the center of a spherical shell was previously made by 
Keller! and Mirimanov.? The results obtained by the 
present investigation are expected to give a very good 
idea of the influence of a spherical shell on the perform- 
ance of a linear antenna of finite length. The method 
used for solving the problem is the classical one of ex- 
panding fields in spherical harmonics.? This method is 
particularly useful when the radius of the shell is at 
most a few wavelengths, the necessary number of terms 
of the infinite sum representing the radiated field thus 
being reasonably small. However, when the primary 
field is being expanded, difficulties often may be ex- 
pected. But if the interest is confined to the far-field 


1 J. B. Keller and H. B. Keller, “Reflection and transmission of 
electromagnetic waves by a spherical shell,” J. Appl. Phys., vol. 20, 
pp. 393-396; April, 1949. ; 

2 R, G, Mirimanov, “Radiation resistance of a dipole placed in the 
center of a thin spherical envelope,” Dokl. Akad. Nauk, vol. 71, pp. 
1061-1064; 1950. ‘ 

3J. A. Stratton, “Electromagnetic Theory,” McGraw-Hill Book 
Co., Inc. New York, N.Y.; 1941. 
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only, a method‘ based on the reciprocity theorem may 
be used with advantage, the expansion of the primary 
field being reduced to the expansion of a plane wave. It 
should be noted, however, that by using the former 
(direct) method the antenna reactance may be deter- 
mined; this will not be the case when the latter (indirect) 
method is used. 

Because of the general character of the present in- 
vestigation the expressions obtained for the radiation 
through a lossless spherical shell may in a wider sense 
represent also the influence of an absorbent shell, the 
real dielectric constant of the lossless shell being re- 
placed by a complex dielectric constant. This formal 
procedure leads, however, to difficulties in the numerical 
calculations because, as far as is known to the author, 
no tables of spherical wave functions of complex or even 
just purely imaginary argument exist. 


SIMPLIFIED BOUNDARY CONDITIONS AT A THIN SHELL 


In order to determine the exact field transmitted 
through a dielectric shell, in addition to the primary 
field, a secondary field representing a sum of standing 
waves must be assumed to exist inside the shell. Outside 
the shell the field is represented by a sum of diverging 
waves, whereas in the shell the field may be represented 
by a sum of standing and diverging waves. Thus, there 
_are four secondary fields and four sets of constants to be 
determined by satisfying the boundary conditions at the 
two surfaces of the shell. This procedure is necessary if 
the electrical thickness of the shell is great. If the shell 
is electrically thin, however, which is assumed in the 
present investigation, the shell may, to a good approxi- 
mation, be replaced by an infinitely thin shell, the di- 
electric constant of which approaches infinity so that 
the product of the infinitesimal thickness and the di- 
electric constant of the shell approaches a constant 
value. Thus, only two secondary fields need to be de- 
termined. To illustrate this, consider a dielectric shell of 
dielectric constant €, permeability uw, and a thickness d, 
which is assumed to be negligible as compared to the 
smallest radius of curvature of the shell. The shell is 
placed in free space. Assuming the electrical thickness of 
the shell to be small, the discontinuity of the tangential 
magnetic field across the shell is approximately equal to 
the displacement current in the shell. Thus, if we use the 
time factor e#! 


ft X (Hy — Hi) = jwed(n X E,) X f, (1) 


the shell dividing the space into the regions 1 and 2, % 
being a unit vector normal to the shell from medium 1 
to medium 2, and (#X E,) XA~$0 X (Li + FE) XA being 
the mean tangential electric field in the shell. Similarly, 
the discontinuity of the tangential electric field across 
the shell is approximately equal to the magnetic cur- 
rent in the shell 


fi X (FE, — E:) = — joud(h X Hi) X A, (2) 


4P. S. Carter, “Antenna arrays around cylinders,” Proc. IRE, 
vol. 31, pp. 671-693; December, 1943. 
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(AXH,) Xa~tnX(H2+H) XH being the mean tan- 
gential magnetic field in the shell. Now, letting the shell 
shrink to a shell of infinitesimal thickness 6, dielectric 
constant e’, and permeability uy’ the boundary conditions 
(1) and (2) at the original shell are preserved by this 
process when in these conditions we replace ed by 
€)(d—5)+e'5 and ud by po(d—5)+p’d, thus setting 
(e’ —€))5 = (€—€0)d and (u’ —po)6 = (u—mo)d. The bound- 
ary conditions at the infinitesimal shell are then, as we 
pass to the limit 6-0, e’ > ~ and p’—> », 


h X (He — Hi) = jole — eo)d(h X E.,) Xn (3) 
f X (E2 — E1) = — jw(u — po)d(a X H.) X #. (4) 


In what follows we assume u=po. Thus the magnetic 
current in the infinitesimal shell is zero, and the bound- 
ary conditions are reduced to 


AX (H.-W) = VAX E) XA= VAX Ex) XA (5) 
where 


4 _B or \d (a =A) ke 
TSP er RON G5) Pa aes 
us ¢ Jj 0 i ¢ 
is the admittance of the shell, x, being the relative di- 
electric constant of the original shell, and ¢ and k being 
the specific impedance and propagation constant re- 
spectively of free space. The simplified boundary con- 
ditions (5) have been shown to be less approximate, the 
larger the dielectric constant of the actual shell is and 
the smaller the gradient along the surface of the shell of 
the electric field component normal to the shell is when 
compared with the value of the tangential electric field 
at the shell. 


(6) 


PRIMARY FIELD 


For the purposes of the following calculations a spheri- 
cal coordinate system (r, 0, ¢) (Fig. 1) is introduced, the 
dipole being placed at the point (x, y, s)=(0, 0, 70) in- 
side an equivalent infinitely thin shell r=7,. The fields 
produced by the radial dipole are evidently of the TM- 
type and are independent of @. Thus a magnetic vector 
potential® ( 


A = rA? (7) 


may be introduced appropriately, VX A being equal to 
the magnetic induction B. The scalar potential A must 
satisfy the wave equation (V?+2)A =0, the solutions 
of which determine the electromagnetic field vectors 
through the spherical components 


d°(r A) 
jopel, = + k?(rA) 
or? 

. 1 d*(rA) 
Jupelg = — 

r  oroé 
ea 1 0d%(rA) 
wey, = 8 
q * sin 6 ordd (8) 


°P. Frank and R. v. Mises, “Die Differentialgleichungen der 
Physik,” Friedr. Vieweg & Sohn, Braunschweig, 2nd ed.; 1935. 
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(r, 2,9) 


Fig. 1—Coordinate system used for the calculations. 


HoT 
1 ao A 
pH, = aaa 
sin @ d¢ 
= 0A (9) 
BK = —_—_ -—— se 
“ a6 . 


_ The scalar potential A>, of the field radiated by a 
radial current distribution in free space must satisfy the 
_inhomogeneous wave equation 


boy 
? 


r 


Yih i) App (10) 


J, being the current density. Integrating this equation 
for the elementary dipole we obtain for the potential at 
the point (7, 0, @) (see Fig. 1) 

po a PAL se s%* 


ee 11 
ss 4r To R ( ) 


_I being the prescribed current of the dipole and A/ the 
length of the dipole. This expression is expanded in 
spherical harmonics,’ giving 


cS 
Apr = — ho (RR) 


To 


to} 


is; 
= >> — (2m + 1)jn(Rro)hn (kr) Pn(cos 6), (12) 
n=0 10 
which is valid for y>70; and where for simplicity a dipole 
constant C has been introduced, 
oR 
ea a 
4a 


(13) 


The notations used here for the spherical Bessel function 
jn and for the spherical Hankel function of the second 


kind h,2 are in accordance with those used by Morse.° 


P,,(cos @) is the Legendre polynomial of mth order. 


6 P. M. Morse, “Vibration and Sound,” McGraw-Hill Book Co., 
Inc., New York, N. Y.; 1936. 
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TRANSMITTED FIELD 


The magnetic potential A, of the field inside the shell 
is taken to be the sum of the primary potential and a 
secondary potential representing standing waves 


A, = A pr = A sex 


ino} 


C 
2 E (2m + 1)jn(Rro) ln (kr) + anja) | 


n=0 


- P,(cos 6). 


(14) 
Outside the shell the potential must represent diverging 


waves, thus 


Ay = D5 dahn® (kr) P,(cos 6). 


n=0 


(15) 


The constants a, and 6b, are to be determined so as to 
make the resultant field satisfy the boundary condi- 
tions at the shell. Inserting (8) and (9) into (5), these 
conditions are seen to be given by 


DG) a). dennis Gee Dee Phan S we oe 
Making use of the Wronskian 
_,« Aphn®(p)] dl pjn(p) 
Di Og ie Ca) Heino) cate epee) 
dp dp 


We derive from (14)—(16) the constants a, and 6, to be 


G 
bn = — (20 + 1)jn(kro) Rn (18) 
To 
dl kryhtn® (kr1)| i 
= — j(x. — 1)kd bn 19 
where for simplicity a constant 
R : (20) 
rerio d[krijn(kri)] dlkriln®(Rrs)] 
1+ j(ke — 1)kd 
d(kr1) d(kr;) 


depending only on the parameters of the spherical shell, 
has been introduced. 

The field transmitted through the dielectric shell may 
now be found by introducing (15) into (8) and (9). Mak- 
ing use at the same time of the differential equation for 
spherical Bessel functions,’ 


p°Zn!’(p) + 2pzn'(p) + [o? — n(w + 1)]2n(o) = 0, (21) 


the nonzero components of the transmitted field are 
found to be 


Ge hn (Rr) 
Eyy= a > dpn(n + 1) P,,(cos @) 
JHo n=0 kr 
s 1 dlkrh»®(kr)| dP,(cos 6) 
Pees 3s 
Jeo n=0 (2 ii dé 
P,,(cos @) 
Hy = — — B53 b In ®(br) = —. (22) 
Ho n=0 dé 
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This represents the field at any point (7, 6, @) outside 
the spherical shell. For practical purposes the interest 
attaches primarily to the far field. Inserting the asymp- 
totic expression for the spherical Hankel function of the 
second kind,’ 


—ip 


he? (o) = grt" p>n, (23) 
we find the far field (22) to be given by 
pom aaa dP,,(cos @) 
imi 24 Ra eaeaaeaneanse 
bi Ho kr di dé 
5 a dP,,(cos 0 
Hy = —-— YS jd, dP, (cos 8) : (24) 
Ko Kf nm d0 


The r component of (22) has been neglected as this com- 
ponent is small of second order. In order to discuss these 
results the antenna gain will be calculated in the next 
section and will be made the object of numerical com- 
putations so as to give an idea of the influence of a 
spherical shell on the radiation pattern of an actual 
dipole. 


ANTENNA GAIN 


The antenna gain of an antenna system of rotational 
~ symmetry is defined by 


P09) _ 


Ku 


Gib) = 4m P(, ¢) 


Z » (25) 
if P(0, ®) sin 6d 
0 


P(6,) and Pr/47 being the actual and the mean power, 
respectively, transmitted per unit solid angle in the di- 
rection (0, ). To find the antenna gain of the radial di- 
pole placed inside the spherical shell P(@, ) is replaced 
by P(@, ) =r? Ex9| 2/2¢, Eos being the far-zone electric 
field from (24). Thus 


dD jb, " 
G(6, 6) = 2 ae - (26) 


= ot), = =| _.8in 000 


0 


The denominator of this expression may be transformed 
by the aid of the following orthogonality relations for 
the Legendre polynomials? 


is dP,(cos 6) dPm(cos 6) 
ih ——_—_ ———— sin 60 
0 d6 


d0 
0 for n~Am 
= |2n(n + 1) 
2n +-.1 


Then, performing the integration of (26) and introduc- 
ing 6, from (18), the antenna gain is reduced to 


IRE. TRANSACTIONS ON ANTENNAS AND PROPAGATION 


dS jrt(2n + 1)jn(hro) Rn 


n=0 


dé 


G(6, ¢) = (28) 


> n(n + 1)(20 + 1)jn2(Rro) | Rn |? 


n=0 
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; 


: 


This expression has been calculated numerically; the re- 


sults are presented in Fig. 2. The computations have 


been made for a shell of polystyrene (k.=2.6) with an | 


electrical thickness 0=+/k-kd of 7/16 and 7/4, these 
thicknesses corresponding to the values 0.195 and 0.779, 
respectively, of (x.—1)kd. The numerical results, how- 


ever, may be used for an electrically thin shell of any ke, — 
the value of the susceptance B= (x,—1)kd to be one of © 


those specified above. For comparison the antenna gain 


$ sin? 6, when the shell is removed, is also shown. The — 


radiation pattern is shown for a shell, the radius of 


which is given by k71=9 and kr: =12, 2.e., approximately © 


1.5 and 2 wavelengths, respectively. Two values 3 and 6 
of the electrical distance kro from the dipole to the 


center of the shell have been considered. From the pat-_ 


terns the disturbing influence of the thin shell on the 
radiation from the dipole is seen to be the more pro- 
nounced, the greater the susceptance B = (x.—1)kd of the 
shell is. Thus, in order to obtain a small distortion of the 
pattern for any 7 and 7, the chosen susceptance 
(ke—1)kd should never exceed 0.2. The dipole is seen to 
radiate a relatively large amount of the transmitted 
energy in directions close to 6 = 180°. This type of distor- 
tion is seen to be stronger, if the susceptance of the shell 
and the distance of the dipole from the center of the 
shell are greater. To avoid strong distortion of the pat- 
tern, 7o should not be chosen greater than one wave- 
length (kro~6). The radius 7; of the shell is seen to in- 


fluence the pattern much less than 79. Results that are 


not to be discussed here, have been obtained for a shell 
with a radius of approximately one wavelength and 
show that the shell distorts the pattern only slightly. 
This is to be expected according to Rayleigh’s law of 
scattering. 


ANTENNA IMPEDANCE 


The spherical shell reacts upon the surrounded an- 
tenna, a wave being reflected back along the antenna 


line, thus changing the antenna impedance. To deter-_ 


mine the resistance of the dipole inside the spherical 


shell the Poynting vector method may be used. To de- . 


termine the reactance also, the emf method may be 
used. According to the former method the radiation re- 
sistance R, is expressed by 

Pr 


Rigi os | 
ae 


The total power transmitted by the dipole is found from. 


the preceding section to amount to 


Pr 


rae n=0 ae +1 (30) 


(29) 
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Introducing (18) for 5, and (13) for C, the radiation re- 
sistance of the dipole normalized with respect to the 


self-resistance Ru of the dipole 


67 
is found from (29) to be expressible by 
R ore in(Rro) |? 
ee nin + 2)(2n 4-1) | Ra? | ay) 
Ru 2 n=0 kro 


Next, let us find the dipole impedance by the emf 
method,’ according to which, assuming the dipole to be 


7R. E. Burgess, “Aerial characteristics,” Wireless Eng., vol. 21, 
pp. 154-160; April, 1944. ; 
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Fig. 2—Primary radiation pattern of a radial dipole surrounded by a dielectric spherical shell for different values of the electrical radius kr: 
and susceptance B=(xe—1)kd of the shell and the electrical distance kro from the dipole to the center of the shell. 


matched with the shell removed (it is no loss of gen- 
erality to assume that) 
E,Al 


Li — Ry -— —> (33) 
I 

I being the prescribed current of the dipole and £, the 

radial electric component of the secondary field in the 

position of the dipole inside the shell, which component 

is found from (8) and (14) to be 

jn(kro) 


kro 


ih > a,n(m-+ 1) (34) 


JHO n=0 


use being made of the differential equation (21). We in- 
troduce this into (33) and insert at the same time (19) 


E, = 
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Fig. 3—Relative resistance and increase of reactance of a radial dipole surrounded by a dielectric spherical shell as a function of the electrical 
distance kro from the dipole to the center of the shell, for («e—1)kd =0.779 and for different values of the electrical radius kr; of the shell. 


for a, and (13) for the dipole constant C. Thus, for the 
dipole impedance normalized with respect to the self- 
resistance (31) 

ZA 


3 2) 
Ru 1— ror >> n(n + 1)(2n + 1)Rn 


n=0 


oe | d[krihn® (kr1) | 


kro d(kry) J - Daa oy 


The real part of this expression may be shown to be 
identical to (32) by introducing into the above expres- 
sion the following formula obtained from (32) by setting 
R,=1, which corresponds to the case where there is no 
shell 


3 < jn(Rro) 3 
~ Donn Hn + 0] * ee: (36) 


0 


n=0 


At the same time the antenna reactance (actually the 
reactance increase) may be shown expressible by 


Xa 
ae) 


n=0 


eset 


A | [t. aie, 


d(kr1) 


where 7, is the spherical Neumann function.’ R4/Ru 
and X4/Ri have been made the object of numerical 
computations, the results of which are presented in 
Fig. 3. The computations have been made for a shell, 
the electrical thickness of which is determined by 
(ke—1)kd =0.779. The thinner the shell, the nearer the 
relative radiation resistance and increase of reactance of 
the dipole are to 1 and 0, respectively. The shell is seen 
to influence the dipole impedance strongly as the dipole 


3 co 
= — >) n(n + 1)(2n + 1)| R,|? 


is placed near the center of the spherical shell. By plac- 
ing the dipole more than half a wavelength from the 
center of the shell but not too close to the shell, the 
vswr on the line feeding the dipole is seen to be less than 
1.2, while a vswr of about 2.2 is obtained by placing the 
dipole at the center of the shell. By the way, only the 
term 2 =1 of (32) and (37) contributes to the impedance 
of the centered dipole, the dipole resistance thus being 
Ra/Ru = R,| 2. 


CONCLUSION 


The numerical results obtained for the gain and im- 
pedance of a radial dipole surrounded by a dielectric 
spherical shell seem to show that a compromise is nec- 
essary as to the position of the dipole inside the shell in 
order to obtain a small influence of the shell on the gain 
as well as on the impedance of the dipole. The numerical 


results are also expected to be practically applicable to | 


a dipole of finite length. 


Ec - 1)kd 
d[krijn(kri)| dl krin,(kr1) | i 
d(kr:) d(krn) | 
d[krijn(kri)| dlkrinn(kri) | ; 
rename rae on 
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The Current Distribution and Input Impedance 
of Cylindrical Antennas* 


E. V. BOHNt 
Summary—This paper treats both the formulation of the problem ee : : : : 
and its solution in a fundamental way permitting only approximations T’=unit vector in the direction of current flow 
whi-h appear evident. A system of four integro-differential equations U=scalar pot. 
is derived and these are applied to the cylindrical antenna problem. r? = (x—x')?+ (y—y’)?+(g-2’)? 
— 


An integral equation of Hallen’s type occurs and this is solved exactly 
by transform methods. The input impedance of a full wave antenna Be ao Bee 
is computed numerically to show that the present theory is consistent that H=i1T>N 
_ with other results. 


N=unit vector normal to the surface and directed so 


qg=surface charge density 


INTRODUCTION gmigumg 
HE cylindrical antenna problem is fundamental re ON ay” eR Fijeew ae! Minar bisel ee SSB 
in antenna theory and considerable effort has : ; 
integration xX, y, 2 


been devoted to its solution. For moderately 
thick antennas the integral equation approach is the 0 d d : : 
most suitable, but it has not been fully exploited due Mi cern ay!” wig with respect (ota 
to the mathematical difficulties involved. Hallen’s! 
method of iteration gives asymptotic results only and 
is dependent on an arbitrary choice of expansion param- E;=tangential electric field of the generator. 
eter. It restricts itself to totically thi t = fa : : 
; ire Baebes! grviees coud te een The £ and A fields are determined by the following: 
with a discontinuous ring driving source and is incon- 


pont <7 ) so =u 


‘sistent in that it obtains a finite input impedance. =n isr> 

It does not give a clear distinction between the cylindri- H(p’) = a f as (1) 
cal antenna and the cylindrical cavity and it is incapa- 
ble of computing the free modes of oscillation of the Elo’) st ee VU. (2) 


antenna. 
In order to obtain a solution which satisfactorily To arrive at the fundamental antenna equations, we 


overcomes these objections, a more fundamental ap- let p’ approach a point on the surface. Taking into 
‘proach must be made to the problem. Infeld? has shown account the singularities in the integrals and the bound- 
the importance of the gap in determining input imped- ary conditions, it is a simple matter to derive the follow- 
ance. However, his method is not suited to the cylindri- ing: 


cal antenna. Synge® has formulated an exact integral pete e fide rs 
equation for the problem which is then solved by itera- tx No) =— f 1X Vwds (1a) 
tive methods. This is basically Hallen’s original method da 
of splitting the kernel into a singular and a finite part >> 4g 5 
‘and has been sufficiently criticized in the past to indi- a Se (2a) 
cate the necessity of a more fundamental approach. 
‘In the interest of simplicity we restrict ourselves to CEP eS E; (3) 
perfectly conducting closed surfaces and introduce the 
following notation: where 
ah 
es Peper sed (oS. ie, if + Was + ee (p’) 
a 4dr J & 2 & 
r —s 
i=surface current density div i = — jug. (4) 
A = vector pot. This is a system of four integro-differential equations 


for the four unknows, 7, g, and the two surface com- 


ae . . . . 

* Manuscript received by the PGAP, July 9, 1956. ponents of T. The antenna problem distinguishes itself 

+ Canadian Aviation Electronics, Montreal, P. Q., Canada. : : . Bae : 

1 For a thorough discussion of this problem and literature refer- from the cavity problem in the choice of N which is 
ences, see S. A. Schelkunoff, “Advanced Antenna Theory.” => 

* L. Infeld, “The influence of the et a ee gee upon the theory directed toward the field. The antenna N is directed 

f ant ,” Quart. Appl. Math., vol. 5; July, 1947. 3 , ‘ 

i. G. BAibere and on Synge, “The general problem of antenna outward. A general discussion of this system appears 
radiation and the fundamental integral equation with application to oyt of the question, but if we restrict ourselves to the 
an antenna of revolution,” Quart. Appl. Math., pt. 1, vol. 6. See also, Gindacal Be tie cunt one ae considerabitteime 
J. L. Synge, pt. 2. , 
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plified. We consider a completely closed cylindrical 
antenna (Fig. 1). Instead of applying the group [(1a), 
(2a), (3), and (4)] to this case it is simpler to start with 
a slightly modified but equivalent group of equations. 
For the cylindrical portion of the antenna we have 


dE, 9E, 


(1b) 
Op Oz 
SL, to (2b) 
di 
— = — jug. 4a 
dz a8 ee) 


The first of these equations is equivalent to (1a) in that 
—~ 


it is obtained from the boundary conditions for H. We 
set 


2rat = I 2raq = Q 
and obtain 
ad’I OF, 
eas BL ate ae 2najubo~— (5) 
2rjw&o ('( OE. 
I(z) = Cy, sin B( — z + 6) + ——— ed Ce 
B z dp 
-sin B(z — u)du. (6) 
Eq. (3) is | 
1 ) 
~ Fy = —jo™ f ipds + f ee 
4rJ 4 4no/ ssc O02 


A refers to the cylindrical portion of the antenna and 
C to the end caps. We now restrict ourselves to a center- 
fed antenna so that Q(+/) =—Q(-—1) and assume the 
charge to be continuous. Partial integration and use of 


(4a) and (5) gives us 
gel Ae 


OW 
1 Qr 
asf (iia iat Cicely ae 


Pot —*) yaeds (7) 


foes 


where 


(ult al) — vl t+ al da 


In F, we have taken g=q(+/) over the cap and also p 
to be constant for a fixed value of wu. 
We now set 


dE, ’ ” 
ae waeni la ae one! 
Op 
and consider (7) to be an integral equation for the y’s. 
y’(z) has an impulse type singularity at z= +]; 
y’’(z) has a doublet type singularity at z= +1. 
These functions are defined by 


+1 Q 
Jf s'@pe- nas = 5 {fo - w) + (-1- w} 
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Fig. 1. 


+1 
yl (2) f(z — u)dz 
7 A Qa fe) 4 
Fics 4r&o = ou ft - 


— f(-1— w)}. 
To determine y we have the integral 
V a Qr +1 
bw) =— gt =—f fy |s— ul )asdd ©) 
2d 4rJo ey 


where 
= (z — u)? + 2a7?(1 — cos $) 


and V is the “driving voltage.” From (6) we have 


wQ 
I(z) = Ci sin B(l — z+ @) afueds Fs B(l — 2) 


— jo cos B(1 — 2) 


27jw&o i 
+ ol. y(u) sin B(z — u)du. 


Over the end caps C of the antenna we have 
1 9 ; 
+= "(phiG)) == = Fo0Gry: 
p Op 


We set [i(p) = — Bp and find B =jwQ/2. To determine 
C; and 6 we have 


1a) ol 
+ = 1) =a; 5 Jem +i = ju. 


Neglecting higher powers in a, this leads to 


>; 

Cy a ao } = os 

B 2 

I(z) = Ising (1 — 2+ <) 
2rjw&o (' ’ 
ag am y(u) sin B(z — u)du (2>0) (9) 
where 
jeQ(+D 


In = 
B 


Eq. (9) is correct provided yy contains no further 
singularities. The problem is now to solve (8) for y and 
then to obtain J from (9). The manner in which Jy is 


on — —— 
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determined will be discussed later. We now turn to the 

solution of the integral (8). To obtain a solution we 
begin with the following relation:4 

e—7Br 1 

on jr 

r=(g—u)?+ oP 

p? = 2a7(1— cos ¢). 


[ KeleV=BF= PerP as 
Cc 


We substitute S=S/] and obtain 
e#Br 1 


= f Ko(4 v= = BR) ern 
r jrl Jo l 


We have® 


1 Qr p we Ss 
“4h Ko (4 Vat = )d6 
0 


a La es a = 
= Kp (< / — Bz — =) Ty (é. /— BP — s) 
= Kolo|s|. 


Ij(z) represents a Bessel function and should not be 
confused with the current amplitude Jp. 
Hence 


24 


1 1 
= wd = Gal K olo|s|e#!#“-9)ds 


Qa 0 vhs eC 
b= 7) 
= De peter VE a0) 
n,m p=jrn 
where 
1 sinh (s — 
R(t, p) = Ram = cep Beles] me 
; jr Jo s— >? 
sinh —t 
Ck eee (11) 


Sat 
The integral equation is now reduced to a system of 
linear equations: 


V 
Fes 4 ees LRnm im 12 
ae pe y (12) 


where 
1 +1 


== ye imag, 
21/1 


Vm => 
Fig. 2 illustrates the path c and the two branch cuts 
W;, and W?% to the branch points +j6/. To evaluate (11), 
we take p and ¢ to the left of c and break up the inte- 
grand into two parts. One part can be evaluated by 
deforming c into W2; the second part can be evaluated 
by deforming ¢ into W; and two integrals around the 
poles p and ¢ which can be readily evaluated. The 
integral along the branch cut Wi can then be trans- 
formed into an integral along W2. The result of this 
operation is the following expression. 


4 See S, A. Schelkunoff, “Electromagnetic Waves,” p. 414. 
5 See Watson, “Theory of Bessel Functions,” p. 360. 
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Fig. 2. 

1 e*stptt — cosh (p — 2) 
lk t, (i — A Jo 0 ——_ 
se 5, font} (s— p(s) 


e-?s—p—t — cosh (p — , 
GEpepie jo 
i sinh (p — #) 


Nipemarrels cies etl: 


For Joly we have the expansion 


Te /fan= 
Jolofs] = 1 -— (2) P+ 8) +o. 
EVERY! 

Hence we can set Jolo[s]~1 in the integrand for the 
lower order matrix elements. The integrals can then 
be reduced to the familiar sin-cos integrals by suitably 
choosing W2. The results are given in the Appendix. 

Since 


mee p=t 
as ai re 


we have the result that the lower order of diagonal ele- 
ments are nearly independent of a and the diagonal ele- 
ments have a logarithmic singularity in a. Hence asymp- 
totically as a—0 the diagonal elements become very 
large and this leads to the following approximate solu- 
tion of (12) 


TRnn 


V 
Ran Wie Te. ia ml: 
bai i! oe Ths § } 


The infinite system of linear (12) reflects the proper- 
ties of the integral (8). We can now obtain a clearer 
insight into Hallen’s method of iteration. Hallen’s 
method in effect attempts to compute all the matrix 
elements by means of (12a). Only the lower order diag- 
onal elements are dominant due to the logarithmic 
singularity of Ko at the origin. Hence we can apply the 
iterative solution (12a) to the lower order elements 
only. This is the fundamental reason for the asymptotic 
divergent character of Hallen’s solution. 

Having solved (8) we can substitute y into (9) and 
obtain for the current the following expression 


(12a) 


I(z) = Io sina(? —2Z +<) + V-fi(z). (13) 
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fi is now a known function but the relationship be- 


tween J and V is still to be determined. To obtain this, 
we can derive the following result 


1 Qr +1 
-—f f x@w|s = a] )aeag 


r? = (g — uw)? + a? + p? — 2apcos@ 


Ek, = 


where 


in the same manner as we did (8). 
For «=0 and p— = we must have 


= ih Nina 


Ez 
= = 120, = ——__—_——__ - (14) 
A. years 
— | I(z)dz 
2rd 0 
This leads to 
Al 21 240 _ 
Geo 1 


THE CASE OF THE SINUSOIDAL CURRENT 
DISTRIBUTION 


If we imagine V fixed and let a—0, we see from (12a) 
and (14a) that, in general, J is of the same order as fi. 
- Hence as a0, fi(z)—0, Zo>—0, but the sinusoidal term 
does not necessarily become dominant in the limit. 
Since the sinusoidal current distribution has played such 
a successful role in antenna theory, it is of interest to 
see under what circumstances J» will become dominant 
over fi. As a—0, we can solve (12) asymptotically and 
obtain 


1K (84) Yn = us Shah oie pela 
2d LK (Ba) 
Hence 
rays Iusin (1 Sa +)+ ab NAN ti 
D 1201K (Ba) 
‘sin B(d — 32) z>0 
where we have taken 
1 | a | <d 
gm = 1 Rie d<l 
and 
sin B(d — z) = 0, Shy), 
Since 


I) and —————— 
1201K (Ba) 


are of the same order, the condition that the sinusoidal 
term dominates is that d—0. The sufficient conditions 
for a sinusoidal current distribution are that a—0 and 


that d—0. Actual practical antennas approximately 
meet these conditions. 
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From the theory of equations of type (8) and (12) it i 
can be seen that we must choose a fixed d>0. If we at- 


tempt to choose d=0, (12) becomes invalid since 


Dl esl2= co. The input current is infinite under this — 


condition. 
The sinusoidal current distribution is also very suc- 


cessful in the computation of antenna impedances. This — 


is all the more surprising because of the very great dis- 


crepancies in the tangential electric fields of the sinus- 


oidal and the actual current distributions.-To obtain 
a more fundamental understanding of this situation, we 
derive the following: 


— => = ad 
pa — pr = jono f | He — H|tar — jobs f | E, — Ei |?dr 


~ = > > —-> 
+nf |B, — Bi ltds + [By Git — is*)ar. (13) 


p is the complex power of the electromagnetic field and 
the indices refer to two different field distributions. Eq. 
(15) can be easily proven by the use of the reciprocity 
theorem and the usual vector identities. 

We can now state the following minimum theorem: 


=> 
if 7, is the correct current distribution meeting the 
—_ 
boundary conditions and if z2 is a second current distri- 
=> 


bution such that 7;=72; over the generator region, then 
the correct current distribution requires the least real 
power. 

The condition that the input current be held constant 
is approximately fulfilled in the case of a concentrated 
generator and provided that sin BJ~0. This minimum 
theorem is the reason why the sinusoidal current distri- 


bution can be successfully applied in computing the 


radiation resistance. 


FREE MODES OF OSCILLATION 


These are defined by setting g, =0 in (12). The eigen- 
values are obtained from ||Rnm(8)||=0. Due to the 
asymptotic behavior of the diagonal matrix elements, 
the problem of solving for the eigenvalues may not be 
so formidable as it first appears. This would be of im- 
portance in the study of the transient behavior of an- 
tennas. 


NUMERICAL RESULTS FOR A FULL WAVE ANTENNA 
We take 


(eg — = 37.1 —=5 
2a a 
1 |u| <d 
ew) = 4 
O 5 a dered (16) 


The arbitrary choice of d brings out a fundamental — 
weakness in present antenna theory. E; is determined — 


by the terminal driving conditions from a coax cable 


or transmission line and can not be arbitrarily chosen. — 
The only known property of E; is that it is “concen-_ 
trated” near the driving point. Without further infor- | 
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mation about £; there is no purpose in taking all the 
Fourier coefficients of (16) in the calculation. We there- 
fore compute y, from g,(m=0, 1,---, 5) and set 
¥n=O(n=6, 7,---,). This simplifies the numerical 
work and still leads to a field which is “concentrated.” 
From (12a) we now obtain 


21 

op 30 = 0:2615 + j0.1458 
21 

Fp Nt = 0.2870 + j0.0679 
21 

— yz = 0.3812 — 70.0098 
V 

21 

Fy 98 = 0.4050 + 70.0070 
21 

— 4 = 0.4573 — 70.0052 
V 

21 

Fp 38 = 0.4698 — 0.0036. 


For the function f; in (13) we have 


fe) t 0.0433 + 0.0773) cos —— ~ (= ) 
ieee awa oe Py ag Ne 
=\( =) Jame 
— 1 — —) sin — 
aie 
THES 


ete 7 2b a l 
FFs) n? —1 
From (14a) we get 
Io 
240 a = 0.2237 + j0.2360. 
For the complex power we have 


I 
p= if E,I*dz = 
0 


Y; is the input admittance and Zp is the loop imped- 


LY, | Vi? = 47, | To |?. 


ance based on the sinusoidal current. We find 


1 
27309510; 


t 


Z; = Zo = 464 — 7362. 


The loop impedance calculated from the field of a 


‘sinusoidal current alone is Z) =200+ 120. 


Since the input current is not even approximately 


‘constant, we can not expect the minimum theorem to 


be applicable here. From Schelkunoff* we have 


Z; = 350 — 7660. 


6S, A. Schelkunoff, “Antenna Theory and Practice,” p. 439. 
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The result (17) is seen to be consistent with other re- 
sults. An accurate comparison between the various 
theories and experimental results is not possible due to 
the lack of knowledge concerning E;. 


CONCLUSION 


The method of solving the cylindrical antenna prob- 
lem by the use of Fourier coefficients is superior to 
Hallen’s iterative procedure both theoretically and 
numerically. Improved solutions by Hallen’s method 
would require graphical analysis, while the solution of 
(12) by arithmetic means involves no theoretical diffi- 
culty. The one objection to (12)—that a large number of 
terms may be required—cannot be answered until there 
is more detailed knowledge of E;. Asymptotically as 
a—0 (12) can be just as readily solved for an impulse 
type of solution (requiring an infinite number of terms) 
as the integral (8) directly. The complex power and the 
input impedance can be computed exactly without 
graphical analysis as it would have to be done by Hal- 
len’s method. The various theories give comparable 
results because each uses a “concentrated” type of field. 


APPENDIX 
Matrix ELEMENTS FoR 1/27fo2"pdo 


Ram is symmetrical about both diagonals 
ee ae Ween 
Reon = Ran 


We use the short hand notation 


(n ~m). 


Kolo[s] = Ko (= V=BE= 3) ip (= V=FP= 3) 


where Jp and Ko are the modified Bessel functions of 
the first and second kind, respectively. 


am > Bl 
rm > Bl 
Anj(m — m)(—1)"*"lRnm 
sty (xn — Bl)(rm + Bl) 
(am — Bl)(arn + Bl) 
+ Ci 2(rn + Bl) — Ci 2(am + Bl) 
+ j| —Si 2 (an — Bl) + Si 2(am — Bl) 

— Si 2(an + Bl) + Sil2(rm + Bl)} 
am < Bl 
am > Bl 

Arj(n — m)(—1)"*™"lRam 
nit (61 — wn)(mm + Bl) _ 
(am — 1) (an + Bl) 
+ Ci 2(6l + xn) — Ci 2(61 + 1m) 
+ j{Si2(6l — rn) + Si 2(rm — Bl) 
— Si 2(61 + xn) + Si 2(61 + wm)} 
an > Bl (20) 


nAzAm (18) 


— Ci2(rn — Bl) + Ci 2(am — Bl) 


(19) 


nA~m 


Ci 2(61 — an) + Ci 2(am — Bl) 


(n = m) 
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el 1¢ t 1 
deste Kelulixn a oc ee 
a oleljen| + 4j ta ny Gh eee 


1 
Js - Ci 2(an + Bl) + isi Ci 2(an — Bl) 


+ — {Si 2(rn — Bl) — Si 2(nn + Bi)} 


wv |S. 


mm < Bl (21) 


e?7Bl — 1 1 
lRan =o K v8 j Sites Sage \ a Oe 2 eos ae + ==} 
ofoljxn] + 4j . a oe 


(n = m) 


1 1 
ae % Ci 2(4rn + Bl) + ag Ci 2(@1 — an) 


is = {x — Si.2(6l 4h) — Si 2(81 — xn)}. 


For Bl=7 and n, m= +1 we have 


Hie t (22) 
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Anj(n — 1)(—1)"Uka 


== Ry erased C+ In2 — Ci2a(m — 1) + Ci2a(n+ 1) 
(n + 1) 
— Cid + j{—Si2a(n — 1) — Si 2e(n + 1) + Si 4c} 
m>1 (23) 
—Anj(m + 1)(—1)" Vk_1m 
= op ED see le alal, Ea ote ee 
(m — 1) ~ 


— Ci2a(m + 1) + jf —24 + Si 4m + Si 2x(m — 1) 
+ Si 2x(m + 1)} 


1 a 3 1 
EROS Sg oa ra pubs 3 
a fe Chaar es St te (24) 
2 2 
Anjlky, = In2a +C+1n2—Ci4r+j7Si4r (25) 
—4rjlky,) = In2a + C+ 1n2 — Ci4e 
+ j[—a + Si 4a. (26) 


Scanning Lens Design for Minimum Mean-Square Phase Error* 
E. K. PROCTOR} anp M. H. REESt 


Summary—Microwave lenses of the constraining or binormal 
type having superior wide-angle scanning properties can be achieved 
by suitable choices of lens contours and of the manner in which the 
index of refraction varies from point to point across the aperture 
plane. The phase error (i.e., departure from a plane phase front at the 
exit surface of the lens) can be expressed in terms of three unknown 
functions: the equations of the two lens surfaces and of the index of 
refraction. The square of this error is integrated over the desired 
angle of scan and also over the aperture area of the lens. The three 
unknown functions which minimize the phase-error integral are then 
determined from the Euler-Lagrange equations of the calculus of 
variations. Graphs of phase errors for several typical lens designs are 
included. 

INTRODUCTION 


Merrett lenses in which metallic plates 


are used in such a manner as to constrain the 
internal wave propagation to directions parallel 
to the lens axis have been called “constraining” or “bi- 
normal” lenses. Snell’s Law is not obeyed at the surfaces 


* Manuscript received by the PGAP, August 22, 1956. This paper 
(less the Appendixes) was presented at the WESCON Convention, 
Los Angeles, Calif., August 21-24, 1956. Work based on a study pro- 
gram carried out by Sperry Gyroscope Co., under Contract DA 
36-039-sc-15323 for the Evans Signal Lab., of the U. S. Signal Corps. 

t General Electric Microwave Lab., Palo Alto, Calif. Formerly 
with Sperry Gyroscope Co., Great Neck, N. Y. 

{ University of Colorado, Boulder, Colo. Formerly with Sperry 
Gyroscope Co., Great Neck, N.Y. 


of such lenses and the refractive properties are de- 
scribed most conveniently in terms of the differential 
phase shifts imparted along various ray paths through 
the lens. Since these phase shifts are independent of the 
angles of incidence (neglecting surface-impedance ef- 
fects), constraining lenses have inherent advantages in 
scanning antenna systems. 

Ruze’ in 1950, and Ellis, Fine, and Reynolds? in 
1951, have reported analytical and experimental work 
on such lenses. A brief summary of the Ruze meth- 
od of analysis is in order here by way of introduction 
to the analytical method of this paper. 

The geometry and notation for what will be called a 
“two-dimensional” lens are shown in Fig. 1. It is as- 
sumed that the lens is so constructed that ray paths 
within the lens are parallel to the x axis for all incident 
and emergent rays. The lens is illuminated by a suitable 
point source which can move along a circular feed path 
to cause the emergent beam to scan through a total 
angle of 26,°. The feed circle is centered at the origin 


1 J. Ruze, “Wide-angle metal-plate optics,” Proc. IRE, vol. 38, 
pp. 53-59; January, 1950. 

_? Ellis, Fine, and Reynolds, “A Point Source Binormal Lens,” 
Air Force Cambridge Res. Labs., Rep. No. E5067; March, 1951. 


BEAM 
DIRECTION 


Fig. 1—Schematic diagram for scanning lens. 


and has a radius equal to the focal length, f. The aper- 
ture or width of the lens is 2; and its center thickness 
is Dy. The lens is to be designed to convert the incident 
circular phase front from the point-source feed into a 
plane emergent phase front. The equations of the illumi- 
nated surface, x(y), the lens thickness, D(y), and the 
index of refraction for any ray path through the lens, 
n(y), are to be determined. These three equations are 
‘independent and may be considered as design param- 
eters which are to be so chosen as to minimize the 
scanning aberrations of the lens. 

_ Since three independent equations are to be deter- 
mined, three independent conditions must be specified. 
In Ruze’s method of analysis, two conditions are pre- 
scribed by the requirement that perfect focus be ob- 
tained at two symmetrical off-axis feed positions. This 
leads to the explicit requirement that the illuminated 
surface, x(y), be an ellipse and yields an equation relat- 
ing n(y) and D(y). For other feed positions, Ruze ob- 
tains an expression for the phase or path-length error, 
5, and expands it in a power series in the coordinate y. 
(The power series converges for all feed angles if the focal 
length-to-aperture ratio is greater than 0.8.) In the 
power series the coefficients, which are functions of the 
feed position, 0, and of the lens equations, represent dis- 
tortions of the phase front. The effects on these coeffi- 
cients of several choices of the third condition including 
specification of a third perfect-focus feed position (the 
‘triple-correction-point design), are considered in some 
detail in Ruze’s paper. 

A more general method of determining the three un- 
known functions is to choose them so that the phase 
‘error is minimized throughout the scan angle, 261, and 
over the entire lens aperture, 21. This can be accom- 
plished by using variational methods. 


THE VARIATIONAL SOLUTION 


Consider two rays originating at an arbitrary point, 
P, on the feed circle as shown in Fig. 1. One of these is a 
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reference ray passing through the center of the lens. 
The other passes through the lens at an arbitrary dis- 
tance, y, from the origin. The path-length difference 
is given by. 

b= f+ mDo + s — f’ — ny) Dy) (1) 
in which m) and Do are the values of u(y) and D(y) 
at y=0 and the other quantities are indicated in the 
figure. For a perfect plane phase front inclined at an 
angle 6 with respect to the y axis, 6 is zero. Hence, it is 
desirable to make 6 as small as possible for all feed posi- 
tions and for all ray paths. 

Using the relations 


f’ = [(fcos@ + x)? + (f sin @ + y)?]1? (2) 
and 
s = (x + D— Do) cos0+ ysin 8, (3) 
the error term may be expressed as 
5=f— [(fcos@+ x)? + (f sin 6 + y)?]1? 
+ (mDo — nD) + (x — D— Do) cos@+ysind (4) 
or 
(4a) 


The square of 6 may be integrated over the lens aperture 
dimension, 24;, and over the full scan angle, 26;. Thus, 


yi 
|| 


—Y¥L 


o= 5[x(y), D(y); ny), 9, y, constants]. 


01 
iH 5?[x(y), D(y), n(y), 9, y, constants]dédy (5) 
= 6, 


is an integral expression which is a measure of the cumu- 
lative phase error. In particular, since it comes from a 
definite integral, A is a nwmber which has a specific value 
for each combination of x(y), D(y), and n(y). If three 
explicit functions can be found such that A is minimized, 
they determine a scanning lens which is optimum in a 
mean-square sense. 

Since the functional dependence of 6? on 6 is known 
explicitly from (4), the integration with respect to @ 
in (5) can be carried out to yield 


a= [ Gl), DO), mG), », olay (6) 


=i 


in which 
(6a) 


The explicit form of G is given in Appendix I. Solution of 
(6) for the minimizing functions x(y), D(y), and n(y) is 
a relatively simple problem in the calculus of variations 
since the derivatives of the unknown functions are not 
involved. 

[Though the details are more complicated, the logic 
of the variational method is the same as for ordinary 
maximum and minimum problems in calculus where the 
expression to be minimized is differentiated once with 
respect to each of the independent variables and the 
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derivatives set equal to zero. In the present case, the 
variables are the three unknown functions x(y), D(y), 
n(y) and the procedure is roughly equivalent to setting 
equal to zero the partial derivatives 


LYN oA EAN 
ARR Bs : 
dx(y) an(y) 


) 
dD(y) 

This results in three simultaneous equations from which 
the unknown functions may be determined. It should 
be noted that these equations are only necessary condi- 
tions for the existence of an extremum which may be a 
maximum, a minimum, or a point of inflection. The 
evaluation of sufficiency criteria for the existence of a 
minimum is a problem of considerable mathematical 
complexity. In most engineering problems it is usually 
advisable to rely upon geometric or physical intuition 
to interpret the results. | 

The variational problem represented by (6) is of the 
type involving several unknown functions of a single 
independent variable, y. It can be shown that the three 
equations which must be satisfied are 


Diaz fer a 0, (7a) 
dD(y) 
eee of (0) (7b) 
dn(y) 
and 
pelle (6 (7c) 
dx(y) 


(These equations are derived using somewhat different 
notations in the works of Courant and Hibert® Margenau 
and Murphy,* and Weinstock® and several other texts 
which include chapters on the calculus of variations.) 
Performing the indicated operations (using the expres- 
sion for G from Appendix I) yields, after some manipu- 
lation, 


N — Lsin 61/01 


x — D— Do => 5 ) (8a) 
, 2 sin? 0; 
6, + sin 6; cos 8, — ————— 
1 
AD, (« + D — Dp) sin 6, 
noDy — nD = = — f» (8b) 
264 0, 
(NV — [sin 6,/61) (Zi sin 6;/04 aa N’) BIL 
c 2 sin? 6; 26; 
6; + sin 4; cos 6, — Aleta 
1 


+ 2f sin 0; — yM’ + 226, = 0. (8c) 


3 R. Courant and D. Hilbert, “Methods of Mathematical Phys- 
ics,” vol. I, Interscience Publishers, New York, N. Y., chap. IV; 1953. 

4 H. Margenau and G. M. Murphy, “The Mathematics of Physics 
and Chemistry,” D. Van Nostrand Co., Inc., New York, N. Y., 
chap. 6; 1945. 

® Weinstock, “Calculus of Variations,” McGraw-Hill Book Co., 
Inc., New York, N. Y.; 1952. 
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The expressions L, M, N, and 


OL oN oN 
Ls » M’=—, and NM =— 
Ox Ox Ox 


4 
| 
; 


contain elliptic integrals and are given in Appendixes 
I and II. Eq. (8c) contains x as the only unknown func- 
tion, but in a very involved implicit form since x appears 
in the arguments of the elliptic integrals. However, a 
graphical solution is possible for specific values of 0: and 
the other constants. To avoid double interpolation in the 
elliptic integral tables, it is advisable to change the vari- 
ables to 


a = tan es (9a) 
x 
and 
Ax/ aaah 
K = sin? / Nee S) i ( 
PAR a? he eG 


These relations are illustrated in Fig. 6. Note that the 
positive x and y axes are to the right and upward 
respectively so that x(y) has only negative values. It is 
also desirable to normalize the dimensions to the focal 
length, f. The required transformation is carried out in 
Appendix ITT. 

The values of x for a given range of the independent 
variable, y, may be substituted into (8a) to obtain the 
corresponding values of D(y). Since D(y) is given ex- 
plicitly, graphical solution is not required. It should be 
noted, however, that the two quantities in the numerator 
of (8a) differ only in the fourth or fifth decimal place. 
Thus, an eight- or ten-place table of incomplete elliptic 
integrals is required.® The values of x(y) and D(y) from 
(8c) and (8a) enable one to find the desired values of 
n(y) quite readily from (8b). 


NUMERICAL EVALUATION 


Eqs. (8a) through (8c) are not very amenable to 
interpretation except on a numerical basis. Accordingly, 
calculations of lens contours, index of refraction, and 
phase error were performed for two versions of the LMSE 
(lease mean-square error) lens. Corresponding calcula- 
tions were carried out for lenses whose design equations 
are determined by conventional methods in order to 
have a basis for comparison. 

Consider first a lens having a plane illuminated sur- 
face, x(y)=0. This will be called a two-parameter lens 
since only n(y) and D(y) can be varied to optimize the 
phase error. Such a lens is not especially good for scan- 
ning, but it offers the advantage that setting «(y)=0 in 


6 “Table of Complete and Incomplete Elliptic Integrals,” reissued 
from Legendre’s “Traité des Fonctions Elliptiques,” Paris, 1825; 
issued by the Biometrika Office, University College, London; printed 
at the Cambridge University Press, 1934. 


--—PLANE ILLUMINATED FACE LMSE 
——-— --SINGLE- POINT- FOCUS, CONSTANT n(=0.6)- 
SINGLE-POINT FOCUS, CONSTANT THICKNESS 


Fig. 2—Phase-error data for two-parameter lens designs. 


e 


--—PLANE ILLUMINATED FACE LMSE 
=—---- SINGLE-POINT-FOCUS, CONSTANT n(=0.6) 
SINGLE -POINT FOCUS, CONSTANT THICKNESS 


Fig. 3—Phase-error data for two-parameter lens designs. 


Fig. 4— Phase error of triple-correction-point lens. (b) phase 
ay =a “a error of LMSE lens. 
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the equations eliminates (8c) and considerably simpli- 
fies (8a) and (8b). For comparison purposes consider 
two other lenses having plane illuminated surfaces and 
designed by the usual methods to give perfect focus at 
#6=0. One of these has a constant thickness, D(y) =Dy, 
and a varying index of refraction. The other has a con- 
stant index of refraction, n(y)=0.6, and a varying 
thickness. For a scan angle of 26,=60° for the LMSE 
lens, quite reasonable ranges of n(y) and D(y) are ob- 
tained. The phase errors of the three lenses for several 
feed positions are shown in Fig. 2 and Fig. 3. For con- 
venience, both 6 and y are normalized to the focal length, 
f. Values of 6/f are shown only for positive values of 0 


_since the lenses are symmetrical about the x axis. Except 


for the negligible error at @ =0, it can be seen that for y>0 
the LMSE design gives the smallest error. For y<0 the 
error of the LMSE lens is greater than for one or theother 
of the comparison lenses but in a direction which tends 
to make the error curve symmetrical. This can be inter- 
preted as meaning that the LMSE approach has redis- 
tributed the error so as to reduce the peak value. This 
comparison is not a real test of the effectiveness of the 
variational approach since the LMSE lens has two pa- 
rameters, D(y) and 2(y), while the comparison lenses 
have only one each and are not good scanning lenses 
in any reasonable sense. However, these calculations 
indicate that the extremum obtained from the Euler 
equations is certainly not a maximum. 

For the three-parameter lens, where x(y), D(y), and 
n(y) are simultaneously optimized, the obvious com- 
parison lens is the Ruze triple-correction-point design 
since the feed path is circular in both. To this end the 
values of «(y), D(y), and m(y) as given by (8) and by the 


‘appropriate equations in Ruze’s! paper were calculated 


for a total scan angle of 26,=60°. The results are shown 
in Table I and reveal a marked similarity. As might be 


TABLE I 


COMPARISON OF x(y), D(y), AND n(y) FOR LEAST MEAN-SQUARE 
Error (LMSE) Anp TRIPLE-CorREcTION-Pornt (TCP) LENSES 


D—D» nD—noDo 
—x/f os te ee 
Sy A ft be 

LMSE TCR LENSE TCP | LMSE TCP 

0 0 0 0 0 0 0 
00.171 | 00.013 | 00.014 | 00.028 | 00.028 | 00.013 | 00.014 
00.258 | 00.031 | 00.031 | 00.063 | 00.064 | 00.029 | 00.030 
00.329 | 00.052 | 00.052 | 00.106 | 00.104 | 00.050 | 00.049 
00.477 | 00.112 | 00.114 | 00.216 | 00.213 | 00.096 | 00.093 

00.877 | 00.479 = 00.739 = aa — 


expected, the calculated values of 6/f as shown in Fig. 4 
for the two lenses also are nearly identical. It is thus 
apparent that the variational approach leads to essen- 
tially the same lens design as Ruze’s power series meth- 
od. However, the result strongly supports the belief 
that a true minimum mean-square-phase-error design 
was obtained. Moreover, the variational method ap- 
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pears to have no limitations on the focal length-to- 
aperture ratio. 

The “two-dimensional” design equations considered 
thus far are of doubtful practical value except insofar 
as they serve to place Ruze’s results on a firmer theoret- 
ical foundation. The following extension to the three- 
dimensional case is thought to be somewhat more im- 
portant. 


EXTENSION TO THREE-DIMENSIONAL CASE 


The method of analysis can easily be extended to 
the three-dimensional geometry shown in Fig. 5. As 
before, the feed path is a circle in the xy plane with its 
center at the origin. The equations to be determined are 
x(y, 2) for the illuminated surface, D(y, z) for the 
thickness, and n(y, z) for the index of refraction. The 
path-length difference between an arbitrary ray and 
the reference ray through the origin is given by 


6=f+tmD+s—f’ — nD. (10) 
From the figure, 
s = («+ D-— Do) cos6+ y sin @ (11) 
and 
f’' =[(fcoso + x)? + (fsind + y)?+ 2]! ~— (12) 


so that 
b=f—(P+27+ y? + 2 + fe cos 6 + 2fy sin 6)1/? 
+ mDyo — nD + (x + D— Do) cos6+ y sin @ 
= 4[x(y, 2), n(y, 2), D(y, 2), 9, 2, 4]. 
The integral to be minimized is 


Z1 yl 64 
A= fi ih a 6°d0dydz 
=O, 


~— 21. = Y4 


(13) 


= {fo m0. Zz), Diy,2), ny, 2), ¥, 2, 6,|dydz (14) 


ey ra 


in which 
(15) 


The three unknown functions x(¥y, z), n(y, z), and D(y, z) 
can be found by simultaneous solution of 


oH 
aa, (16a) 
dx(y, 2) 

oH 
= = 0, (16b) 
an(y, 2) 

and 

OH 
ah (16c) 
dD(y, 2) 


Carrying out the indicated operations and rearranging 
terms lead to a set of equations which may be made 
identical to (8) by adding terms in 2? to the arguments 
and coefficients of the elliptic integrals as shown in 
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aty.e) 
INDEX OF REFRACTION 


uly,2) 
ILLUMINATED FACE 


Fig. 5—Three-dimensional lens. 


Appendix I. Thus, the algebraic form of the original 
solution is retained. However, graphical solution of (8c) 
for x(y, 2) is considerably more laborious than for the 
two-dimensional case. It should be noted that in the 
three-dimensional case none of the equations have 
rotational symmetry about the x axis since the feed 
path is a curve, not a surface. 

Numerical calculations for the three-dimensional 
case have not been carried out. However, in view of the 
similarity of (8) for the two- and three-dimensional cases. 
and the small phase errors indicated by Fig. 4, it seems 
reasonable to assume that the phase error is, in fact, mini- 
mized by the variational method. 


CONCLUSION 


In general, the analytical method described in the 
foregoing is of greater interest than the specific results. 
In the two-dimensional case Ruze’s far simpler equa- 
tions should be used for design purposes since they 
apparently lead to the same result. In the three-dimen- 
sional case, however, it would probably be impractical 
to determine an optimum design by means of a power- 
series analysis. A double series would be required and 
it is usually quite difficult to evaluate the coefficients 
and to associate them with the physical system being 
analyzed. 

Certain additional extensions, which the authors con- 
sidered but lacked time to carry out, should be noted for 
the benefit of others who might be interested in pursuing 
the method further. 

1) A circular feed path most certainly is not ideal. It 

should be possible and desirable to carry through 
a similar analysis with the feed path as a fourth 
unknown equation to be determined. | 

2) The numerical similarities between the Ruze 

triple-correction-point design and the least-mean- 
square-error design indicate that suitable ap- 


1957 
proximations would reduce (8) to the Ruze equa- 
tions in the two-dimensional case. 

3) In the three-dimensional case it is possible to 
write three independent algebraic equations, each 
of which imposes a requirement for perfect focus at 
one of three distinct feed points. A triple-cor- 
rection-point design can be obtained by simul- 
taneous solution of these equations. It is suspected 
that the resulting equations for x(y, z), 2(y, z), and 
D(y, 2) would be easier to work with than (8) 
and that they would determine a lens which closely 
approximates the optimum design. Numerical 
calculations to verify this hypothesis for a few 
representative cases should prove to be profitable 
and illuminating. 

4) Nonconstraining variable-index-of-refraction ma- 
terials have recently been developed for micro- 
wave lenses. In principle, the path length through 
such a medium can be expressed as a line integral 
involving the unknown functions. The expected 
result of a variational analysis would be a set of 
simultaneous integral or differential equations for 
the determination of the unknown functions in- 
stead of the algebraic equations obtained for con- 
straining lenses. The mathematics would be much 
more difficult to carry out, but the resulting opti- 
mized nonconstraining scanning lens might have 
some useful properties. 

5) The mathematical solution obtained in this paper 
implicitly assumes that all phase errors are equal- 
ly undesirable regardless of their location in the 
radiating aperture. However, since most lenses 
have a tapered energy distribution, it might be 
preferable to modify the phase error by a weighting 
function related to the distribution of energy over 
the aperture. This can be accomplished by multi- 
plying the integrands of (5) and (14) by suitably- 
chosen weighting functions w(y) and w(y, 2) re- 
spectively. This will not change the method of 
determining the lens equations but should alter 
the results considerably. 


APPENDIX | 
EVALUATION OF G= [%,,5°d8 
Squaring (4) gives 

82 = —2(mDo—nD+f)(fP-+x?+y’+ 2fx cos 0+ 2fy sin 6)? 
—2y sin O(f?-+2°+y?+ 2fx cos 6+ 2fy sin 6) 1/2-y? sin? 8 
—2 cos 6(x-+ D— Do) (PP+2?+4°+ 2fx cos 0+ 2fy sin g) 1/2 
+(x-+ D— Dy)? cos? 0+ 2y(x+ D— Do) sin 6 cos 6 
+2y sin 6(noDo—nD-+ 2f) 
+ [2(moDo—nD) («+ D— Do) + 2f(2«x+ D— Dy)| cos 6 
+ (moDo—nD4tfy?t+Ptety’. (17) 


Hence, 
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1 


61 
6°d0= Ay (a+b cos 6+c sin 6)1/2d9 


—4, —o, 


G= 


61 
+As (a+6 cos 6+¢ sin 6)1/? sin 6d0+ A; fe sin? 6d0 
=) 


—6, 1 


01 


1 
+43 f (a+6 cos 6+¢ sin 6)!/2 cosddb-+As f 


cos? 6d@ 
6, —), 
61 41 
+45 f sin 6 cos odd Ax f sin 6d0 
6, —b; 
61 01 
+45 f cos 0d0+ Ag dé (18) 
—O) —O, 
where the upper case letters Ai, As, +--+, represent 
constant coefficients and ; 
a=f?+toet y? (19a) 
b = 2fx (19b) 
c = 2fy. ‘(19c) 


The several integrals, after some manipulation, reduce 
to the following: 


61 


p= (a + 6 cos@ + ¢ sin 6)'/2d6 
—6} 
= 2a + m[E(8, k) — E(—v, k)] (20) 
in which 
m = 2f/x? + ¥? (21a) 
a= taint (21b) 
x 

hp / ls (210) 
t atm : 
ie (214) 

Cee 
ae oe fs (21e) 


and E is the elliptic integral of the second kind. 


61 
(a+ 6 cos 8+ sin @)1/? sin 6d0 
04 


M 


2 cosa [(a + m cos 28)! — (a + m cos 2y)8/] 


ae 1 1 
+4/a+m sina] = = — |e, k) — E(—7, k)] 


1 1 
if E 2 =| [F(8, #) — F(—, ®)] 


vee 2 4/1 a (1 26) 
ES —i—<—- — cos 
+ é sin 26 A 


1 R? 
si es sin 24/1 5 (T= «os m> (22) 
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in which F is the elliptic integral of the first kind. 


41 
(a + bcos 6+ csin @)'/? cos 6d0 


—Oy 


NM = 


2 sin a 
= [(a + m cos 26)?/? 
3m 


— (a + m cos 2y)*!?| 


+ 4\/a + m cos a | 2 —|ee, k) — E(—7,8)] 


3k 


1 1 
+[ 5-5, |1F@» - 8-7) 


~ sin 2384/1 el 28) 
cheat ce ra cos 28 


+ ae 24/1 - (1 Z )} 
Se a —aCOSEAY))i 7a 
6 4 Zz 


The remaining integrals are simple. 
The final expression for G is 


G = — 2L(mDo — nD +f) — 2yM — 2N(x + D— Dp) 
+ y?(@; — sin 61 cos 6) 
+ («© + D — Do)?(61 + sin 61 cos 61) 
+ 4 sin 6,[(moDo — nD)(x + D — Dp) 
Eitan D2.) 
+ 26:[(mDo -nD+fPr+ f+ 2 + 9’. 


The expression for H= [%,,67d0 in the three-dimen- 
sional case is identical to (24) if (19a) is changed to 


(23) 


(24) 


a=Ptart yt a, (25) 
in which case (21c) becomes 
= 
n= 4/ 
a+m 
Afx/x2 v2 1/2 
-( 1 an =< shal (26) 
foe ee te ae DEN Pa? 


APPENDIX II 
EVALUATION OF L’, M’, AND NV’ 


Using the expressions for L, M, and N from Appendix 
I, one obtains 
Go, EE We 


E(y, k) — E(—B, k 
= aa GD 28, 8] 


CSG aren 
ie a a 
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2 2 
es Pe Mae, IR — 5]: 
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f coinage 

3(x* + 7) 
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[R-S] 


[((R + S)x sin 4; 


(29) 
In these equations, 


F=the elliptic integral of the first kind, 
-E=the elliptic integral of the second kind, 
6,=half the total scan angle, 

f=the focal length of the lens, 

a, k, 8, and y are as given in Appendix I, 


R= V1 — k* sin? B, 


and 
S = V1 — R sin? 7. 


APPENDIX III 


EQUATION. OF THE ILLUMINATED SURFACE IN TERMS 
OF @ AND kK 


The transformation of (8c) to the coordinates a and x 
is effected by using the relations 
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x K 
a = cos a tan? oa (30a) 
eee = sin a tan? va. (30b) 
$i 2 
k* = sin? x. (30c) 


After substitution and reduction of the resulting ex- 
pression, one obtains 


em a 


2 sin? 6, 


LL 
2/04 


6, + sin 4; cos #; — 
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K K 
— 2sin@,+ M’ HUGE eg — 26; Sukie = 0, (31) 


where 
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a: 3 cos? ark ee tant K/2 
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Me | 1 - —— ie) — 2-8, 0) 
3 cos? x/2 tan‘ «/2 
sin @ COS a COSK 2 
eel tan* | 
‘[Fly, *) — FCB, 0] 
(inated [1+ costa = |v -m 
3 cos? x/2 tan‘ x/2 
sin @ sin 6; Iv +i COS @ COS ay so Ww): (36) 
3 sin? «/2 3 sin? «/2 
in which 


V = V1 — sin? « sin? B 
W = V1 — sin’? « sin? y. 


The curves x=constant are a family of circles with cen- 
ters at the origin. For the normalized polar-coordinate 
radius defined by 


OGY» 
it can be shown that 
= = tan? 7 (38) 


The coordinate @ is the angle of p/f measured from the 
positive x axis in a counterclockwise sense. These rela- 
tions are shown in Fig. 6. To compute values of x vs y 
take successive integral values of « and find correspond- 
ing values of a which satisfy (31). Each pair of such val- 
ues gives a point on the contour x(y). 


Fig. 6—Illuminated face of scanning lens in terms of coordinates 
cand a. 
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Second-Order Beams of Two-Dimensional Slot Arrays* 
L. A. KURTZ} anv J. S. YEEt 


Summary—Second-order beams in the space patterns of two- 
dimensional slot array antennas may far exceed the sidelobe level 
and therefore increase the susceptibility of a search radar to jamming 
or to spurious target echoes. These beams arise from the nonuniform 
orientation of successive slots;i.e., alternating inclination or displace- 
ment with respect to the waveguide axis. In this paper, expressions 
are derived for the positions of the second-order beams and for their 
amplitudes relative to the main beam. Experimental verification of 
these expressions is presented for two-dimensional arrays of edge, 
series, and shunt slots. Methods for suppressing the second-order 
beams are discussed. 


INTRODUCTION 


HE RADIATION PATTERNS of two-dimen- 
altstea slot arrays frequently contain second-order 
beams whose magnitudes far exceed the low values 
demanded of the sidelobes. Lying outside the principal 
coordinate planes of the main beam, these beams may 
be easily overlooked during radiation pattern measure- 
ments. Because of the present-day emphasis on achiev- 
ing low sidelobes to reduce the susceptibility of a radar to 
jamming, it is imperative that these second-order beams 
be suppressed. 
The source of the second-order beams lies in the 
-small differences in the orientation of successive slots. 
For instance, edge slots and series slots are inclined to 
the waveguide axis in order to couple to the fields with- 
in; however, successive slots have opposite inclinations 
in order to preserve the proper phase relations. Similar- 
ly, shunt slots, which are parallel to the waveguide axis, 
are displaced on alternate sides of the center line in order 
to couple with the proper phase relations to the fields 


within. As a result of this nonuniform orientation of 


successive slots, there are usually certain directions in 
which their radiated fields fail to cancel as expected so 
that second-order beams appear. 

The term “second-order beams” was used by Gruen- 
berg! in a study of linear arrays of shunt slots. He con- 
sidered the slots on one side of the centerline as forming 
one array and the slots on the opposite side as forming a 
secondary array. In the present analysis, only one array 
is assumed, but groups of slots instead of single slots are 
considered as elements of the array. Of course, either 
method should lead to the same results. 

The following study was made to determine the posi- 
tion and magnitudes of the second-order beams as func- 


* Original manuscript received by the PGAP, October 31, 1956; 
revised manuscript received, April 19, 1957. Presented at the Spring 
Meeting of URSI, Natl. Bur. of Standards, Washington, D. C., May, 
1956. This work was performed while the authors were members of 
the Technical Staff, Microwave Lab., Hughes Aircraft Co., Culver 
City, Calif. 

+ Rantec Corp., Calabasas, Calif. 

t Boling Aircraft Co., Seattle, Washington. 

1H. Gruenberg, “Second-order beams of slotted waveguide 
arrays,” Can. J. Phys. vol. 31, pp. 55-69; January, 1953. 


tions of slot geometry and relative phase. Three types 
of slots and two types of arrays are discussed: edge, 
series, and shunt slots with adjacent linear arrays that 
are either identical or mirror images. Experimental 
measurements on arrays of each type of slot are pre- 
sented along with means for suppressing the second- 
order beams. 


ANALYSIS 


The radiation pattern of a two-dimensional array is 
generally described by the product of an array factor 
and an element factor. These factors express the sum of 
the contributions from each of the elements to the field 
in a given direction in space. The element factor de- 
scribes the radiation from a single element, and the 
array factor accounts for the relative amplitudes, 
phases, and positions of the elements in the array. Sepa- 
rating the element factor from the array factor is per- 
missible only if all elements of the array are identical. 

Often, for the purpose of radiation-pattern calcula- 
tions, the slots in an array are assumed to be identical. 
This assumption cannot be made, however, if second- 
order beams are to be taken into account. Instead, it 
becomes necessary to group the slots into the smallest 
number so that successive groups appear identical. The 


radiation pattern of the group then becomes the ele- 


ment factor for the array. The slots within an element 
group, however, are assumed to have equal amounts of 
inclination or displacement and equal amplitudes of 
excitation. These assumptions may be made safely for 
large arrays of many elements where the change in aper- 
ture illumination between successive slots is small. For 
small arrays the assumptions are still adequate if the 
aperture distributions are symmetrical and there are 
even numbers of rows and columns. To allow for the 
variations due to tapered aperture distributions, an 
averaged value of inclination or displacement should be 
used as discussed in the experimental portion of this 
paper. 

A group of four slots has been chosen to form the ele- 
ments in the present analysis. This choice was made be- 
cause four slots are required in the element group when 
adjacent linear arrays appear as mirror images of each 
other (left-hand group at the bottom of Table I). Only 
two slots per group are required when adjacent linear 
arrays are identical; however, this case is included in the 
more general analysis of the four-slot group. 

With four slots combined into one element, the ele- 
ment spacing becomes twice the slot spacing and the 
phase shift between elements becomes twice the phase 
shift between slots. The radiation pattern of a plane 
two-dimensional array may then be written as 
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TABLE I* 
SUMMARY OF SECOND-ORDER BEAM EXPRESSIONS 


Inclined Slots 
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F,(0, 6) =element factor for group of four slots, . 
An, Am=relative amplitudes of the groups. In practical arrays, recurrence of the main beams is 
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avoided by restricting? the values of d, and dz to lie be- 
tween \/2 and X, or by limiting the values of 7: or y2 to 
less than their full range of —7 to +7. Observing these 
restrictions, solutions to (2) can be found for at most 
only three values of p and g, namely, 0 and +1. Thus, 
there are nine possible directions of principal maxima 
of which only one, p=q=0, corresponds to the main 
beam. The other values for p and g correspond to sec- 
ond-order beams. 

Since the array factor portion of (1) has principal 
maxima in the directions of the second-order beams as 
well as the direction of the main beam, the magnitudes 
of the second-order beams are determined almost en- 
tirely by the group element factor. The method followed 
in calculating these magnitudes is outlined in the Ap- 
pendix, while the results of the calculations are pre- 
sented in Table I. Two expressions giving the relative 
magnitudes of the second-order beams have been de- 
rived, one for inclined slots and one for displaced slots. 
The values of » and q corresponding to the types of 
slots are shown at the bottom of Table I. A general pic- 
ture of the positions of the second-order beams can be 
obtained if the slots are assumed to be in phase and 
spaced 0.7. Then for the mirror image array, there are 
four second-order beams approximately centered in the 
four sectors of space in front of the array. For identical 
arrays there are only two second-order beams and they 
‘lie in a plane perpendicular to the plane of the array 
and parallel to the linear arrays. They lie at 45° relative 
to the broadside direction. 

The expressions for Sin Table I contain only the para- 
meters of the array and the values of p and g. The para- 
meters are the slots spacing, the relative phase shift, 
and the angle of inclination or the displacement. Once 
these are specified, the magnitudes of the second-order 
beams can be calculated directly. If the positions of the 
second-order beams are required they may be obtained 
from (2), again in terms of the parameters of the array 
and the values of p apd q. 


EXPERIMENTAL VERIFICATION 


Measurements of second-order beams have been per- 
formed using three two-dimensional arrays consisting 
of identical linear arrays of edge, series, and shunt 
slots respectively. The slot spacings were one-half guide 
wavelength so that the relative phase along the linear 
arrays was zero. The spacings between linear arrays 
were also one-half guide wavelength, and the relative 
phase between linear arrays was varied by changing 
the path length between the coupling slots in the feed 
array. Fig. 1 is a photograph of the edge-slot two di- 
mensional array and Fig. 2 is a radiation pattern, y;= 
72=0, illustrating second-order beams. 

The second-order beams of the edge-slot array were 
measured for several values of phase shift y2. The meas- 
ured values together with values calculated using the 


2S. Silver, “Microwave Antenna Theory and Design,” McGraw- 
Hill Book Co., Inc., New York, N. Y., ch. 9, sec. 9: 19, p. 318; 1949. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


Fig. 1—Sixteen-by-sixteen array of edge slots. The path length be- 
tween feed slots can be varied by means of waveguide inserts to 
obtain a variable phase shift, ye. 
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Fig. 2—E-plane radiation pattern of 16 X16 array of edge slots show- 
ing cross-polarized second-order beams. There is a gap of 0.380 
inch between waveguides. 


expressions from Table I are shown in Fig. 3. Since the 
aperture distribution of the array is tapered, an average 
value of 13.6° was used for B. In order to choose this 
average value, the conditions of broadside operation 
were considered. For these conditions (y= y2.=0, 
qg=0, p= +1, L1=L), the expression for second order 
beams reduces to 


1 \?2 1 
1 — oe Joos: (sin B) sin’ 8 


S = 10 logy (3) 


For small angles the cosine terms approximately cancel 
and the denominator is nearly unity so that the ratio 
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d =1,259 INCHES 
4) = d= 0.880 INCH 
nro 
Bove=!3-6° 
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INTER-ELEMENT PHASE SHIFT, Yp X 180° 


Fig. 3—Second-order beam magnitude vs phase shift for a 
16X16 array of edge slots. 


Fig. 4—Seven-by-ten array of series slots showing the variable 
path length feed. All slots are inclined 11.1°. 


of second-order beam to main beam is a constant times 
sin B. Thus, by adding the second-order beam and main 
beam fields of all the slot groups an average value was 
obtained for sin B and hence for B. 


>> A, sin Ba 


Fiaaie ipa 2 


Eine Dave = 


All the linear arrays are identical, hence B varies only 
along one coordinate of the array and the averaging 
need be done for only one linear array. 

The second-order beams of the series-slot array in 
Fig. 4, were measured, and their values together with 
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Fig. 5—Second-order beam magnitude vs phase shift for a 
6X10 array of series slots. 
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Fig. 6—H-plane radiation pattern of the 7X10 array of series 
slots with uniform illumination: f=8820 mc. 
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Fig. 7—Second-order beam magnitude ys phase shift for a 
6X10 array of shunt slots. 


values calculated using the expressions of Table I are 
plotted in Fig. 5. A radiation pattern showing the cross- 
polarized, second-order beams is given in Fig. 6. Only 
six of the seven rows of slots were used in the second- 
order beam measurements to avoid having a side lobe 
of the main beam lie in the same direction as the second- 
order beam. 

An array of shunt slots, similar in appearance to the 
series slot array, was used in a third measurement of 
second-order beams. The results of these measurements 
are shown in Fig. 7. The differences between measured 
and calculated values can be explained largely by the 
fact that the second-order beams and the side lobes of 
the main beam have the same polarization. It is note- 
worthy that a two-dimensional array of shunt slots 
with broadside beam does not produce second-order 
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beams provided that the beam widths are sufficiently 
narrow and that the linear arrays are identical. 

The displacements of all the slots in the edge slot ar- 
ray are one-tenth inch. If the displacements were un- 
equal, however, an average value for 6 could be cal- 
culated in a manner similar to the calculation of B 
average. Adding the contributions to the second-order 
beam field and main beam field of each element group 
leads to the relation 


oA a SiD & rin) 

m 1 
SS (5) 
ot 
STA, cos i Om 

m 1 


SUPPRESSION OF SECOND-ORDER BEAMS 


The methods used for suppressing second-order beams 
may be simple or complex, depending on the design and 
scanning requirements of the two-dimensional array. 
Perhaps the most obvious method is to eliminate the 
cause of the second-order beams by aligning the slots so 
they are colinear. Colinear alignment of the slots may 
be obtained by using a corrugated-wall waveguide,! 
probe-fed slots,’ or iris-fed slots.4 For the three types of 
slots considered here, such alignment is not possible. 
It was seen, however, that operating near y1=0 effect- 
ively suppresses second-order beams of identical linear 
~ arrays of shunt slots. 

Another method for suppressing second-order beams 
is to restrict the spacing between slots so that the beams 
are outside visible space. This method, when practical, 
serves only for near-broadside operation of large arrays 
with narrow beamwidths and it requires that adjacent 
linear arrays be mirror images. The limitations on spac- 
ing can be readily determined from (2) which gives the 
directions of the second-order beams. For example, con- 
sider the case of broadside operation, yi1=7y2=0, and 
equal slot spacing, d:=d,=d. Then, 


: ait 
sin 6 sin @ maT 
1 

Sree ae (6) 


and these relations can be satisfied only as long as 
|sin 6| is larger than or equal to | cos @). If d is chosen 
such that (6) cannot be satisfied, there will be no second- 
order beam; therefore, d should be less than 0.707. 
(This may require loading a waveguide to reduce its 
physical size sufficiently.) An easy way to visualize this 
situation is to note that (2) represents two sets of cones 
centered, respectively, on the y and gz axis. The direc- 
tions of the second-order beams correspond with the 
intersections of the cones determined by p= +1, 
g=2+1. With the above restrictions imposed on d, 


3 Tbid., p. 301. 

4R. Tang, “A Slot with Variable Coupling and Its Applications to 
a Linear Array,” Hughes Aircraft Co., Report RLM (M) 56-15; 
July, 1956. 
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Fig. 8—Cones represented by (2) with m1=72=0, di=d,=0.6, © 


p=q=1. Since d, and d are less than 0.707, the cones do not in- 
tersect and there are no second-order beams. 


these cones do not intersect, and there are no second- 


order beams. Such a situation is depicted in Fig. 8. 


A method that is suitable for suppressing second- — ; 


order beams of edge slots is to suppress the z component 
of the surface current generated by the slots. The gap 
between the waveguides supporting the edge slots 
forms an effective choke for this purpose if its depth 
(dimension c in Fig. 9) is suitably chosen. A curve show- 
ing the effectiveness of chokes in the array in Fig. 1 is 
given in Fig. 10. In addition to its critical dependence on 
choke depth, the effectiveness of chokes depends also 
on the phase shift between slots, 72, in the manner shown 
in Fig. 11. The reason for this dependence is that the 
field in the choke, E,, consists of approximately equal 
contributions from the slots on either side of the choke. 
Therefore, as the relative phase of these contributions 
tends toward 180° the value of E, tends toward zero 
with the result that there is very little suppression of the 
second-order beams. The situation is reversed if the 
linear arrays happen to be mirror images because then 
the components of the surface current are 180° out of 


phase when y2=0. The chokes are effective in this case — 


only when 72 tends toward 180°. 

In order to eliminate second-order beams for all 
values of y and d and for all types of slots, it is necessary 
to use baffles between slots, as shown in Fig. 12. The 
baffles form parallel plate horns that force the electric 
field to lie entirely in the y direction and to be inde- 
pendent of the slot orientation. Any amount of sup- 
pression can be obtained by adjusting the width and 
height of the baffles. As an example, baffles one-half 
wavelength high, with a 0.39-wavelength gap between, 
reduced the second-order beams by more than 12 db 
over the entire range of 7. 


CONCLUSION 


The second-order beam magnitudes have been deter- 
mined for two-dimensional arrays with mirror-image 
linear arrays and with identical linear arrays consisting 
of edge slots, series slots, and shunt slots. The second- 
order beams are of seriously high magnitudes for all 
cases except the special one of a broadside array of shunt 
slots in which the linear arrays are identical. Equations 
have been derived for the directions and magnitudes of 
the second-order beams, and they are summarized in 
Table I. The equations have been verified experiment- 
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Fig.49—The choke between these waveguides suppresses second- 
order beams when the depth c is one quarter wavelength. At- 
tenuation of second-order beams is greatest for y2=0 and it is 

_ zero for y2=180°. 


SECOND-ORDER BEAM MAGNITUDE, 08 


Fig. 12—Two-dimensional array of edge slots showing 
baffles that suppress the second-order beams. 


CHOKE DEPTH “C”, INCHES 


Fig. 10—Magnitude of second-order beam vs 
choke depth: \=1.259 inches. 
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This may be done conveniently in terms of the magnetic 
: ve radiation vector L. Using the magnetic radiation vector, 
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Fig. 11—Second-order beam magnitude versus phase shift Gye wertich! 2s 


for a 16X16 edge-slot array with chokes. grike 
} Eko=-j ; (Ly cos $) 
ally by using three two-dimensional arrays. Several Ar 
“means are described that result in reduced levels for the _enite ; : 
second-order beams. So far as it is known, however, Ey = j ae (Z, sin ¢ cos @ — L, sin 8), (7) 
baffles are the only means by which these beams can be Sherk 
satisfactorily suppressed for all scan angles. Baffles may Hi 
increase weight, conplexity, and mutual coupling, but as cos E cos a | 
they have yielded the best results so far achieved. Tak §Eim see 
K sin? O 
APPENDIX 
T 
To illustrate the method used in deriving the expres- cos E cos | 
3 : : i a lb foe 4g Em 2 
sions given in Table I, the derivation will be made for pees op (8) 
the combination of inclined slots shown at the lower K sin? O 
left corner of Table 1. The group element factor is ob- 
: : 4 =S, R d J. R. Whinnery, “Fields and Waves in Modern 
tained by adding together, with proper phase retarda- Rede? fone Wiley and Sons, Co New York, N. Y., sec.12:21 


tions, the fields of the four individual slots in the group. _ p. 535; 1953. 
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The radiated field for the four inclined slots is obtained 

by inserting their combined magnetic radiation vector 
in (7). This combined magnetic radiation vector is 

4gEm Snes Bae 

Lye! (kd sin @ sin $—Y1) ef (kd2 cos @—y2) 


Ly = sin B[—L,; — 


+ Lei hu sin @ sin ¢—y1 + Lge? ha2 cos 6-72) | 


4gEm ; : : 
cos BL; + Lye? (kd1 sin 6 sin $—71) ef (kd2 cos 6—~y2) 


(9) 


LI, and L2 are obtained by substituting —B and aa as 
indicated in Table I into the expression 


Tv 
cos E cos 6, | 
2 


L; ee ee a te? a wee 
sin? @ 
where cos 0;=cos (6 cos B;)+sin 0 sin ¢ sin B,. 

It should be noted that the algebraic sign of sin B has 
been accounted for inside the brackets of (9). Hence sin 
B is always positive. Effectively sin B=|sin B,|. 

The expressions for the magnetic radiation vector can 
be simplified for the purpose of calculating second-order 
beams since only directions given by (2) are of inter- 
est. Eq. (2) can be solved for sin @ sin ¢ and cos @. Let- 
_ ting) be the unit measure of d and7z the unit measure of y 


A 2 


28 Lye a1 sin @ sin ¢—y1) -+ Loe? be cos o—y2) | 


(10) 


sin sin ¢@ = = 0, +1 
: 2d, , 
v2 +4¢ 
cos 6 = , = 0, +1. 11 
2d, q (11) 
With the substitution of (11), (9) and (10) become 
4gEn 
L, = cm Bi —Li[t + tor] + Lo[eirm 4 etar]} 
AgEmn 
ja _ cos B{ Lilt + eto] + Zoleior 4 eiar]} (12) 


sin? BL, + 12)*] 1 = ( 


S = 10 logio 
| sin? BL, ae 
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Vuk Pp 
2d, 


12)') 1 _ ( Qf 
2d, 


and 


Odober 
| 
. 


¢ 
> (13) 


In (12) it is apparent that only for p and g not equal to 
zero do the second order beams appear for this slot con- 
figuration. 

The magnitude of the second-order beam relative to 
mainbeam is found by taking the ratio of the electric 
fields radiated in their respective directions. This ratio, 
expressed in db and denoted by S, is 


|Z? + £2|p=+1,g=41 
| Ee + Eo] p = 9 = 0 


S = 10 log 10 (14) 


The expression | E% +E?,| may be simplified by again 
using (11) and neglecting constants that will eventually 
cancel in the ratio. 
| Ee + E,?| 

~ L,? cos? @ + (Ly sin ¢ cos 6 — L, sin 6)? 

~ L,(1 — sin? ¢ sin? 6) + L.(1 — cos? @) 


— 2L,L, cos @ sin 6 sin @ 
Vist y | (* sil ‘| 
~ LPL ——_ J OF 2 P= 
: | ( 2d 2ds 
+ q 
es 2yts(~ *) € as ‘). 
2d i 2d2 
Finally, then, the magnitudes of the second-order beams | 


for the four inclined slots are obtained substituting (12) 
and (15) into (14): 


)] + cos? B(Li — r,)*| 1 S ( 


2 
) | 4 eos? BU Desh 12)'] 1 2 ( A: 
2d» 


(15) 


> 
Nw 
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S| + 
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sin 2B(Ls¢ — Leh) (MP) (ES YI p= Sree 


2d, 
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sin 2B(L.2 — L,”) (=) (= 


(16) 
eee 


Eqs. 13 and 16 appear in Table I. 
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20-70 MC Monopole Antennas on Ground-Based Vehicles” 
ROBERT E. WEBSTER} 


Summary—Antenna design for ground-based vehicular applica- 
tion is discussed with particular emphasis on omnidirectional, verti- 
cally polarized types. Limitations of vehicle size and shape in con- 
trolling the efficiency and radiation patterns are considered. Pro- 
nounced interferences are shown to arise from the usually nonsym- 
metrical currents to the earth. Modifications of vehicle and antenna 
proved effective in reducing this interference. The radiation char- 
acteristics of vertical monopoles excited against the vehicle are shown 
over the frequency band. 


INTRODUCTION 


HIS paper is concerned with communication an- 
_tennas on ground-based vehicles in the frequency 
range where the vehicle dimensions become ap- 
preciable in terms of a wavelength. The vehicle of pri- 
mary interest was the 3-ton army weapons carrier, the 
body of which is elevated several feet above the earth 
and insulated from it with lossy dielectrics (rubber tires). 
Profiles of this truck with the presently used 20-70 mc 
whip antenna are shown in Fig. 1. This antenna and the 
associated radio equipment are designed for short-range 
ground wave communication among a number of such 
units. Vertical polarization and omnidirectional radia- 
tion patterns were therefore of particular interest in the 
investigation. 

The two factors which render the design of antennas 
for this type of ground-based vehicle somewhat less than 
straightforward are the relatively complex geometry of 
the vehicle, and the presence of an imperfectly conduct- 
ing ground. This first difficulty has received considerable 
attention in investigations of aircraft antennas.! Cur- 
rent distribution measurements have revealed the 
nature of favored current flow paths on such struc- 
tures.2-4 Since the contours of truck vehicles are even 
more complex than those of most aircraft, many dis- 
crete paths with high current densities are to be ex- 
pected. The second complication of a lossy material 
beneath the behicle is an added problem for ground- 
based structures. The effect of lossy ground on antenna 
impedances has been treated for certain symmetrical 


* Manuscript received by the PGAP, January 30, 1956; re- 
vised manuscript received, March 15, 1957. This article is based 
on work done under a contract between the Ohio State Univ. Res. 
Foundation and the Signal Corps. Eng. Labs., Fort Monmouth, N. J. 

+ Antenna Lab., Ohio State Univ., Columbus, Ohio. 

1J, V. N. Granger and J. T. Bolljahn, “Aircraft antennas,” 
Proc. IRE, vol. 43, pp. 533-550; May, 1955. } 

2 J. V. N. Granger and T. Morita, “Radio-frequency current dis- 
tributions on aircraft structures,” Proc. IRE, vol. 39, pp. 932-938; 
August, 1951. ; : i 

$1. Carswell, “Current distribution on wing-cap and tail-cap 
antennas,” IRE TRANS., vol. AP-3, pp. 207-212; October, 1955. 

4A. R. Ellis, “Methods of Improving Tail-Cap Antenna Pat- 
terns,” Tech. Rep. No. 35, Contract AF 19(604)-1296, Stanford 
Res. Inst., Stanford, Calif.; January, 1955. 
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Fig. 1—Profiles of ?-ton army weapons carrier 
with standard whip antenna. 


configurations.°* Obstacles to the determination of 
ground current distribution and resultant losses and 
radiation patterns for the vehicular antennas of interest 
here are the elevation of the bodies above the earth and 
their pronounced asymmetries. 

The success achieved in this type of antenna design 
is then dependent on the extent to which the current dis- 
tribution over the vehicle surface can be controlled in 
the presence of those contours which are considered un- 
alterable, and an imperfectly conducting earth. The 
available design parameters are controls of the driven 
elements (number, size, shape, and location) and prac- 
ticable modifications of the vehicle. 


LIMITATIONS IMPOSED BY THE VEHICLE 


In view of the sharply discontinuous metal contours 
of these vehicles, it would seem desirable to avoid the 
problem of controlling the body surface current distri- 
bution. At the high-frequency end of the 20-70 mc band 


5F. R. Abbott, “Design of optimum buried-conductor rf ground 
system,” Proc. IRE, vol. 40, pp. 846-852; July 1952. 

6 J. R. Wait, “Impedance of a top-loaded antenna of arbitrary 
length over a circular grounded screen,” J. Appl. Phys., vol. 25, pp. 
553-555; May, 1954. 
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a choke system, such as loaded dipoles with sleeves, 
for minimizing the antenna-vehicle coupling is feasible. 
Such systems, of course, are frequency sensitive and 
often are not easily adjusted. At the low frequencies, 
however, purely mechanical considerations limit the 
protrusion of the radiators from the vehicle. Then it be- 
comes necessary to use the vehicle as a counterpoise in 
order to obtain reasonable efficiencies, and no effective 
isolation of vehicle and antenna is possible. The pri- 
mary design problem then becomes that of arranging 
the driven element to produce the desired current dis- 
tribution on the vehicle. 

Regarding impedance limitations, it can be shown 
that no great difficulty exists in this frequency range if 
full use is made of the available volume within the 
vehicle contours (assuming, of course, that no attempt 
is made to obtain a highly directive antenna). The rela- 
tively large electrically open regions enclosed only by 
canvas and wood structure on these vehicles permit this 
freedom in impedance design. With compensation at 
the low-frequency end we can expect a standing-wave 
ratio less than 5 to 1 over the entire band. Such imped- 
ances were actually obtained with fan-shaped screen 
antennas placed along the canvas contours.” Impedance 
bandwidths of the presently used whip antennas, which 
consist of thin linear elements protruding to a maxi- 
mum height of approximately seven feet above the 
vehicle, also fall within these limits. 

The rather obvious mechanical limitations of antenna 
elements on ground-based vehicles may be briefly sum- 
marized as follows. 

1) Flexibility and light weight are always desirable. 

2) If the structure is confined within the vehicle con- 
tours, however, it may be reasonably rigid. 

3) Protruding elements, such as whips, must be of a 
smooth, thin shape and highly flexible to prevent 
damage from passing obstacles. 

The tendency of the complex geometry of the vehicle 
surfaces to control the surface current distribution is 
the limitation of primary importance in antenna design 
for this class of vehicles. The results of simple modifica- 
tions of the vehicle contours are shown in the following 
sections, in addition to improvements in radiation char- 
acteristics obtained by major contour changes. 


EFFECT OF THE EARTH 


The close proximity to the earth of current-carrying 
surfaces of ground-based vehicles suggests an influence 
on the efficiency and current distribution of vehicular 
antennas due to the earth’s behavior as an imperfect 
conductor. Measurements of earth constants by Mc- 
Petrie’ yielded the values o = 0.10 mho/meter and €/€) =7 


7D. R. McCoy and R. E. Webster, “Model Techniques in 20-70 
MC Vehicular Antenna Research,” Rep. 522-7, Ohio State Univ. 
Res. Foundation, Antenna Lab., prepared under Contract DA 
36-039 sc42548, U. S. Army Signal Corps Eng. Lab., Fort Mon- 
mouth, N. J.; April 30, 1954. 

8 J. S. McPetrie, “A determination of the electrical constants of 
the earth’s surface at wavelengths of 1.5 and 0.46 meters,” Proc. 
Phys. Soc., vol. 46, pp. 637-648; September, 1934. 
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Fig. 2—Horizontal plane relative electric field intensity patterns for 
standard whip antenna on 3-ton weapons carrier vehicle (see 
Fig. 1). Full-scale and 1/15-scale model patterns. Vertical polari- 
zation. 28 mc full-scale frequency. 


to 16 for the conductivity and relative dielectric con- 
stant, respectively. These values are known to vary 
with soil types and moisture content. Choosing ¢=0.10 
and €/é) = 15 as typical values and computing the earth’s 
intrinsic impedance, we obtain 71:=39.5 |40° as com- 
pared with the value for free space above the earth 
nu=377 | 0°. This impedance discontinuity results in a 


reflection coefficient p =0.85 | 172° for a plane wave in- 


cident normal to the earth from the space above. The 
ratio /we over the 20-70 mc range is seen to be approxi- 
mately 9 to 2.6. These measurements depict the earth 
as a lossy conductor in this frequency band. Thus, its 
presence is commonly expected to produce pronounced 
changes in the free-space radiation properties of vehicu- 
lar antennas. 

The patterns of Fig. 2 indicate this effect to some ex- 
tent although at that frequency the modification was 
not critical. Striking effects are seen in Figs. 3 and 4, 
which show radiation patterns and current density dis- 
tribution measurements obtained for the previously de- 
scribed whip antenna and for a symmetrically located, 
inverted-L monopole on the 1/15-scale model vehicle 
at frequencies where the currents to ground were criti- 
cal. Notice that the vehicle was modified by removal of 
the top cab structure and closure of the gap behind the 
cab so that the vehicle would approximate more closely 
the idealization of an electrically smooth body. The 
notation of the current distribution plots is as follows. 

1) The dot denotes the point of measurement. 

2) The arrow stem shows the axis of current flow. 

3) The arrow head indicates the assumed positive 
sense of current. 

4) The numbers give the amplitude and phase of the 


voltage induced in a small loop held perpendicular to the 
surface. 
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Fig. 3—Horizontal plane radiation pattern and current density dis- 
tribution for whip antenna on modified 1/15-scale model of 2-ton 
weapons carrier vehicle over ground plane. 44 mc equivalent full- 
scale frequency. 


Fig. 4—Horizontal plane radiation pattern and current density dis- 
tribution for inverted-Z antenna on modified 1/15-scale model 
of 2-ton weapons carrier vehicle over ground plane. 50 mc equiva- 
lent full-scale frequency. 


Because of the finite loop size and the very complex 
surface geometries, the loop voltage vs current density 
correspondence is not perfect. Nevertheless, these meas- 
urements have been invaluable in establishing the 
nature of current distributions on these vehicles. 

The deep null off the left front (Fig. 3) is character- 
istic for the whip antenna. Its depth and azimuthal loca- 
tion change with frequency and ground conditions, as 
will be seen in the next section, but its presence and pre- 
sumably its cause is invariable. It is seen from the cur- 
rent distribution plots that the only currents capable of 
producing this destructive interference with the driven 
element currents are those flowing to the ground 
through the wheels and adjacent regions. The largest 
current densities existing in the region of the wheels 
were those shown on the axles and springs. The very 
complex geometry around the wheels and the finite size 
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of the probing loop prevented a more detailed distribu- 
tion plot in these regions. It is clear, however, that there 
must exist vertical components of currents on the vehicle 
surfaces leading to the wheels, springs, and axles. Dis- 
placement currents between the vehicle extremities and 
the ground plane, and currents along the ground plane 
beneath these extremities, must also exist. Taking into 
account the phases, amplitudes, and spacings of the 
axle currents and the driven element currents, one can 
obtain a qualitative horizontal plane pattern of the 
vertically polarized radiation which has the character- 
istics of the measured pattern. The vehicle body thus 
seems to behave somewhat as a linear element loaded 
at its lower extremities. It is apparent, then, that di- 
rectional effects will occur when the spacing between 
the currents in the driven element and the currents to 
ground through these extremities is sufficiently large to 
produce significant phase differences. These model 
measurements were made to simulate 44 mc (the lowest 
frequency of severe destructive interference). The full- 
scale patterns of the next section indicate that over 
actual earth the interference is not as severe at this 
frequency as the model measurements show. The im- 
portant point demonstrated by the scale model meas- 
urements here, however, is that in this frequency range 
there do exist current distributions on the vehicle similar 
to those on linear elements, and more or less discrete 
paths of current flow to the earth. These currents to 
ground are important in determining the radiation pat- 
terns, and their effect on efficiency may be anticipated. 

Similarly the measured pattern of the inverted-L 
monopole may be surmised from the current densities 
measured on the driven element and the region of the 
wheels (see Fig. 4). Consider the geometry to be ap- 
proximated by three point sources in free space as in 
Fig. 5. The end currents represent the radiation from 
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Fig. 5—Point source approximation. 
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front and rear extremities of the vehicle, while the 
center point source accounts for the driven element cur- 
rent. A figure-of-eight type pattern with filled-in nulls 
clearly results from this idealization. This qualitative 
synthesis is not proposed as a substitute for direct pat- 
tern determination. Our choice of the point-source am- 
plitudes, in fact, was prejudiced by the knowledge of deep 
minima in the radiation pattern. These distributions are 
of greatest value, however, in revealing the specific sur- 
faces associated with the interference currents. The 
lowest frequency of major interference is higher than 
for the whip location, as would be expected with smaller 
spacings of the driven element and wheel regions. This 
relatively symmetrical geometry shows also the linear- 
element behavior of the vehicle. At this frequency the 
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body is near a wavelength long, and with the excitation 
employed the phases of the body currents are seen to 
correspond approximately to those of a self-resonant 
linear element. 


SINGLE-FEED VERTICAL MONOPOLES 


Consider now the characteristics of a single monopole 
excited against the 3-ton weapons carrier body, where 
vertically polarized omnidirectional radiation is desired. 
The standard whip antenna is a particular case of this 
single-feedpoint geometry. Radiation patterns for this 
antenna over the 20-70 mc range are shown in Fig. 6. 
The vehicle was located on central Ohio clay earth. 
Note the characteristic minimum discussed in the previ- 
ous section. Its position and depth are seen to change 
with frequency. Minima of less importance are also pres- 
ent, and their effects become more pronounced at the 
higher frequencies. Differences in the grounds and im- 
perfect scaling of the model vehicle tires are offered as 
reasonable explanations of the discrepancies between 
model and full scale results. The observed dependence 
of minima depths on ground conditions is in agreement 
with this explanation. 

To improve the pattern circularity it is necessary to 
reduce the currents to ground in the wheel regions or to 
change their relative phase and/or amplitude with re- 
spect to the driven element currents so that the de- 
structive interference is diminished. A change in length 
of the driven monopole or relocation of it are two obvi- 
ously convenient controls. Fig. 7 shows the results of 
lengthening the whip antenna with a view toward chang- 
ing the relative phases and amplitudes of the interfering 
currents. An improvement is obtained in pattern cir- 
cularity at the previously discussed critical model fre- 
quency. Model pattern measurement at adjacent fre- 
quencies and over the entire frequency range showed 
that this improvement was actual and not merely a shift 
in the critical frequency. As mentioned before, the fre- 
quencies at which destructive interference is critical 
do not correspond exactly between the full-scale and 
model measurements. Effects similar to these produced 
by changing the antenna length on the model, however, 
should occur with modification of the lengths on the full- 
scale vehicle at corresponding critical frequencies. 

The alternate choice of relocating the monopole was 
investigated by using an inverted-L radiator sym- 
metrically located with the feed point near the center 
of the vehicle. This configuration, which fits within the 
vehicle contours, is shown in Fig. 4. As noted in the pre- 
vious section, the lowest frequency where destructive 
interference occurred was raised to 50 mc by this change. 
Noticeable deviation from acceptable circularity through- 
out the entire 20-70 mc range was restricted to the re- 
gion around this critical frequency. The 1/15-scale 
model patterns are shown in Fig. 8 (opposite). Attempts 
to improve the relationship between the driven element 
and lower extremity currents by increasing the length 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


Fig. 6—Horizontal plane radiation patterns for whip antenna (see 
Fig. 1) on weapons carrier vehicle located over central Ohio clay 
earth. Vertical polarization. 


ast 


Fig. 7—Horizontal plane radiation patterns for whip antennas on 
1/15-scale model weapons carrier over ground plane. 44 mc equi- 
valent full-scale frequency (the frequency of critical interference). 


of the monopole resulted in the patterns of Fig. 9. This 
type of monopole thus shows promise of operating as 
an effective omnidirectional antenna over a broad fre- 
quency band. The improvement was again found to be 
actual and not merely a small shift in the critical fre- 
quency. For the longer monopoles (J=136 inches) 
other frequencies of critical interference were noted 
above 66 mc. Two or three different lengths would thus 
be required to produce omnidirectional patterns over the 
frequency band. 

Notice that these results were obtained on a vehicle 
with considerably modified contours, the most notable 
of which were removal of the top cab structure and 
closure of the gap between the cab and rear bed. An 
additional modification of the lower surfaces of the 
vehicle showed appreciable effects—all cylindrical ele- 
ments and gaps (axles, drive shaft, springs, etc.) were 
bonded to the body resulting in improvements of the 
patterns of Figs. 3 and 8. The minima were generally 


October 


60 me 65 me TO me 


4 Fig. 8—Horizontal plane radiation patterns for inverted-Z mono- 
pole (Fig. 4) on 1/15-scale model weapons carrier vehicle over 
ground plane. Equivalent full-scale frequencies. 
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_ Fig. 9—Horizontal plane radiation patterns for symmetrically located 
vertical monopole on 1/15-scale model weapons carrier vehicle 
over ground plane. 49 mc equivalent full-scale frequency (the 
frequency of critical interference). 


filled in and shifted to higher frequencies, but omni- 
directional patterns were obtained again only over a 
portion of the frequency band when the monopoles were 
restricted to one length. The vehicle in this stripped- 
- down condition corresponded more closely to the ideali- 
zation of an electrically smooth body, although high 
current densities were still noted along the outer edges 
of the body surfaces and on the lower extremities. In 
the final design of such a monopole or other antenna on 
this type of vehicle all of these modifications may not 
be necessary. The changes mentioned, however, were 
found to have critical effects on radiation patterns of 
various radiators on the vehicle. Extensive bonding of 
all joints, gaps, and cylindrical elements separated from 
the body has proved generally desirable. Behavior of 
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the steering column and adjacent body surface as a 
resonant circuit producing nulls in the pattern in one 
instance is an example of the local resonances which 
must be suppressed in this type of vehicle. 


MULTIPLE FEEDS 


The radiation characteristics shown for vertical 
monopoles excited at a single feed point indicate that 
this configuration is at a basic disadvantage. Numerous 
interference possibilities exist among such a discrete 
driven element and the various current paths to ground 
along the bottom of the vehicle. It should be possible to 
resolve this interference problem by modifying the dis- 
creteness of either or both of these current regions. 

This solution has not yet been applied to the ground 
current paths. Mechanical problems are immediately 
presented in this situation. Extensive bonding of the 
bottom surfaces is a practicable first step toward redis- 
tribution of the ground currents, but it proved to be in- 
sufficient with the antennas discussed before. Efficiency 
and interference considerations, however, make redis- 
tribution of ground currents particularly attractive if 
practicable means can be found to change the bottom 
surfaces of these vehicles to more symmetrical struc- 
tures. 

Means of “spreading” the feed system and the result- 
ing effects are more obvious. An example of this ap- 
proach is the structure sketched in Fig. 10. Here all of 


Fig. 10—Sketch of multiply fed “spread” monopole 
on weapons carrier. 


the vehicle contours have been metallized, and the en- 
tire top half of the truck is excited uniformly as an ex- 
tended monopole. The phase difference between the 
distributed feed currents and currents to ground should 
be minimized. Radiation patterns obtained for this 
configuration before and after bonding the lower vehicle 
surfaces are shown in Figs. 11 and 12. The characteristic 
figure-of-eight pattern found at 50 mc for both this 
distributed monopole and the inverted L singly-fed 
monopole is seen to be eliminated here by the lower sur- 
face bonding. These patterns do not yet attain the de- 
sired omnidirectional qualities over the entire frequency 
range, but further modification of the vehicle’s lower 
surfaces may extend the pattern bandwidth. 
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Fig. 11—Horizontal plane radiation patterns for multiply fed mono- 
pole (Fig. 10) on 1/15-scale model weapons carrier vehicle over 
ground plane. Equivalent full-scale frequencies. 
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Fig. 12—Patterns for antenna of Fig. 10 continued. Lower extremi- 
ties of vehicle model extensively bonded to body. 


CONCLUSION 


Some aspects of antenna design for ground-based 
vehicles in the high-frequency range have been discussed 
In this range, where the vehicle dimensions vary from 
an appreciable fraction of a wavelength to several wave- 
lengths, current distributions on the vehicle surfaces 
similar to those of linear elements in free space have 
been shown to exist. These distributions and their effects 
on the antenna characteristics were seen to depend on 
the surface contours of the vehicle, properties of the 
adjacent earth, and design and location of the driven 
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Ober 
elements. The only practicable means of predicting an- 
tenna behavior of these complex geometries is through 
systematic measurements of surface current densities, 
radiation patterns, and impedances for driven elements 
which appear likely to give the desired properties in the 
presence of the vehicle body which is considered as an 
isolated surface with loading corresponding to ground 
effects. 

Omnidirectional horizontal plane patterns were found 
to be obtainable with a series of vertical monopoles sym- 
metrically located on vehicles modified to approximate 
an electrically smooth body. Persistence of the omni- 
directional property throughout a large frequency band 
is dependent on control of interference between vertical 
currents on the monopoles and on the vehicle surface in 
the regions of favored flow paths to ground. Current 
paths to ground were shown to be more or less discrete 
with relatively high densities in the regions of the 
wheels. Interference between these currents and those 
of the driven elements is pronounced in some cases. 
Length and location of the driven elements are conven- 
ient controls over these effects. Two alternate systems 
for reducing this interference problem were suggested: 


1) Isolation of the antenna from the vehicle in the 
higher frequency region by “choke” arrangements 
to impede the antenna current flow onto the vehicle 
surfaces. 

2) Modification of the lower extremities of the ve- 
hicle and use of multiply fed driven elements to 
destroy the discreteness of the interfering currents. 


Patterns obtained with a multiply fed monopole indi- 
cate the improvement possible with the latter system. 

Modification of the vehicle surfaces was shown to be 
necessary in certain cases. Presence of the top cab struc- 
ture and the centrally located gap between the cab and 
rear bed of the #-ton army weapons carrier vehicle was 
particularly detrimental to the omnidirectional charac- 
teristics of a symmetrically located vertical monopole. 
Extensive bonding of the complex lower surfaces 
(springs, axles, drive shaft, etc.) also produced improve- 
ments. All practicable changes in contours to eliminate 
resonant effects and render the vehicle an electrically 
smooth symmetrical body are considered generally de- 
sirable. 

Results of these investigations apply generally to 
other ground-based vehicles. The basic problem with 
any radiator excited against the body is the nonsym- 
metrical nature of the currents to ground. Thus, it is 
unrealistic, perhaps hopeless, to exclude consideration of 
the vehicle’s lower surfaces and currents to ground from 
the design of any hf antenna for ground-based ve- 
hicles. 
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New Method of Antenna Array Synthesis Applied to Generation of 
Double-Step Patterns* 
C. J. SLETTEN}, P. BLACKSMITH, JR.t, anv G. R. FORBES, JR.t 


Summary—A method of synthesizing linear antenna arrays uti- 
lizes the natural phase distribution that exists on transmission lines. 
The far-field pattern function obtained with equispaced elements is 
decomposed in terms of the feeding coefficients realizable along the 
feeding guide. Such functions are sufficiently general to permit syn- 
thesis of any arbitrary pattern with arbitrary phase. This theory is 
applied to an asymmetric pattern. The engineering techniques 
needed to design this pattern are described. Experimental results 
are in good agreement with the theory. 


INTRODUCTION 


EVERAL methods of designing linear antenna 
S arrays for achieving asymmetrical functions such 

as csc? @ patterns are in use.!~* For deriving feeding 
coefficients that will produce a given power plot, a 
method developed by Woodward? is often applied be- 
cause the feeding coefficients can be realized on a feeding 
waveguide. The Woodward technique, however, in- 
volves approximating a continuous aperture distribution 
by means of discrete radiating elements and, since phas- 
ing is approximately obtained by varying the space be- 
tween elements, pattern calculation is difficult. 

The method presented here is based on a distribution 
of discrete radiators whose elements are uniformly 
spaced about ),/4 apart.* For the double-step pattern 
design described in this paper, almost half of the pos- 
sible elements can be omitted, leaving most of the spac- 
ings A,/2 apart. 


METHOD OF PATTERN SYNTHESIS 


We shall restrict our discussion to a treatment of the 
problem in one dimension, the @ plane. Fig. 1 shows 6 
measured from the normal to the array. The radiation 
pattern produced by a linear array of omnidirectional 
elements will always have symmetry that is even with 
respect to the line of the array unless the elements are 
made directive by a ground screen or some other device. 

The required power pattern E(0)-E*(6) is specified 
in the interval —90°<@<90°. Weshall assume that indi- 


* Manuscript received by the PGAP, February 27, 1956. 

+ Air Force Cambridge Research Center, Bedford, Mass. 

- 1S, A. Schelkunoff, “A mathematical theory of linear arrays,” 
Bell Sys. Tech. J., vol. 22, pp. 80-107; January, 1943. 
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Fig. 1—Linear array with ground screen. 


vidual radiating elements produce identical radiation 
patterns and the same polarization at each point in 
space; i.e., E(0) =f(6)£, where £& is a unit vector in the 
direction of polarization, and the complex function f(6) 
is the scalar function to be generated. The complex 
function f(@) is ordinarily |f(@)|e7™, with the phase 
function 7(8) being adjusted for different applications or 
to simplify design. 

The far-field pattern of a linear array of equally 
spaced radiators, where the directivity of each element 
is D(@), is given by . 

n=N 


(0) = D(@) > Ryei entnd(2m/d)¥) 
n=—N 


(1) 


where 


R, is the amplitude of excitation 
gn is the phase of excitation 

d is the wavelength in free space 
d is the element spacing 

y is equal to sin 0. 


Fig. 2 illustrates how f(@) can be resolved into com- 
ponent functions that are simply related to phase on the 
feeding line or guide. On a transmission-line shunt 
loaded with real admittances at nodes or antinodes in 
the standing-wave pattern, the phase increases 90° 
from the feeding end at each successive maximum and 
minimum of the transmission-line fields. When 180° 
phase reversal of the radiator is possible, by crossing the 
leads on a dipole or, as shown here, by placing longitu- 
dinal shunt slots on opposite sides of a rectangular wave- 
guide in the TE: mode, the values of e’# can be 1, —1, 
7, and —j. 

In Fig. 2 the phase reference is the slot in the center 
of the array. The feeding coefficients An, Bn, Cn, and D, 
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—=7 zt ei (amY/21) 
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VAG = > Cy cos ou 


n=0 1 


N 
fal) = — 52 Dy sin ™ 


fu) = ~ 4X An cos (n+5)™ 


fl) = 2 B, sin (1 +5) = 
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Fig. 2—Resolution of radiated pattern in terms of phase on transmission line. 


are assumed to be real, and are adjustable by the cou- 
pling conductance g of the radiator. The phase of the ele- 
ment is obtained from its longitudinal position and its 
location right or left of the center line. Only the first few 
elements (terms) are shown, but it is clear that more can 
be added, in pairs, on the ends of each elementary array. 
The term under each radiator represents the pattern 
contribution by that element. Elements are identified 
and summed to illustrate how they can be combined 
into four simple sine and cosine series. Each of these 
four series is in the form of a Fourier series in y=sin 0 
in the interval —/ tol. It is clear that superimposing the 
four elemental slot arrays that correspond to the four 
series will result in longitudinal slot separations of about 
d,/4 where X, is the guide wavelength. 

Any ideal required pattern f(@) can be decomposed 
into even and odd parts since electromagnetic fields are 
linear and can be superposed on one another. Hence 


f(@) ideal — Lccamewen(a) + freat,oaa(6) | 
=F A Liisag’ even) oe Fitsageaaie) (2) 


feven = 3[f(0) + f(—8)]; foaa = 3[f(6) — f(—8)]. (3) 
We can let 
ae ee 2, C, cos — (4) 
n=N 1 
Jrol@) = full) = By sin (> i =)= (5) 
n=N 1 
fie) = jful) = v2, A, COs (» + =)= (6) 
fuo(@) = ifel) = LD Da sin (7) 


The degree of approximation to f(@), of course, wil 
depend on the number of terms in the series, that is, on 
the number of radiators used. Of the various types of 
approximations that could be used over the interval —/ 
to 1, approximation in the least square sense is the most 
convenient. For this type of approximation the coeffi- 
cients A,, Bn, C,, and D, are merely Fourier coefficients 
for which we have the following explicit formulas: 


1 . 1 
A,.= Sf fo cos (» = > dy, 
1 d 1 
BS —f Fr,o(@) sin (» + = dy, 
MS ae. 
1 l 
Cn = Fil ve 98) COS os dy, 
IL J_j l 
nay 


1 l 
DR = rs Fi,o(8) sin dy. 
—l 


To check orthogonality, we note that 


_ (m+ 3)m) © 
sin ————_—— 
i [eae OER iil ak 
oh 1 (m+ arb = 0, mA Nn. 
Sede? 4 Cha 


It is evident that this decomposition is sufficient to ap- 
proximate any arbitrary pattern f(0) in the least square 
sense. 

We have developed f(@) in a least square approxima- 
tion in Y=sin 6. If we want a least square approxima- 


tion in 6 we must evaluate a system of equations where 
integrals like 
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sin 17 
nN 
C f _, £08) sin E sin | dp 
sin —l 


are typical. These can usually be expressed in Bessel or 


Struve functions when f(6) is simple. 


Only when /=1 or d=)/2 will the expansion interval 
agree with the physical interval —7/2<@<7/2. When 
L>1, d<X/2, part of f(@) will not be visible. When J< 1, 
d>x/2, f(@) will start repeating when @=sin-/, and 
this will generally be the case when slotted waveguide is 


used. 


APPLICATION OF METHOD TO THE DOUBLE- 
STEP FUNCTION 


Our experimental efforts have been limited to gen- 
erating the ideal function shown in Fig. 3. We have used 
both waveguide and two-wire line-feeding systems. For 
most of the work presented here we used 1/16-inch slots 
on 1 inch by 3} inch waveguide, at a wavelength of 3.2 
cm. The double-step function was selected because it 

“minimized the number of \,/4 spaced elements and 


; their resulting mutual couplings. 


The pattern in Fig. 3 can be approximated by a co- 
sine symmetric main lobe and a step function on one 
side of the 0=0 axis. To simplify the design we let 
-fw(@) produce the main lobe and let a combination of 
f.(@) and f.(@) produce the step function. We thereby 
obtain the odd and even patterns in phase quadrature 


where they do not destructively interfere with each 
_ other. 

Asa first step, suppose we make f, produce an odd 
function as shown in Fig. 4. Then 


n=N 1 my 
ri, = D Basin (n+ —)™ a, 
< n=0 2 l 
1 aay 1\ ry 
Pe netoTt fis ima SID (n+ 5)" a 
hale e a2 Pe Dy! 
1 eri ( ~)= 1 
= — sin (n+ —)— dy = ——— : (8) 
aeisthis2 2) 1 (n + dv 
n=0 Bo = 2/r = 0.6366 
nm = 1 By, = 2/3% = 0.2122 
n= 2 By = 2/50 = 0.1273 
i — Bz = 27 1z-= 0.09095: 


Note that when f,,(W) generates a constant function 
equal to ?, the series reduces to one term, Co =F. 

Now, to produce a step function of unity amplitude 
‘in the interval between 0<0<z/2 and approximately 
zero in the interval —71/2<0<0, we can write for the 


first seven radiators: 


bry 3k 
= fot fu = 0.5 + 0.637 sin — + 0.212 sin 


teten 


(9) 


Skary 
+ 0.127 sin 


Sletten, Blacksmith, and Forbes: Antenna Array Synthesis 


371 


| 
| 
f(@) | 6:3 


1 y= sin 6 ! wet 
Fig. 3—Ideal field pattern. 


Fig. 4—Desired characteristics of odd function. 


where k= 2d/\ and d=)g/2. The first slots on each side 
of the phase reference center are spaced ),/4, and all 
other slots are spaced A,/2. 

The peak quadrature pattern can be derived by using 
Tchebycheff coefficients to get low side lobes. We com- 
puted these and for a 32-db side lobe level obtained 


Ao = 3.755; A; = 2.502; and A» = 1.000. 


These values result in an equation of the form 


k 3h Sk 
FGAY =2SeTRS cos =P42502 codes hooness ae 


Since the peak value of the pattern is to be 16 db 
above the step function, the coefficients for f..(¥) have 
to be modified by a multiplicative constant. By setting 
¥=0 in fu(W) and adjusting the sum at 6.3, we get the 
desired 16 db above the 1.0 value of the function fetep (W). 
The resulting equation for f,(W) is then 


kaw 3k 5k 
fowl) = 3.259 cos wy + 2.170 cos ue + 0.868 cos = 
We now calculate the field (£) and power (P) excita- 
tion for the final array (see Table I). With the distribu- 
tion given in Table I, the patterns of the individual 


TABLE I 

FIELD AND POWER EXCITATION OF DOUBLE-STEP ARRAY 

Slot Step Peak Relative 

No. Eraa Ered Total Erna Phea 
1 —0.0636 0.434 0.370 0.137 
2 —0.1061 1.085 0.979 0.860 
3 —0.3183 1.625 1.307 1.708 
4 0.5000 — 0.500 0.250 
5 0.3183 1.625 1.943 Baths 
6 0.1061 1.085 1.191 1.418 
7 0.0636 0.434 0.498 0.248 
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radiating elements combine to produce the total field: 


k 3k 
fiy) = D(®) | 0.500 + 0.637 sin = + 0.212 sin ey 


Skrp 


+ 0.127 sin 


k 3k 
+5 (3.259 cos ee aes 


“)] 0 


where D(@) is the pattern in the @ plane of the individ- 
ually radiating elements, and y=sin 6. For more ef- 
fective coupling in the actual array, the slot spacings 
were reduced so that 


+ 0.868 cos 


i,= 0.9/2 
and 


0.9X¢ 


k = 1.26. 


A theoretical plot of the function is shown in Fig. 5 
with D(@) =1. Fig. 6 identifies the slots and shows their 
positions. 


METHOD OF DETERMINING SLOT COUPLINGS 


Several difficulties hinder realization of the desired 
feeding coefficients. First we must have room to accom- 
- modate slot no. 4 (see Fig. 6) without having to cut 
into slots no. 3 and no. 5. This is easy since slot no. 4 
is coupled low enough to be cut beside either no. 3 or 
no. 5, regardless of the type of load at the terminal of 
the array. Slot no. 5, however, couples considerably 
more power than the other elements. We must be care- 
ful to achieve its relative power of 3.775 because realiz- 
able values of slot conductance g are less than 1.0. The 
gain of the array will be reduced if we dissipate enough 
power in a matched load to reduce all the couplings. 
It is a simple matter to find the values of g, if we assume 
a matched terminal load and no reflected power at any 
slot. This is a good approximation on multislot arrays 
where individual couplings are low but a more general 
approach is necessary when the power to be coupled is 
large and the number of slots to do the job is small. Our 
design formulation must be general enough to allow for 
a change in frequency, which makes the slots reactive 
and effectively changes the slot spacings. 


General Procedure 


The general procedure (see Fig. 7) for finding the 
value of the mth slot admittance (yn) is to determine the 
total admittance (yz,_,) at the position of the preceding 
radiator (aa’) and reflect this admittance to position 
bb’, by using the following expression for normalized 
values: 


Yin=1 7 % tan Bh 
1b 71,1 tanpn 


Vn in = 
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RELATIVE POWER ONE WAY (aD) 


Fig. 5—Theoretical pattern of seven-slot array double- 
step function with D(@) =0.45 \, spacing. 


Slot Slot Slot 

No. g Coupling | Disp Length | Width | Remarks 
(in.) (in.) (in.) 

gi | 0.0457 | 0.1953db | 0.0630 | 0.621 | 0.0625 | load end 

ge | 0.3250 | 1.2222db | 0.1360 | 0.627 | 0.0625 

gs | 0.5170 | 1.8099db | 0.1980 | 0.632 | 0.0625 

gs | 0.0249 | 0.1072db | 0.0480 | 0.620 | 0.0625 

gs | 0.8070 | 2.5696db | 0.2580 | 0.637 | 0.0625 

ge | 0.1850 | 0.7372db | 0.1160 | 0.625 | 0.0625 

gz | 0.0240 | 0.1030db | 0.0475 | 0.620 | 0.0625 | feed end 


Fig. 6—Slot location of seven-slot array double-step function 
with 0.45 \, spacing and matched load, \=3.2 cm. 


8 =27/,, h is distance between shunt elements, and the 
subscript 7m means the admittance looking toward the 
right before y, is added. The power dissipated in y, is 


P, = Eg,. 


The power dissipated to the right of yn is 
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n=l 


by P; = Eg, in: 


7=0 


Eliminating E? we see 


£n in 


8n in are 


This procedure can be started at the load end of the 
array and repeated until all the values g, are found. The 
_ computation of gin at each junction is simplified by the 
use of a Smith chart. 


EXPERIMENTAL RESULTS 


___In our study of the loading problem it was seen that a 
- short at \,/4 from slot no. 1 tends to establish all slots 
- except no. 4 in regions of relatively low guide conduct- 
ance. This gives maximum radiated power with mini- 
- mum slot conductance (small transverse slot displace- 
ment) on the highly coupled slots but throws slot no. 4 
- into a region of maximum slot conductance so that its 
phasing in relation to the rest of the array becomes 
more critical. Several such models tested gave poor pat- 
_ tern results. A five-slot array with resonant spacings 

and matched load absorbing about 50 per cent of the 
input power produced the pattern shown in Fig. 8. 
- Although patterns were best at \=3.4 cm instead of at 
\=3.2 cm, they were fairly good over a wide band of 
- frequencies. 

In order to obtain a beamwidth of approximately 20° 
and deliver a relatively small fraction of the power in the 
matched load, an array was constructed to the following 
specifications. A matched load was used but the spacings 
were 0.225\, and 0.45,. This placed slot no. 5 in a region 
of lower slot conductance where it was possible to get 
high power out of it. Fig. 9 is a photograph of the array, 
and Fig. 6 is a diagram of the slot locations. Fig. 10 
shows the patterns obtained experimentally at A= 3:2 
cm. The agreement with the calculated pattern in Fig. 5 


is good. 
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Fig. 9—Seven-slot array for double-step function. 


RELATIVE POWER ONE WAY (db) 


Ord ANGLE 36° 72 108" 


Fig. 10—H plane of seven-slot array with matched 
load on the end, 0.45 \, spacing. 


DISCUSSION 


An important consequence of this method of array 
synthesis is that arbitrary patterns can be achieved by 
varying the amplitude of radiations at each element 
while maintaining fixed phase relationships between 
radiators. 

The synthesis of complex antenna patterns by a 
method that utilizes the properties of transmission lines 
admittedly presents experimental difficulties. But these 
same difficulties exist in any method of pattern synthesis 
used to determine power and phase distribution on an 
array. Our method requires uniform spacings of ap- 
proximately \,/4. The bandwidth will probably be less 
than with arrays having \,/2 spacing although experi- 
mental results indicate that 10 per cent bandwidths can 
be achieved. 

Mutual couplings between elements create the big- 
gest problem in realizing feeding coefficients. It is known 
that the curve for intercouplings between parallel elec- 
tric dipoles monotonically decreases as the spacing in- 
creases.’ Conditions at \,/4 are not much worse than 
those at A,/2, and it seems plausible that suitable corner 
reflectors behind the radiating dipoles would diminish 
exterior couplings. 

We believe that a joint experimental and theoretical 
attack on the combined problem of both the exterior 
radiation couplings and the interior transmission line 
couplings will produce solutions that will enable us to 
maintain precise control of the radiation pattern.® 


5 R. W. P. King, H. R. Mimno, and A. H. Wing, “Transmission 
Lines, Antennas, and Waveguides,” McGraw-Hill Book Co., Inc., 
New York, N. Y., p. 125; 1945. 

6B. J. Starkey and E. Fitch, “Mutual impedance and self-im- 
pedance of coupled parallel aerials,” J. TEE, pt. 111, vol. 97, pp. 
129-137; May, 1950. 
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Convergent Representations for the Radiation Fields from Slots 
in Large Circular Cylinders* 
L. L. BAILIN} anp R. J. SPELLMIRE}{ 


Summary—Expressions are derived which will approximate the 
radiation fields from slots in large circular cylinders. This derivation 
is accomplished by applying a Watson transformation to the slowly 
convergent harmonic series representation and then evaluating the 
resulting integral to third order with the techniques of Van der Pol 
and Bremmer. Highly convergent representations are obtained which 
are valid in different azimuthal sectors around cylinders where ka is 
in the range of 10 to 1000 or more. 


I. INTRODUCTION 


T WAS pointed out in a recent report by Sensiper! 
if and by others? that cylindrical wave representa- 
tions in terms of harmonic components had very 
limited applicability in the case of large cylinders. The 
reason for this limitation upon the applicability of the 
harmonic representation is contained in the fact that 
these representations are slowly convergent series as 
functions of ka where k=27/) and a is the radius of the 
cylinder. 
The basis for this present study of more convergent 
representations has been described by Van der Pol and 
Bremmer?® in their classic studies of terrestrial radio 
waves. Sensiper! applied their techniques to cylindrical 
waves and investigated the radiation fields from both 
axial and circumferential slots in large cylinders (see 
Fig. 1 for geometry). However, the mathematical tech- 
niques described therein may be applied with slight 
modifications to the analogous dipole and scattering 
problems. In the principal plane, 0=7/2, he derived 
the second-order approximations (sometimes referred to 
as the geometrical optics representations) in the so- 
called optically illuminated region where the azimuthal 
angle @6<7/2, and the residue series representation in 
the shadow region where r>¢>7/2. Although these 


* Manuscript received by the PGAP, July 7, 1956; revised manu- 
script received, June 28, 1957. The work described in this report was 
sponsored by the Air Force Cambridge Research Center under Con- 
tract AF19(604)-262, Change C. Much of the material discussed was 
presented at the URSI Symposium on Electromagnetic Wave Theory, 
University of Michigan, Ann Arbor, Mich.; June, 1955. 

+ Dept. of Elec. Eng., Univ. of Southern Calif., Los Angeles, 
Calif.; formerly at Hughes Aircraft Co., Culver City, Calif. 

t Hughes Aircraft Co., Culver City, Calif. 

1S. Sensiper, W. G. Sterns, and T. T. Taylor, “A further study of 
the pattern of single slots on circular conducting cylinders,” presented 
ee URSI-IRE Convention, NBS, Wahsington, D. C.; April, 
1952. 


S. Sensiper, “Cylindrical radio waves,” IRE TrRans., vol. AP-5, 
pp. 56-70; January, 1957. 

2 W. Franz and K. Depperman, “Theories der beugung am zylin- 
der unter beriicksichtigung der kriechwelle,” Ann. Phys., vol. 10, 
pp. 361-373; 1952. 

W. Franz, “On the Green’s functions of the cylinder and the 
sphere,” Z. Natur., vol. 9a, pp. 705-716; September, 1954. 

I. Imai, “The diffraction of electromagnetic waves by a circular 
cylinder,” Z. Phys., vol. 137, pp. 31-48; 1954. 

3B. Van der Pol and H. Bremmer, “The diffraction of electro- 
magnetic waves from an electrical point source round a finitely con- 
ducting sphere,” Phil. Mag., vol. 24, pp. 141-176; July, 1937. Also see 
H. Bremmer, “Terrestrial Radio Waves,” Elsevier Press, Inc., New 
York, N. Y.; 1949. 


Fig. 1—Geometry of the cylinder and slot. 


second-order results would prove to be adequate when 
compared with the harmonic series representation of the 
fields from slots on cylinders where ka = 100, they would 
deviate noticeably for smaller values of ka or for the 
fields off the principal plane which depend on ka sin @. 
Consequently, more accurate third-order results would 
be immediately useful in the calculation of element fac- 
tors for linear slot arraysin a range of ka which is cur- 
rently practical. These results would also be useful for 
a complete description of the three-dimensional field. 

The general method to be used in the derivation of the 
third-order approximations will be to choose the prin- 
cipal plane so as to simplify the geometry to two dimen- 
sions (r, @) (this method was used by Sensiper et al.1). 
This simplification is not, however, an essential feature 
of the procedure. The three-dimensional behavior can 
be obtained by replacing ka with ka sin @ whenever it 
appears in the argument of a Hankel function and multi- 
plying by a factor which is a function of 6 only. Detailed 
descriptions of slot shapes, methods of excitation, and 
all but the principal field component are nonessential 
to this study. For instance, the convergence difficulties 
encountered in the cross-polarized component of the 
radiation field from arbitrarily shaped slots on large 
cylinders‘ are the same as those described in Section III 


4S. Silver and W. K. Saunders, “The radiation from a transverse 


rectangular slot in a circular cylinder,” J. Appl. Phys. 
745-749; August, 1950, P Pe Oo eer 
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for the narrow half-wavelength axial slot. Thus, this 
paper will be concerned primarily with the principal 
components of half-wavelength delta slots which are 
parallel and perpendicular to the axis of a large cylinder. 
By an appropriate superposition of these results, the 
radiation characteristics of an arbitrary slot may be 
found. 


II. CIRCUMFERENTIAL SLOT 


It can be shown‘ that in the principal plane the 
main component of the radiation characteristics from a 
narrow circumferential slot may be derived entirely 
from the electric Hertz vector given by 


Hy) (kr) 


T,=K aera Ss 
rata En” Ra) 


p [—jno]S(n, $0) (1) 
where K=|Eo| /2¢ok? (| Eo| = maximum electric field 
_ strength of a center-fed slot), and the function S(m, ¢o) 
depends on the source characteristics. For example, 
the source function for a half-wavelength (¢) =7/2ka) 
delta slot (J—-0) excited by a sinusoidal field E(¢) = Eo 
cos (kad) is (see Bailin®) 


(ka)? cos (=~) 


(ka)? — n? 


S(n, cD) = 


Eq. (1) leads to the harmonic representation for the 
radiation fields in the principal plane. It may be modi- 
fied to include the behavior of the fields off the principal 
plane by multiplying by (1/sin @) and inserting ka sin @ 
as the argument of the Hankel function in the denomi- 
nator. Such modifications, however, will not affect the 
general procedure which follows. 

When the technique of Van der Pol and Bremmer? is 
used, the summation in (1) is converted to an integral 
by a Watson’ transformation. Eq. (1) becomes 
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U- PLANE 
Ce 


Fig. 2—y plane with contours C, and C. 


v=n for n=0, +1, 42, +3,--- then (2) yields the 
harmonic representation of (1). However, if the integral 
is evaluated at the poles determined from H, (ka) =0, 
then the residue series representation which will be 
discussed in Part B is obtained. The third possibility is 
to evaluate the integral along the path C, without re- 
course to the theory of residues. This method will lead 
to the geometrical optics results. 


A. Third Order Geometrical Optics Representation 
(|| <z/2) 
To obtain the geometrical optics results, the contour 


C, of Fig. 2 was considered and the integrals modified as 
follows: since in the upper half-plane 


= — 27 © exp [jv(2m + 1)z], 


sin vr m=0 


and in the lower half-plane 


= 27 Dd) exp [—jo(2m + 1], 


m=0 


sin yr 


(2) can be written as 


exp (jv[2(m + 1) — ])dv 


co+je 
| [see 
Il. a K 2 —o+je 


oes er f S(», $0) 


—w—je 


H,(kr) exp [jv(r — 4) 
H,® (ka) sin vr 


K 
iy = ay S(y, $0) dy (2) 
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where C, is as shown in Fig. 2. The value of this inte- 
gral may be obtained in three different ways. If the 
integral is evaluated from the residues at the poles 


6 L. L. Bailin and R. J. Spellmire, “Convergent Representation for 
the Radiation Fields from Slots in Large Circular Cylinders, Tech. 
Memo. No. 386, Hughes Aircraft Co. Res. Labs., Culver City, Calif., 
July, 1955. This report is an expanded version of the present paper 
and contains many of the details omitted here. [ } 

6 L. L. Bailin, “Radiation field produced by a slot in large circular 
cylinders,” IRE Trans., vol. AP-3, pp. 127-137; July, 1955. 2 

7G. N. Watson, “The diffraction of radio waves by the earth, 
Proc. Roy. Soc., London, vol. A., 95, pp. 83-99; 1918, and vol. A., 95, 
Se BAL ELA 1010 : 


H,® (ka) H, (ka) 
H,® (kr) A, (ka) 
H,® (ka) HH, (ka) 


(3) 
exp (—jv|2mm + ¢])dy 


where the factor H,(ka)/H,“ (ka) makes the de- 
nominators real and nonexponential for large ka. To 
obtain the geometrical optics approximation, the Han- 
kel functions in the numerators are replaced by their 
Sommerfeld integral representations: 


1 —ntjo 
H,®(z) = —{ exp liz cos 71 + jv(ri — a/2)|dr1 


T 


n—Joo 
and 
1 m—7+jo 
HON a) = exp [—jz cos r2 — jv(r2 — 1/2) |dre 
T % —r+n—Jjoo 


where —arg z<74<a—arg z. When these equations are 
inserted into (3), the result in terms of amplitude and 
phase functions is 
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Sea S(v, bo)dv 
tle H,® (ka) H,™ (ka) 
where 
Smt = — kr cos to + ka cos 71 
+ r[2(m + 1)r — o + (71 — 72)] 
and 
Sm- = — kr costs + kacosr3 + v| —2mr —¢ + (73 — 74)]. 


Only the m=0 wave is considered, since, as Bremmer? 
points out, the m0 terms correspond to waves re- 
flected from infinity and are small compared with the 
m= terms. It is shown! by a consideration of the saddle 
points in (4), that the integral with phase factor Sot 
corresponds to rays reaching the point of observation 
in the region of space where 27 >¢> 32/2 and that the 
integral with phase factor S )~ corresponds to rays 
reaching the observation point located in the region 
of space where 1/2 >@>0. If the source is symmetrical 
about ¢=0, then these two integrals should yield the 
same result for symmetrically located observation 
points. 

Thus, to obtain the desired third-order approxima- 
tion for r/2>@>0, the second term in (4) is considered 


I. = Seo S(v, bo)dv 
_ oie Hy (ka) H,™ (ka) 


—ntyeo 7—7+7joo 
. i} dr3 i drs exp [So] (5) 


n—Jo 1+—Jo 


and the function So~ is expanded in a Taylor series to 
third order in 73, T4, and v. It can be shown‘ that the 
resulting integrals may be evaluated by saddle point 
methods to yield a far field result valid for large ka or 
small ¢. Thus, (5) becomes 


exp |—7(kr—7/4 
eae Co) ee eames ee Gsariatlea 
(kr) 1/2 
and is valid where 
ka cos® @ 
3 sin? cy 
is large. Since #,=k?II, and the form factor for a half- 
wavelength delta slot is given by 
cos (1/2 sin 
S(ka sin 4, 40). = esos) 
cos? f 


(6) becomes 
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—ntjo T—7+joo 
if drs f dr2 exp [jSm*] 


n—jeo —t+4—joo 


—1+7J00 T— 9+ Ie 
f drs f dr, exp Lisa] | (4) 


n—joo —r+n—Jjo 


a i — i(kr — 1/4 
E, = B?K2(ka) |= ie ra a |= an : ‘| 


23/24 /3(Z) exp E (a cos¢ — Z+ =) 


(7) 
ZH )3(Z) His (Z) 


which is the principal component of the radiation field 
Ep at 0=a/2. 

Eq. (7) is now the third-order geometrical optics re- 
sults which should agree precisely with the second- 
order results as 60. When ¢—0, Z—~ and the po- 
tential tends to 


exp [—j(&r — 1/4)] 


Tl, = 2K 
(kr) 1/2 


exp [jka cos 4] 
(2/m) 1/2 


which leads to the second-order result! given by 


ka cos ¢ S(ka sin ¢, do) (8) 


E. = BK(2m)!"ka E (7/2 sin ] 
cos ¢ 


xp |—7(kr — 7/4 
-exp [j(ka cos dill = PI an [Ot (9) 


To check the accuracy of the third-order representa- 
tion from (7) and the second-order values given by (9), 
the far field harmonic series result for the half-wave slot 
can be obtained from (1) and the asymptotic behavior 
of H, (kr) as 


(ka cos 6)Z"/Hy5°(Z) exp] j (i Hee ae eps =) 


ZH1)3(Z)H1)3(Z) 
Howe BK (2) "Gy: cep ey eg 
1 (kr) 1/2 
a cos Mie exp [—jno] 
. ina / 2 ae 
2 exp [jnn/2] fine ee eee 


If the common factors are canceled from (7), (9), and 
(10) and if they are all renormalized so that the co- 
efficient of the series in (10) is unity, then the geometri- 
cal optics approximations are 
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ee |= (r/2 sin = 


ka cos ¢ 


5 
Z2Hy)(Z) exp E (40 conan 2 =) 


12 
ml 
ZH 1/3" (Z) H1/3°(Z) wy 
to third order. 
m |cos (1/2 sin ¢) 
E,~ hee exp [j ka cos | (12) 
ka cos } 


to second order, and of course the harmonic series be- 
comes 


: 5 egestas 
ee > exp [jnx/2 — jno| 2ka (13) 
aw H(Ra) —— [(ba)* = v4] 


which are the results tabulated by Bailin® for the prin- 
cipal H plane (@=90°). 

Fig. 3 shows a numerical comparison for ka = 12 of the 
second- and third-order geometrical optics results with 
the previously computed harmonic series for both mag- 
nitude and phase in the principal H plane. For this value 
of ka, the harmonic series has been accurately computed 
for a reference. It is evident from Fig. 3 that the third- 
order geometrical optics values are an improvement 
over the second-order values since they approximate 
the harmonic series over a larger azimuthal sector of 
space in both amplitude and phase. This improvement 
is due to the third-order terms which were taken into 
account in the evaluation of the integral in (5). 

For example, both amplitudes are about 97 per cent 
of the harmonic series amplitude at ¢=50°, but at 
= 80° the second-order approximation is 58 per cent 
while the third-order approximation is 83 per cent of the 
harmonic series. In the phase, the second-order approxi- 
mation is 98 per cent of the harmonic series at ¢=50° 
and drops to 97 per cent at ¢=80° while the third- 
order approximation stays above 99 per cent of the har- 
monic series in this range. However, both of these 
approximations begin to deteriorate badly, especially 
in amplitude, as ¢—7/2 from below. There are several 
possible reasons for this deterioration. In the third- 
order approximation, it is thought that the primary 
source of difficulty lies in the saddle point evaluation of 
the v-plane integration in (5). Here the technique (de- 
tails are given by Bailin and Spellmire®) required the 
splitting of M13“ (x) into an amplitude and phase func- 
tion. For large x, this division was accomplished with 
the aid of the asymptotic behavior of the Hankel func- 
tion. However, it can be shown that x is large if either 
ka is large or ¢ is small. In the region where ¢—7/2, x 
is small regardless of the value of ka. Consequently, a 
different method must be used in the splitting of 
H,3(x) in this azimuthal sector, generally referred to 
as the transition region. 
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Fig. 3—£, for circumferential slot ka =12. 


B. Residue Series Representation (r> | | >7/2) 


As stated in Section I, an evaluation of (2) at the 
poles, determined from H,,,°)(ka) =0, gives the residue 
series representation. When a Watson transformation 
is used on (2) in which the contour C, in Fig. 2 is dis- 
torted to the contour C2 to enclose the complex poles 
Ym, the following is obtained: 


Il,= —j2rK a SPY Sa ee ta? 14 
mine XD TPR / Al > «sll foyel Sebi” 
where 
7 SID YmT 0H, (ka) 
M41 a . | | =] 
exp [jrmm/2] Ov V—Vm 


For the far field (large 7), if the asymptotic approxi- 
mation 


H,,.(kr) 2\12 
Bae prem ape (—) exp [—j(kr — 1/4)] 


and the source function® for the half-wavelength delta 
slot are used, FE, is obtained from (14) as 


Be s()! (2) em ste 270) 


Tv 


(=) 
ae 2ka cos [¥m(x — ¢)| 


15 
M., (15) 


Ds 


Amon 
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378 IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 
TABLE I 
COEFFICIENTS FOR THE COMPUTATION OF Um 
m Aim Aom Agm Aun Aim Agm Bim Bom Bam Bam 
1 0.9278790 | 0.0573973 | 0.0182598 | 0.0029389 | 0.0012396 0.0012127 1.6071338 | 0.0994151 | 0.0050902 | 0.0021471 
2 1.6223057 | 0.1754581 | 0.0975930 | 0.0274627 | 0.0202530 0.0346404 | 2.8099161 | 0.3039027 | 0.0475667 | 0.0350793 
3 2,.1908369 | 0.3199841 | 0.2403540 | 0.0913382 | 0.0909658 0.2101110 | 3.7946412 | 0.5542293 | 0.1582025 | 0.1575575 
4 2..6933084 | 0.4835936 | 0.4465596 | 0.2086204 | 0.2554220 0.7252791 4.6649475 | 0.8376095 | 0.3613413 | 0.4424042 
5 3.1526334 | 0.6626059 | 0.7162120 | 0.3916574 | 0.5613002 1.8656465 | 5.4605218 | 1.1476682 | 0.6783707 | 0.9722010 
6 3.5806434 | 0.8547331 | 1.0493115 | 0.6517136 | 1.0607993 4.0045622 | 6.2018571 | 1.4804426 | 1.1288013 | 1.8373595 
7 3.9844482 | 1.0583875 | 1.4458583 | 0.9992755 | 1.8099590 7.6032218 | 6.9012675 | 1.8331829 1.7307961 3.1349431 
8 4, 3687383 | 1.2723905 | 1.9058524 | 1.4442309 | 2.8681905 | 13.2106676 | 7.5668775 | 2.2038472 | 2.5014816-4 4. 9678550 
9 4. 7368138 | 1.4958256 | 2.4292941 | 1.9959871 | 4.2979301 | 21.4637892 | 8.2044031 | 2.5908486 | 3.4571516 | 7.4442385 
10 5.0911066 | 1.7279561 | 3.0161831 | 2.6635522 | 6.1643701 33.0873220 | 8.8180562 | 2.9929108 | 4.6134084 | 10.6770096 
11 5.4334743 | 1.9681744 | 3.6665195 | 3.4555958 | 8.5352427 | 48.8938486 | 9.4110546 | 3.4089815 | 5.9852685 | 14.7834842 
ie 5.7653767 | 2.2159692 | 4.3803035 | 4.3804952 | 11.4806429 | 69.7838007 | 9.9859265 | 3.8381751 7.5872415 | 19.8850705 
1 6.0879863 | 2.4709030 | 5.1575349 | 5.4463702 | 15.0728680 96.7454525 | 10.5447027 | 4.2797338 9 .4333913 | 26.1070066 
14 6.4022640 | 2.7325966 | 5.9982137 | 6.6611135 | 19.3863367 | 130.8549281 | 11.0890478 | 4.7330009 | 11.5373888 | 33.5781433 
15 6.7090091 | 3.0007177 | 6.9023400 | 8.0324133 | 24.4973875 | 173.2761977 | 11.6203459 | 5.1974007 | 13.9125502 | 42.4307491 
Eq. (15) is now compared with the harmonic series 0H, (ka) 
of (13) by canceling the common factors and renormal- ap =e 


izing to give 


Qe 2 COS Vm(r — ¢) 


ka 


To evaluate the residue series of (16), numerical val- 
ues of the complex roots, v, and the associated function 


ee] 
Ov V=Vin, 


E, 


must be obtained. The complex roots of H,,,(z) =0, 
where z=ka is large, can be obtained in several ways.? 
The Sommerfield integral representation for the Hankel 
function leads to a second-order approximation for Ym, 
and the Langer’ representation for the Hankel function 
of large argument leads to a third-order approximation. 
To facilitate computation for different values of ka, 
both the second- and third-order approximations for 
Vm Were put into the following form: 


Dae Ane 1 — Ame! a Aga 8A a Ol 
Ame Ame 98 ee pie 
—j{ Bint"? + Boma"? — BymZ®/?-+ Baz t3-+ » » - } (17) 


where z=ka. 

The third-order approximations were used in prefer- 
ence to the second order, because better results were 
obtained with the same amount of effort once the co- 
efficients Anm and Brym were known. The coefficients for 
the third-order approximations are as accessible as those 
for the second order, since they depend upon the roots 
of the transcendental equation J1/3(x) +J_1/3(x) =0, 
which are tabulated.’ 

Similarly, the associated function 


_ §R.E. Langer, “On the asymptotic solutions of differential equa- 
tions with application to the Bessel functions of large complex order,” 
Trans. Amer. Math. Soc., vol. 34, pp. 447-480; 1932. 


was cast into the form: 


Pier 
Ov V=Vm, 


= Cot VEL Dine? = Dont 8 Das tl8 pol aie} 
— jE me-*!*{1 — Dim?!* + Damz8!8 + +++}. (18) 


Detailed methods for the numerical evaluation of the 
constants Anm, Bum, Cm, Dnm, and E, which are shown 
in Table I and Table II are given in a previous report.® 

To check the validity of the equations, sample compu- 
tations were made using ka=12 and the third-order 
approximation for vy, and 


0H, (ka) | 
0, v=Vm 


in (16). Four terms in the residue series were found to 
be sufficient to approximate accurately the fields in both 
magnitude and phase for about 70° on either side of 
¢ = 180°, compared with 25 terms in the harmonic series. 
The results of (16) for 1<m<4 and ¢=180°(10°)90°, ex- 
pressed as a magnitude and phase, are shown in Fig. 4 
and are compared with the harmonic series results of 
(13). The residue series curves have no noticeable devi- 
ation from the harmonic series curves for both amplitude 
and phase in the range 180° >> 100°. Since the residue 
series begins to deteriorate as 690°, check computa- 
tions for 6=90° and ¢=95° were made to indicate the 
sum of the residue series after each of fifteen terms was 
considered. The result is poor but recognizable for $ 
= 95° after fifteen terms; however, for ¢=90° the series 
did not converge sufficiently to indicate even remotely 
the result obtained from the harmonic series. If the sum 
of the series is computed using the Euler® transforma- 


® K. Knapp, “Theory and Application of Infinite Series,” Hof 
Publishing Co., New York, N. Y., 2nd Eng. ed., pp. 244-246. 
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TABLE II 
0H,® (zg) 
COEFFICIENTS FOR THE COMPUTATION OF ————— 
Oy v=vm 

Mm Ga Dig Dom Dam. y Des 

1 1.1131027 0.3711514 0.4045138 0.1269169 1.9279502 

2 —1.2748603 0.6489219 1.2365613 0.6783325 —2.2081226 

3 1.3734258 0.8763343 YP SOA) 1.6706108 2.3788431 

4 —1.4458853 1.0773228 34081819 3.1038693 —2.5043466 

5 1.5038019 1.2610527 4 6697918 4.9781220 2.6046611 

6 —1.5523556 1.4322566 6.0238305 72933722 —2.6887585 

7 1.5943383 1.5937785 7.4591089 10 .0496210 2.7614748 

8 —1.6314333 1.7474944 8 .9673195 13 .2468689 —2.8257252 

9 1.6647605 1.8947246 10.5420045 16 .8851162 2 .8834495 
10 —1.6950300 2 0364416 12.1779711 20 .9643627 —2.9358779 
11 1.7228219 2 .1733886 13 .8709372 25 .4846086 2 .9840148 
12 —1.7485417 2 .3061495 15.6172996 30 .4458543 —3.0285628 
13 1.7724935 2 4351933 17 .4139749 35 .8480994 3.0700485 
14 —1.7949326 2.5609043 19 .2582913 41 6913441 —3.1089142 
15 1.8160517 2 .6836023 21.1479054 47 .9755884 3.1454936 

/ . 
K H,® (kr) exp [jv(r — ¢)] 
I,* = oe ee ee 
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Fig. 4—E, for circumferential slot ka =12. 


tion, excellent results are obtained in the range 105°>¢ 
>90° using four terms. 


II]. AxIAL SLOT 


For the radiation characteristics of an axial slot, it is 
also possible to start with the equation’ for the mag- 
netic Hertz vector 


27 Jo H,'(ka) sin vr 


where C; is again the contour of Fig. 2. For an infinite 
axial slot or the principal H-plane fields of a half-wave- 
length slot, 


Vo 
j2mopka 


/ 


where Vo=voltage across the gap for a delta source. 
As in the case of the circumferential slot, this integral 
may be evaluated in three different ways. If the residue 
theorem is used at the poles y=n, then (19) becomes 
the harmonic representation, 


ia ky 


H,)(kr) 


H.'(ka) exp [—jng]. 

If the integral is evaluated at the poles given by 
H.,,,’(ka) =0, then the residue series representation 
(which will be discussed in Part B) is obtained. Of 
course, the third possibility is to evaluate the integral 
along the path Ci, without using the residue theorem. 
This method, again, will lead to the geometrical optics 
results. 


A. Third Order Geometrical Optics Representation 
(|| <7/2) 

To obtain the geometrical optics results, the contour 
C, is considered, and the integral is modified as before 
so that (19) can be given as 


exp [jv(2m + 1)x — ¢]dv 


: tie H,2)(kr)H,’ (kr) 
"EL 
: erie H,'(ka)H,' (ka) 


wo—je H,® (kr) H,’ (ka) 
fea H,”’ (ka) H,)' (ka) 


m=0 = 


exp [—jv(2mmr + 6)|e | (20) 
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From a consideration of the saddle points of the inte- 
grals in (20), it can be shown! that the first integral 
gives the fields in the region of space where 27 > $ > 32/2 
and that the second integral describes the fields in the 
region of space where 7/2>@>0. Thus, if again only 
the m=0 term for 7/2 >¢=>0 is considered, then 


H,® (kr) H,’ (ka) 


co—je 
_'K (. Penk ie Sa (21) 
je Hy®" (ka) H, (ha) 


exp [—jro]dy 


This integral may be evaluated to third order by the 
techniques described by Bailin and Spellmire!® and 
yields a far-field result which is valid for large ka or 
small ¢. Thus (21) becomes 


j2'K exp [—j(kr — 1/4)] ae) 
= (kr) 1/2 


cx -s( 2) sty a (3) on 


where 


Hise = 


a 
[1 + (1/2ka)?]*"2 
H,,)" (ka) Hy, (ka) 


sin y.\!/? 
Ai) 
) ( ka ) 


f(vs) = [(ka)? — (v,)2] "2 + y, sin“ (=) 


and 


iy Ee 1 13», | 
rie Wa eamuale eae 
The value of v, is determined from 


fon —aes(2) +B one 


and 
a= 8") B= Hf"). 


It can be shown!‘ that to third order 


H,,’ (ka) A,, (ka) = oe aN B?) 
where 
X1/2 y) 
iene | tag exyecetitsT 4) ee eipintx 
oe —9/3(X) we sin? YJ 1/3(X) + Wz 1/3( | 


ees (as 
GS eras) SS 29g X) SIT a ed ae A 
Peay | 2/3(X) n? yoJ1/3(X) 


2 
= Jar(X) | (23) 


/3 
with 


10 Op. cit., Appendix B. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


October 


ka sin® y 


Since the far-field components are obtained from the 
Hertz vector 


rs) Hig 
or 
Ee 4 } 218 Wig 


Eg = jou 
(24) 


then to third order as r> ” 


k= (——) (=) "(=)" exp Her ofA) 


2ra T wkr. 


“exp E eal (= ie | 


* oS oo epee ee eer “o een 


(25) 


Now as ¢—0, v.—0 and 
ay 
> 0, 
27a? 
nro pZoalle-8) 
me 2 12 


2 
H,," (ka) Hh, (ka) > — 


ar 


and 


When these asymptotic behaviors are used, (25) re- 
duces to 


sage) hs 
-exp [—j(kr — 7/4) }rbo| 1 


; 1 
“exp E (+0 cos @ + =) 


which is a better approximation in both amplitude and 
phase than the second-order result! given by 


rg (5) (=) ° exp [=i(ér — /4)] nko 


2ra/ \rkr 


‘. aia 


(26) 


-exp [j(ka cos ¢)]. 


To make a numerical comparison of the results from 
these approximate equations and of the harmonic series 
for the far field which can be obtained from (20) as 


V 2 1/2 
Ey = (-) (—) exp [—j(kr — 7/4)] 

exp Ijnr/2 — jno| 
| 2 ebay 


the common factors are canceled out, and the harmonic 
series is renormalized to give 


-exp [—jr/2 5 (27) 
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: exp linx/2 — jno] 
Ey = exp [—jr/2] - (28) 


H,," (ka) 


Then the modified second-order geometrical optics 
approximation becomes 


Ey = rha| 1 — ——. 

12(ka)? 
and the third-order geometrical optics approximation 
becomes 


Ey = (=) 400 exp Ee ie Ga f a 


3 
es (= “). 

3V3 a 21 oF 
Numerical results of (28) (harmonic series)", (29) 
(modified second-order geometrical optics approxima- 
tion), (30) (third-order geometrical optics approxima- 
tion), and Sensiper’s second-order approximation are 
shown in Fig. 5 for ka=12 and in Fig. 6 for ka=40. 
The gain of accuracy of the third-order geometrical 
optics approximations over the second-order geometrical 
optics approximations is negligible for the phase; how- 
ever, the third-order approximations give better values 
for ¢>0°, since the amplitude curves for the second- 


(30) 


1 The numerical results for the harmonic series for ka =12 were 
taken from Bailin, op. cit., and those for ka =40 from computations 
by the NBS Computation Lab. presented to us in a private communi- 
cation from Dr. W. Saunders of the Diamond Ordinance Fuze Lab- 


oratory. 
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Fig. 6—£g for axial slot ka =40. 


order approximations are constant (that is, independent 
of @) while the amplitude curves for the third-order ap- 
proximation tend to follow the harmonic series curves. 


B. Residue Series Representation 


When the method of the circumferential slot is fol- 
lowed (19) becomes 
eo 4H ge) (kr) 


— 72n’'K 
sho p> exp [j’tmx/2] 


COS '¥m(a — o) 
My 


where 


j sin Umm 


8H,’ (ka) 
M rym \ = | . 


exp [j’¥mm/2| ov 
Since 
orl,* 


) 


Es = jwu 
tous or 


then in the far field (,( ~), E, becomes 


V /2\32 
ae —(=) ey oe weer 


2mra \rkr 


2a ‘ 32 

eT MY onl 3) 
For comparison purposes, the common factors of (27) 
and (32) were canceled, and (32) now becomes 


COS '¥m(ar — ) 


33 
M| vn re 


Ey ~ 2x exp [jr] >> 


m=1 
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TABLE III 
"Wm 
COEFFICIENTS FOR COMPUTATION OF ie 
a 
M Kom Kim Kam Kom Kam Kem 
1 1 | 0.40430823 (1—j./3) | 0.07273174 (—1—jx/3) | 0.06525312 | 0.01188593 (—1+j/3) | 0.00105769 (—1—jv/3) 
2 1 1.28904796 (1-7/3) | 0.13017057 (—1—j/3) 0.10353615 0.04777835 (—1+jV/3) 0.05356433 (=1—7\/3) 
3 1 | 1,91285748 (1—x/3) | 0.25700457 (—1—j/3) | 0.21435659 | 0.10779227 (—1+7,/3) | 0.15174500 (—1—j/3) 
4 1 2.44590997 (1—j4/3) | 0.40905317 (—1—j+/3) 0.38878826 0.21188267 (—1+jv/3) 0.34684651 (—1—jv3) 
5 1 2.92565024 (1—j»/3) | 0.57917419 (—1—j/3) 0.62669711 0.37421154 (—1+ 4/3) 0.63461976 (—1—7j/3) 
6 1 3.36865792 (174/83) | 0.76394571 (—1—j+/3) 0.92806346 0.60700924 C—-1-F4/3) 1.11126040(—1 —jv/3) 
7 1 | 378414515 (1-3/3) | 0.96125758 (—1—74/3) | 1.29288129 | 0.92138696 (—1+jx/3) | 1.81438363 (—1—J4/3) 
8 1 4.17790446 (1-7/3) | 1.16964384 (—1—j+/3) 1.72114855 1.32768068 (—1+j V3) 2 .80118528 (1 ==fa/a) 
9 1 455387890 (1-7/3) | 1.38801180 (—1—j/3) 221286439 1.83563957 (—1+jV/3) 4.13218591 (—1—ja/3) 
10 1 | 4.91490612 (1-7/3) | 1.61550800 (—1—jx/3) | 2.76802826 | 2.45454479 (—1+j/3) | 5.87087372 (—1—jv/2) 
11 1 5 26311466 (1—Jx/3) | 1.85144325 (—1—jx/3) | 3.38664078 | 3.19329156 (—1+j/3) | 8.08345789 (—1—jx/3) 
12 1 5 60015281 (1—74/3) | 2.09524661 (—1—7/3) 4.06869950 4 06044434 (—1+jv/3) 10.83502643 (1 74/9) 
13 1 592733011 (1—jr/3) | 2.34643577 (—1—j/3) 4.81420662 5 .06428846 (—1+jV/3) 14.20739296 (—1—J4/8) 
14 1 6.24570886 (1—jx/3) | 2.60459678 (—1—j1/3) 562316141 6.21285803 (—1+jv3) 18.26292494 (--Tjy/a) 
15 i 6.55616578 (1—jx/3) | 2.86939695 (—1—7+/3) 6.49556373 7.51396679 (—1+ 74/3) | 23.08049493 (—1—j+/3) 
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for ‘v.12 As before, ‘vm was cast in the form 


le 
Vin 


= > Kym(1/ka)??!? 


p=0 


ka Se 


where the coefficients K,, are given in Table III. The 
associated function 


0H, '(z) | 
Ov V='Vin 


12 Details of the numerical techniques involved are shown by 
Bailin and Spellmire, op. cit. 


and f,,,,(z) is obtained from the Langer representa- 
tion 


tan am — 'a&m\1!/2 
H,,,2(2) = (a) 
tan Qm 


“exp [—jr/6]H1/3 [yn (tan ‘Om = ‘am) | (37) 


where '¥,,=2) COS ‘Atm. 

Results of (33) and (28)! are shown in Fig. 7 for ka 
= 12 and in Fig. 8 for ka=40. For ka=12 up to fifteen 
terms (m= 15) were used, but for ka = 40 the use of only 


a 


nea 


om ge 


—s 
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four terms was considered adequate. The residue series 
curves have no noticeable deviation from the harmonic 
series curves in both amplitude and phase in the range 
180°>@>100°. In the range 100°>¢>90° the residue 
series curves for amplitude deviate quite noticeably 
from the harmonic series curves, while the residue series 
curves for phase are well within graphical accuracy. 
If, however, the Euler summation technique is used on 
the residue series in the region around ¢=90°, good 
results are obtained. It is worth noting that the residue 
series can be summed, in the Euler sense, in the geomet- 
ric optics region 90°>¢>80° using eight terms. 


IV. CONCLUSION 


The approximation formulas described present a 
tractable, and in some cases simple, representation for 
the radiation fields from slots on cylinders. The accuracy 
of the various approximations discussed varies with the 


‘type of slot and with the region of space which is of 


interest. In the geometrical optics region, for the cir- 


-cumferential slot, the third-order values are an im- 


provement over the second-order values, since they 


_ approximate better the true field over a larger azimuthal 


sector of space in both amplitude and phase. For the 
axial slot, the third-order approximation is the lowest 


_ order that can be used to obtain the magnitude, unless 


one is willing to accept a constant (that is, independent 


of @). 
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In the shadow region, where the residue series is 
applicable for both the circumferential and axial slots, 
it is easier to compute from the third-order approxima- 
tion formulas than from the second-order formulas if 
the tables in this paper are used. Unfortunately, no 
computations using second-order approximation for- 
mulas were made to test the gain in accuracy. 

In general, the third-order approximations are ade- 
quate to approximate the fields except in the transition 
region (that is, the region around 90°). In this region, 
the approximations involved in the geometrical optics 
formulas are no longer valid, and the residue series re- 
quires too many terms to be practical when ordinary 
summation is used. However, the Euler transformation 
gives good results with the residue series even into the 
geometric optics region around ¢=90°. Further study 
is being given to extending the range of the residue 
series by using a generalized Euler transformation. !° 
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Some Observations on Scattering by Turbulent Inhomogeneities* 
MARTIN BALSERt 


Summary—Several topics in the theory of radio wave scattering 


by dielectric inhomogeneities due to turbulence are discussed. A 


corrected derivation of the scattering formula is presented (which 
does not change the result), using the proper dyadic Green’s function 
and examining the terms in the wave equation caused by the non- 
vanishing divergence of the scattered field. The statistics (distribu- 
tion and correlation functions) of the received signal are considered. 
In particular, the significance of the space correlation function in 


3 testing some theoretical results is questioned. It is suggested that 


the frequency dependence of the time correlation function may give 


- more information about the nature of the turbulence. 


INTRODUCTION 
ete process of scattering radio waves by di- 


electric inhomogeneities has been of interest for 
some time in connection with the anomalous long- 
distance propagation of vhf and uhf (and higher fre- 
* Manuscript received by the PGAP, July 25, 1956; revised 


i i in this docu- 
manuscript received, January 29, 1957. The research in t 
ment was supported jointly by the Army, Navy, and Air Force under 


contract with the Massachusetts Institute of Technology. 


+ Lincoln Lab., Mass. Inst. Tech., Cambridge, Mass. 


quency) waves. Despite this, a number of points about 
the electromagnetic theory involved seem to have 
never been adequately considered. It is the purpose of 
this paper to examine critically some of the assumptions 
and steps in the development of the formalism of scat- 
tered fields, as well as to discuss some inferences pertain- 
ing to these formulas that may be made from elementary 
observations on the nature of turbulence. 


DERIVATION OF THE SCATTERED FIELD 


The Quasistatic Approximation 


The wave equation for the electric field, E(r), is 
easily derived from the pair of Maxwell’s equations 
connecting the electric and magnetic fields, and may be 
written 


os 


— po—- 


ot 


e2 
Marcie hiner kek) = (1) 
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It is already assumed that the magnetic susceptibility 
has its free-space value, mo, everywhere in the region 
to which (1) applies. 

The basic assumption we shall make (in various 
forms) several times in the derivation is what we may 
call the quasistatic approximation; that is, since one 
variation takes place very slowly compared with an- 
other, we regard the problem as though the first varia- 
tion did not exist and simply add its effect to the solu- 
tion of the “static” problem. The justification for such 
an approximation must, of course, be drawn from ob- 
servations (whether theoretical or experimental) of the 
physical situation. 

In this case, the approximation is actually applied 
twice. Firstly, since the highest frequency with signi- 
ficant power in the spectrum of e(t) is very small com- 
pared with w, the radio frequency, we approximate 

0? ri 
— (eH) >~e—-: 


0? ot? (2a) 


Further, we assume that there are no impressed cur- 
rents in the region, so 
os OL 


ee fi 


ot Ot 


(2b) 


The “static” (or monochromatic) problem may now 
be solved by replacing 0/d¢ by —ww (this convention 
will be used throughout), so that (1) may now be 
written 


VXVX £E— warek = 0 (3) 


where by € we now mean the effective (generally com- 
plex) dielectric constant e+i0/w. 

There is abundant theoretical and experimental evi- 
dence that dielectric configurations in the atmosphere 
do not change in the order of microseconds or less, so 
this first assumption is easily accepted. At this point, 
however, a second assumption is made. We are ordinar- 
ily interested in the statistical properties of the scattered 
field, and hence in the statistics of the turbulent field 
which causes it. These properties are studied by taking 
time averages, such as the average power, distribution 
of signal level, and correlation functions. The assump- 
tion made is that the statistics themselves are not chang- 
ing while the measurement is being made, or in other 
words, that the process is what we might call “locally 
stationary.” Thus we may solve the statistical problem 
as though the process were truly stationary, then ac- 
count for the slowly changing statistics afterward. As 
an example, we note that the received signal level is 
expected to vary (cf. our first assumption) according (as 
will be shown) to a Rayleigh distribution about a fixed 
mean determined by (Ae?). (Ae?) must, however, remain 
constant during the measurement or the observed dis- 
tribution will be changed. Itis, in fact, observed that for 
periods of the order of a minute the distribution of signal 
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level is Rayleigh, but over many minutes it is not.'? 
Thus we may expect that our assumption of “local 
stationarity” may apply as long as the time required 
for our statistical analysis is not more than a minute 
or two. (In cases of doubt, more involved statistical © 
tests may have to be applied to determine, for instance, — 
the confidence we may place in the hypothesis that the | 
mean of the process has not changed during our meas- | 
urement.) . 


Discussion of the Wave Equation 


Letting ko?=wuoeo and e(7) =e+Ae(7), we can re- 
write (3) as 


Ae 
—-VXVXE+ 27°F = — kh? — £. 


(4) 
€0 ; 
When we expand the left side of (4), we get a term 
V(V-E) which, we must note, does not vanish for the 
very reason that we get scattering; 7.¢., € is not con- 
stant. Thus, 


A i% 
2k AB Bie ee v(— veB), (5) 


€0 € 


The last term in (5) has been considered in a recent 
paper® and neglected on the grounds that | Vel ~he/Lo, 
where Lp is the “scale length” of the turbulence; whereas 
the term with which it is compared is koAe=27Ae/X. 
It is, however, not these gross characteristics of the 
turbulence which are relevant, but rather particular 
Fourier amplitudes [see, for instance, (17)]. These 
Fourier amplitudes, as we shall see, are of the same 
order of magnitude for the two terms, but for an en- 
tirely different reason the contribution of the additional 
term to the first-order scattered field is not small, but 
precisely zero. 

We wish now to solve the vector (5) and at the same 
time to satisfy certain auxiliary conditions on E.To 
do this, it is not adequate to use a scalar Green’s func- 
tion. We require rather a dyadic Green’s function; 
1.e., a solution of 


(Vig'=b ey?) I(r, or )eaied (ete tpl) (6) 


where e€ is the unit dyadic, with the property that ¢-A 
=A-e=A for an arbitrary vector, A. The standard 
solution to (6) is 


erkol r—r'| 


1 
T(r, f= (eS Say) See ee 
cry=(s ke Vesa 


1 J. H, Chisholm, P. A. Portmann, J. T. deBettencourt, and J. F. 
Roche, “Investigations of angular scattering and multipath proper- 
ties of tropospheric propagation of short radio waves beyond the 
horizon,” Proc. IRE, vol. 43, pp. 1317-1335; October, 1955. 

? R. A. Silverman and M. Balser, “Statistics of electromagnetic 
radiation scattered by a turbulent medium,” Phys. Rev., vol. 96, pp. 
560-563; November 1, 1954. 

8 A. D. Wheelon, “Near-field corrections to line-of-sight propa- 
gation,” Proc. IRE, vol. 43, pp. 1459-1466; October, 1955. 
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which automatically gives, at large distances from the 
sources, transverse outgoing spherical waves. The temp- 
tation is now simply to write the usual solution for the 
scattered field 


E;.(r) = E(r) Fay Eixnc(r) 
= fre, 7) be Ac(r’) E(r') 


€0 


0(2 veey-nie) Jo 


where the integral is taken over all space, but in effect 
only where Ae#0. This solution clearly satisfies (5), 
but is not the solution we seek. The reason is that the 
solution for I’, (7), has the additional property V-T'=0, 
which it imparts to the solution for E,,. As we have seen, 
V-E,-~0. We may now either seek another Green’s 
function which gives the proper divergence or, more 
simply, use the same Green’s function to solve the wave 
equation for D, which is divergenceless. This latter 
equation can easily be seen to be 


(V2 + ko?) D 


Ae rie 
_ c[ tot SB + v(—ve-z) 
€ 


€0 


1 i 
a 2(v —-v) D+ pvt—| 
€ € 


= — F(r) (8) 


and the solution of (8) is 


Dey) if T(r, 1')-F(r’) dr’. (9) 


Finally, in still another quasistatic approximation, . 


we may separate Ae, hence F(r), into the sum of two 
terms. The first describes the deviations of ¢ from €o 
which are fairly constant in time, in effect the average 
over, say, many minutes or an hour. It thus includes 
the normal lapse rate of air density in the troposphere 
or electron density in the ionosphere, as well as meteor- 
ological effects such as persistent elevated layers and 
ducts. The second term, which is the one we consider, 
describes the short term fluctuations in e, which we 
‘ascribe to turbulence. The solution of (9) is just the sum 
of the individual solutions. Since turbulence theories 
imply that the characteristics of the turbulent flow 
do not depend on the boundary conditions, we may 
conclude that these contributions will be incoherent. 


The Far-Field and Single-Scattering Approximations 


We may immediately simplify (9) by substituting 
the far-field approximation 
etkok1 


47 Ry 


T(r, 1’) ~ (e — aa) (10) 


where Ri=|r—r’| and @=(r—r’)/R:. (a is the unit vec- 
tor from the source to the field point.) Eq. (10) neglects 
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terms 0(\/R) compared with the term retained. As 
already noted, the form of I insures that the far-field 
contribution from each scatterer is transverse. Thus 
(8) becomes 


etkoR1 


An R, 


Dlr) ~f (e — at): F(r’) dr’. (11) 


The principal difficulty in the mathematical analysis 
is the solution of this integral equation. Fortunately 
we do not have to solve (11) in full generality. Rather, 
we use the fact that F(r)ade(r)<e) and assume that 
the scattered field may be adequately represented by 


Dy(r) 4 pitt r')-Fo(r')dr’ (12) 


where Fo(r) is the function F(r) in which we neglect 
terms 0(Ae?). In particular, D and E are replaced by 
the incident fields. This is the single-scattering approxi- 
mation. If desired, we may iterate the process to find 
the result of 2 scatterings, 


D,(r) = nt r')- Fa—s(r')dr’ (13) 


where F,,, comprises terms 0(Ae”) and is determined 
by the previous solutions. 

It should be emphasized that the adequacy of the 
single-scattering approximation is an assumption and 
does not follow from the statement that Aeé. For 
ionospheric scattering, it appears quite clearly justified 
since the inhomogeneities involved are fairly well re- 
stricted to the height of the common volume of the two 
antennas. Thus, radiation can arrive at the receiver 
only from this region. In the troposphere the case is 
not so clear-cut, but rough considerations support the 
adequacy of the approximation. A more detailed study 
of the relative importance of single and multiple scatter- 
ing is being carried out.4 We shall henceforth assume 
that (12) is an adequate approximation in the cases of 
interest. 

Since the incident fields we shall be interested in are 
due to a transmitting antenna, Fo(r) will be of the form 


etkoRo 


Ro ie) 


Fo(r) = fo(r) 


where Ro is the distance of the scattering point from 
the antenna. The space dependence of the incident 
field, except for the (relatively) slow variations as the 
direction from the radiating antenna changes, is thus 
explicitly exhibited. It is to be understood that the 
operator V in fo(r), when operating on the incident 
field, is to be replaced by iko, where ko is the propaga- 
tion vector of magnitude & and direction perpendicular 
to the wave front (here a sphere). Let us pick an arbi- 


4H. Feshbach, private communication. 
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trary origin, 0, in the scattering region (see Fig. 1) 
and rewrite 


ko(Ro + Ri) = ko(r)-ro + Ki(r) 71 + K(r)-r 


where ky and k; are the propagation vectors from the 
transmitter, JT, to the scattering point, P, and from P 
to the receiver, R. K(r) = ko(r) — ki(r) and | K (r)| =2Ro 
sin 6/2, where @ is the scattering angle. Combining (15) 
with the far-field and single-scattering approximations, 
we find 


D, = fe — at) -folr’) 


(15) 


eilko(r’) ror k1(7’) + rt] 
eiK(r) “1dr! (16) 
ArRoR1 

For the case of a very limited scattering region, far 
removed from both transmitter and receiver, the propa- 
gation vectors can be considered constant over the 
region; 7.e., the rays are parallel. This is the plane 
wave, or Fraunhofer, limit. Eq. (16) then takes on the 

simple form 
ei(ko-ro+k1* 71) 


Di = (e — na)- f folr’)e®r'dr’. (17) 


Anr 071 

It can now be shown (see Appendix) that the contri- 

butions to the first-order scattered field of all but the 

first term in F(r) [see (8)] are strictly zero. Using this 
result, we may finally write 


e* (ko -rotki +71) 


k 2 
Ey = (e — ad) By [ Ae(r’)e***"'dr' (18) 
€0 


Anrory 
where Ep is defined by 


etkoRo 


yy (19) 


Wetae = Eo 
0 

(e—ada)-Eo is, of course, that component of Eo per- 
pendicular to k;. Eq. (18) is the familiar result that the 
scattered field is proportional to the component of the 
Fourier transform of Ae corresponding to the wave 
number K. The polarization of the incident field is pre- 
served (to first order) since the unpolarized term gives 
no contribution. 

The approximation to Fraunhofer scattering is appli- 
able as long as the length L of the wavefront over the 
scattering region (see Fig. 2) as seen by either the trans- 
mitter or receiver is plane, say to within a quarter 
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wavelength. This implies (since y is a small angle) 


For a reasonable ionospheric circuit, \=6 m and d= 800 
km, and for a tropospheric circuit we may have A=1 m 
and d=150 km. Both cases require y~0.25° or less. 
Present antenna patterns are not so narrow. Taking 
W~/a, where a is the relevant dimension of the an- 
tenna, (20) becomes 


(21) 


Thus for the 150-km troposphere case with a 20-m 
dish, we would need \<5 mm to satisfy the Fraunhofer 
condition. 

In general then, an accurate solution requires con- 
sideration of Fresnel scattering, characterized by non- 
plane wavefronts or, in another way, by differing 
angles of arrival. Let us suppose that the transmitting 
antenna causing the incident field radiates Pr watts 


with a vector voltage gain function gr(r), which has the 


magnitude of the voltage gain and the direction of the 
electric vector in the far field. gr(r) is actually a func- 
tion of two angular variables only and is taken to be 
normalized; 2.e., 


1 T Qr 
— 0 
ap i | er(0, 4) 


If we now assume that the receiving antenna has a 
voltage gain function gr(7) and a radiation resistance R, 
then it can easily be shown using (16) or (18) that the 
current in a load matched to the receiving antenna is 


I ae, Ae(r’) gr(r’)-gr(r’) 
dor 2R €0 RoR, 


et (Ok -ro+1 “T1) etK (rk’) “1 dy! 


2 sin odddd = 1. (22) 


(23) 


The scalar product of the gain functions picks out just 
that component of the effective incident field [see (18) ] 
which is accepted by the receiving antenna. 

In a final application of the quasistatic approximation, 
we note that the quantities in the integrand of (23), 
with the exception of Ae and exp [iK-r], vary only 
slightly in distances which are large compared with the 


(20) 
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correlation distance of Ae. The scattering volume may 
thus be broken up into smaller volumes, many cor- 
relation distances on a side in each of which the Fraun- 
hofer approximation, (18), applies. We may, in fact, 
completely neglect the term exp {i[ko-rotki-ri]} 
since the phases in different subvolumes are uncorre- 
lated in any case. The contributions of the subvolumes 
to the received current are, of course, incoherent. 


STATISTICAL CONSIDERATIONS 
Distribution of Received Signal 


We have assumed that the scattering is taking place 
in a region which is large compared to the largest 
eddies, which are roughly regions with dimensions the 
size of a correlation distance of Ae. (These largest eddies 
are the “blobs” of heuristic turbulence theories.) We 
have also seen, at the end of the last section, that the 
other factors in the received field (or current) change 
slowly over many correlation distances, insuring that 
there are many eddies in the same region having the 
same characteristics. It is further reasonable to conclude 
from turbulence theories® that, since the motions well 
within the large eddies do not depend on the boundary 
conditions, the states of similar largest eddies, though 
statistically identical, are statistically independent. 

The received signal may thus be thought of as the 
sum of a large number of independent electrical signals 
(each corresponding to a different part of the scattering 
volume) with the same statistics. If we consider the 
(mn)-dimensional random variables consisting of the 

values of the independent signals at m different points 

at each of m arbitrary times, the generalized Central 
Limit Theorem® states that the sum (i.e., the received 
signal at those points and times) is asymptotically 
normally distributed. Stated differently, the received 
signal is a Gaussian process in space and time. (The 
process is, in fact, a narrow-band Gaussian noise, since 
the variations in amplitude and phase are slow com- 
pared to the rf.) 

If the scattering region has several subvolumes differ- 
ing in some characteristic such as intensity of turbulence 
or antenna pattern, the signal from each subvolume is 
Gaussian, as above. The total signal is therefore also 
Gaussian and its spectrum is the sum of the individual 
spectra. The characteristics of the received signal are 
thus completely determined by the mean power and 
correlation functions (or spectra) in space and time. 

As already noted, an immediate implication of this 
conclusion is that the one-dimensional distribution of 
the signal is Gaussian and that of its envelope is Ray- 
leigh. Considerably more complex statistical tests would 
entail finding joint distributions at a number of differ- 
ent times and places, all of which must conform to the 


5 See for instance C. F. v. Weizsicker, “Das spectrum der tur- 
bulenz bei grossen reynoldsschen zahlen,” Z. Phys., vol. 124, pp. 


614-627; 1948. oy" : 
6H. Cramér, “Mathematical Methods of Statistics,” Princeton 


University Press, Princeton, New Jersey, p. 316; 1946. 


Balser: Some Observations on Scattering by Turbulent Inhomogeneities 


387 


Gaussian character of the signal. A note of caution may 
be required in that the signal is not strictly a Gaussian 
process for any finite number of components, and a 
sufficiently fine analysis of the received signal will pro- 
duce some departure from the Gaussian character. 


Mean Received Power 


Using the approximations already stated, and defining 
a (normalized) correlation function of Ae by 


(Me(7!)Ac(r’’)) = (Ae(r’)) Rr’ — 77), (24) 
we find from (22) that the mean received power 
p, - MIEPR 
2 
_ ke Pr ¢ (gr(r)-gr(r))? (Ac?(r)) 
ab oar nn a Bion S(K(r))dr (25) 
where 
S(Ki= f Raper au, (26) 


S(K) is thus the component of the spectrum of Ae cor- 
responding to the wave number K (which may be a 
function of position). These expressions could be some- 
what generalized by including a dependence on position 
of the correlation function (hence the spectrum). It is 
however, precisely the conclusion of the dimensional 
theories of tubulence that the spectrum is a universal 
function in the inertial range to which the theory ap- 
plies. The only dependence on position to be expected in 
the correlation function of Ae is the normalizing param- 
eter, related to the intensity of turbulence, shown 
explicitly in (24). 

The expression for Pz exhibits the contributions from 
different regions and, in particular, the factors which 
effectively limit the “region of integration;” 7.e., the 
region from which most of the power is scattered into 
the reciever. The first term in the integrand gives the 
common region of illumination of the two antennas 
(together with the attenuation due to the expansion of 
the wavefronts with distance). The case where this 
factor limits the region of integration is treated by 
Booker and de Bettencourt? and implicitly by Villars 
and Weisskopf* and other Fraunhofer treatments. The 
last term gives the decrease of scattered power with 
increasing angle of scattering (recall that K is deter- 
mined geometrically by the angle of scattering). This 
may be the limiting factor if the antenna beam is fairly 
wide. Finally, (Ae?(r)) may be effectively nonzero only 
in a limited region of space; e.g., a layer of turbulence. 
The extent to which the region of integration is deter- 


7H. G. Booker and J. T. deBettencourt, “Theory of radio trans- 
mission by tropospheric scattering using very narrow beams,” 
Proc. IRE, vol. 43, pp. 281-290; March, 1955. 

8 F, Villars and V. F. Weisskopf, “The scattering of electromag- 
netic waves by turbulent atmospheric fluctuations,” Phys. Rev., vol. 
94, pp. 232-240; April 15, 1954. 
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mined by these three (or four) factors must, of course, 
be determined in each case. Obviously, if one factor 
strongly predominates, the others may be taken to be 
constant, thus simplifying (25) still further. 

The scattering process can, alternatively, be de- 
scribed by the mean scattering cross section per unit 
volume; %.¢e., the average power scattered per unit 
volume into a unit solid angle when the incident field 
has unit power per unit area (thus the dimensions are 
those of inverse length). This formulation has the ad- 
vantage of depending only on the characteristics of the 
scattering region, and not on the antenna patterns, 
transmitted power, etc. (The received power is then de- 
termined by the incident power, scattering volume, and 
receiving antenna gain, so that these quantities must, of 
course, be known to extract the scattering cross section 
from any experimental data.) Note that the introduc- 
tion of a quantity which, when multiplied by the rele- 
vant volume, gives the total power scattered in a given 
direction presumes that the process is incoherent. 

From (25) we can see that the mean cross section per 
unit volume is 

4 2 
dyiy a ME see) 


167? €0” 


MAS eM) cK). 
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Finally, in a volume V, in which (Ae?(r)) and K(r) are 

effectively constant, (27) may be rewritten for any 
point in V as 

(K) w 1 

Q 4 Veo 


> {| Ae(X) |?) (28) 


where 
Ae(K) =f Actryetrar. 
Vv 


V in this case need not be identified with any physically 
significant volume. It may be picked arbitrarily, as long 
as its dimensions are large compared with a correlation 
distance of Ae, yet small enough for the indicated con- 
ditions. 


Correlations of the Received Signal 


For a given (fixed) transmitter, we note that the 
only quantities in (23) which are at our disposal to vary 
are the time in Ae(7, ¢) and the vector K(r). The general 
(normalized) correlation function between the currents 
resulting from two such situations is 


psedleg 
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where AK = Ki—K>. The assumption made is that the 
change between the two conditions is so slight that h 
and J2 have the same mean amplitude, and the scatter- 
ing volume is substantially the same for both cases. 
(Insofar as the volumes are different, those contribu- 
tions are clearly uncorrelated, but the effect is negligible 
for practical systems.) 

The expression in (29) (actually only the real part 
of it) is the rf correlation of the received signals. It may 
be related to the correlations ordinarily measured; 7.e., 
those of amplitude and phase.? In particular, 


A - f Lanti ) 
canis CENT Ret Ur ee 


is the expression for the amplitude correlation function 
Tats 

If the correlation is simultaneous, =f, and (29) 
may be rewritten 


gr(r): gn(r))? 


pa f SADE (acl) S(Kle)) eH (31) 


(30) 


In particular, for frequency correlation 


Aw 
AK(r) = — K(r) (32) 
® 
and for space correlation it can easily be seen that 
Ko 
AK(r) es AK\(r) ie R, aX (a x s) (33) 
1 


where s is the vector displacement of the second receiver 
from the first. 

In effect then, these correlation functions are Fourier 
transforms of the distribution of power contributed 
by different parts of the scattering volume. Equiva- 
lently, they may be thought of in terms of an “angular 
spectrum.”!° In any case, the significant feature of 
(31) is that only insofar as the scattering volume is con- 
trolled by the falloff of the spectrum is the correlation 
function any indication of the characteristics of the 
turbulence. Insofar as the other factors affect this vol- 
ume (for instance, by using high-gain antennas), the 
correlation function is determined by geometrical con- 
siderations. 

In evaluating (31) (even assuming the scattering 
volume is determined completely by the spectrum), 
one cannot expect more than order-of-magnitude agree- 
ment with experimental correlation functions unless one 
uses much greater accuracy than is ordinarily done. Al- 


® R. Price, “A note on the envelope and phase-modulated com- 
ponents of narrow band Gaussian noise,” IRE Trans, vol. IT-1, pp. 
9-13; September, 1955. 

0 H. G, Booker and P. C. Clemmow, “The concept of an angular 
spectrum of plane waves, and its relation to that of polar diagram 
and aperture distribution,” Proc. IEE, pt. III, vol. 97, pp. 11-17; 
January, 1950. 
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so, the spectrum S(K) is itself still a matter of much 
conjecture, so a verification of a hypothesis which de- 
pends to some extent on a particular form of S(K) must, 
at least to that extent, be considered uncertain. A case 
in point is a recent attempt!! to show that experimental 
space correlation curves, which are broader than ex- 
pected, could be explained by assuming that atmos- 
pheric turbulence is anisotropic. In that derivation, a 
more refined volume integration would have produced 
a broader correlation curve without the assumption 
of anisotropy. Some effect would also be obtained by 
taking a spectrum which falls off more rapidly than the 
assumed one. It is also pertinent to note that one of the 
fundamental observations of the statistical turbulence 
theories! is that over distances much smaller than the 
dimension of the largest eddies, the dynamic laws are 
universal and isotropic despite the particular large- 
scale anisotropic source of the turbulence. (A reserva- 
_ tion is required for the ionosphere, where the magnetic 
field continues to affect the dynamics of the flow of 
electrons on all scales and may well produce anisotropy.) 

A test of theory is desirable which does not depend so 
critically on exact evaluation of complicated integrals 
or dubious approximations of antenna patterns. Such 
tests as frequency or angle dependence of received 
power have been proposed and are suitable, but difficul- 
ties of one sort or another have limited the accuracy 
with which such results can be stated. Another test 
which has thus far been largely neglected holds some 
promise not only of giving some insight into the mecha- 
nism of turbulence, but of distinguishing scattering by 
turbulence-induced fluctuations from all other pro- 
posed mechanisms for producing the fields we are at- 
tempting to explain. This test involves the frequency 
dependence of the time correlation function. 

In analogy with (24) we might write 


(Ae(r’, ts) Ae(r’’, t2)) = (Ae*(r’)) Rr’ — r’’, ty — te). 


We should note, however, that this function is not the 
correlation function of the turbulent field since we have 
neglected the drift velocity of the entire field. The rela- 
tion between the fixed coordinates and the moving co- 
ordinate system of the turbulent field is, of course, 
17=tm+Vt, where V is the drift velocity, so that the 
proper analogy with (24) is 


(Ac(1’, t1)Ae(r’’, tz)) 
= (Ae(r’))R(r’ — 7” — V(ti — te), ti — te) (34) 


where the function R(Ar, A?) is the correlation function 
of Ae in space and time in the moving coordinate sys- 
tem of the field and is the one normally discussed in 
papers on turbulence. 


1 H, Staras, “Forward scattering of radio waves by anisotropic 
turbulence,” Proc. IRE, vol. 43, pp. 1374-1380; October, 1955. 4 

2 G, K, Batchelor, “The Theory of Homogeneous Turbulence, 
Cambridge University Press, Cambridge, Eng., p. 110; 1953. 
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It is clear from the geometry [or from an expression 
such as (9) | that a receiver moving with an appropriate 
velocity (2V if the reciever and transmitter are equi- 
distant from the scattering volume) has approximately 
the same effect as if there were no drift in the turbulent 
field. Thus the time correlation function at a fixed loca- 
tion has two influencing factors. First, the space pat- 
tern is being swept by at the drift velocity," causing, 
among other things, the maxima of cross correlations 
between spaced antennas to be displaced from the 
origin. The effect of this motion on the autocorrelation 
function, as well as the effect of all models in which 
waves of all frequencies are reflected or scattered from 
the same obstacles, is that of a Doppler effect and hence 
proportional to frequency. The other influence is the 
change actually occurring in the turbulent field, exhib- 
ited by (34). The crux of the matter here is that the 
“correlation time” (or, equivalently, the fading rate, de- 
fined either as the reciprocal of a “correlation time” 
or in terms of mean crossings) for a given frequency is 
determined by the length of time required for inhomo- 
geneities involved in scattering that frequency to change 
their configuration. The fading at any frequency can 
still be thought of as a Doppler effect from the moving 
eddies, but the eddies of different sizes involved in scat- 
tering different frequencies may move at different 
speeds.“ In particular, the von Weizsicker-Heisenberg 
theory concludes that eddies of length Z move with 
speeds proportional to Z1/%. Thus we can expect that 
the fading rate would be proportional to the 3? power 
of the frequency. (This result is the same whether we 
attribute the fading to changes in position due to the 
random motions or to the actual decay in the turbulent 
field. The dimensional theory states that the lifetime 
of an eddy is the time required to move its own length.) 

The observed fading rate may thus be expected, if this 
picture is correct, to depend on a power of frequency 
somewhere between 2 and 1, closer to the latter if the 
“correlation time” due to the drift velocity is mych less 
than that due to the turbulence, close to the former if 
the reverse is true. Actually, the frequency dependence 
is somewhat clouded by the fact that the eddies involved 
are those of the effective wavelength 1/|K|. Since the 
off-center scattering effectively takes place at a higher 
frequency, the observed correlation function is a 
weighted average of correlation functions corresponding 
to slightly different frequencies. 

Two recent papers-6 have presented data on fading 
rates at different frequencies in tropospheric scatter 


18 This is the model of H. G. Booker, J. A. Ratcliffe, and IDL iB! 
Shinn, “Diffraction from an irregular screen with applications to 
ionospheric problems,” Phil. Trans., vol. 242, pp. 579-607; September, 
1950 


4 This observation has been made independently by A. D. Whee- 
lon (private communication). , 
1 L. G. Trolese, “Characteristics of epee scattered fields,” 
: . 43, pp. 1300-1305; October, : 
Tre KA noe Se L. Rice, H. B. Janes, and A. P. Barsis, “The 
rate of fading through a turbulent atmosphere,” Proc. IRE, vol. 43, 
pp. 1341-1353; October, 1955. 
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experiments. In each case the author attempted to fit 
the data to the linear dependence predicted by the 
simple Doppler picture but found that the increase in 
fading rate appeared to be too slow, as expected from 
the model just presented. 

To verify this hypothesis, experiments should be 
undertaken at widely spaced frequencies (consistent 
with the requirement that the eddies involved be in 
the inertial turbulence range) so that a power law may 
be more critically tested. The records should be taken 
simultaneously since the fading raté changes markedly 
in time’ (although in the ionosphere the fading rate 
seems to be remarkably constant).!7 The relative fading 
times of drift and decay patterns could be judged by 
cross-correlating two antennas separated by a space 
correlation distance in the direction of drift. If the maxi- 
mum correlation (at the displaced time, of course) is 
large, the fading at a point is due mainly to the drift and 
the turbulent effect will be masked. If the correlation is 
very small, the main effect is turbulent decay. (In 
theory, of course, the turbulent effect can always be 
measured by running along with the drift.) Finally, the 
antenna patterns should be as narrow as possible to 
keep the correlation function as free of dissimilar con- 
tributions as possible. 


APPENDIX 


According to (17), the received first-order scattered 
field is proportional to the Fourier transform of fo(r) 
corresponding to the wave number K. (We shall use the 
notation Tx { fo(r) } for this quantity in this Appendix.) 
Thus the statement leading to (18) will be proved if we 
can show [see (8)] that 


(e— aa) TeV (— Vea) =.0 


1 1 
(e — at): Tr {2 (v —- v) Do + DV? =| =0. (36) 


(35) 


and 


We shall retain only first order terms in Ae and replace 
V by 2ko when it operates on Eo or Do =€0Eo. 
To prove (35), we write 


7G, R. Sugar, “Some fading characteristics of regular vhf iono- 
spheric propagation,” Proc. IRE, vol. 43, pp. 1432-1436; October, 
1955. 
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eat [(Eo-V)Ve+ (Ve-iko) Eo + Ve X (iko X Eo) ]. (37) 


€0 
It is easily shown that 


Tx{Vf@)} = — iKTx{fO} (38) 


so that 


Tx {v (= ve fa) 
ial [—(E)-K)K + (K+ ko)Eo + K X (ko X Eo) |Tx{ Ae} 


1 
= — Tx{ Ac} (Eo: K)(ko — K) 


€0 


aga Tx| Ae} (Eo: K) ki. (39) 


€0 
But by definition of #, ki=kow, and since (e—a@a@)-a4=0, 
(35) is proved. It appears then that the additional 
term in F(r), due to the fact that V-E#0, corresponds 
to a longitudinal excitation in the direction of scattering 
and gives no far field. Note, however, that the magni- 
tude of this excitation is not small compared with the 
one which does produce a far field, hence cannot be 
neglected on these grounds. 

The proof of (36) is similar and does not even require 
the use of the dyadic. 


1 1 
Te {2 (v a v) Denies —| 
€ ( 
1 
Tx{2(Ve-iko) Eo + EvV’e} 


aS) 7 ne} (20ers 


since K = k)— k; and | Ro| = | ki = ko. 
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Scanning Characteristics of Microwave Aplanatic Lenses* 
G. G. CLOUTIER} anp G. BEKEFIt 


Summary—The imaging properties of two types of solid dielec- 
tric, circularly symmetric lenses were studied at a wavelength of 
1.25 cm. One of these was a wide-angle aplanat proposed some time 
ago by Martin; the other was a modification of the two-point cor- 
rected lens recently discussed by Sternberg. The work was carried 
out for the purpose of finding their scanning characteristics and to 
compare their ability to produce good pencil beams of radiation. 
Measurements disclosed that in these lenses (having focal number 
of 1.5), a beam, of half-power width 1.25°, could be scanned up to 
about 20° of arc from the optical axis by displacement of the source. 
The presence of astigmatism prevented either lens from being em- 
ployed for wider scanning angles. In view of this fact, there appears 
_ to be little advantage in this case in using the two-point corrected 

lens over the simple aplanat which can be more easily designed. 


INTRODUCTION 


| HE RANGE over which a pencil beam of radia- 
“ee can be scanned by displacing the primary 

source transversely to the principal axis of a lens 
depends upon the degree to which the lens can be cor- 
rected for the off-axis aberrations. The amount of cor- 


and Friedlander.’ They are characterized by the fact 
that they are free from axial spherical aberration and 
that they fulfil the Abbe sine condition. Satisfaction of 
the latter condition implies that the optical system is 
free from ordinary coma;? that is, from all those types 
of coma which vary as the first power of the source dis- 
placement from the optic axis. By ignoring all higher 
powers of displacement, other than the first, it is indi- 
cated that coma will only be corrected for small dis- 
placements of the source. The answer as to how large 
these can be without causing bad deterioration of the 
diffraction image is best found by experiment. 
Aplanatism in the above sense, namely the satisfac- 
tion of axial stigmatism and of the sine condition, may 
be expected to give insufficient correction when the 
source is displaced through quite large angles and if the 
lens has a small focal number. Higher orders of coma 
could then become significant. This fact led Ruze® to 
develop another type of lens system with which we shall 


also be concerned here. He designed two-dimensional 
metal plate lenses corrected exactly for spherical aberra- 
tion and coma—astigmatism being absent in two-di- 
mensional lenses illuminated by line sources—at two 
extra-axial object points. These foci are well displaced 
from the axis of the lens and are situated symmetrically 
on either side of it, e.g., at angular distances +69. To 
insure good behavior of the lens over the whole region 
lying within —0@)<@<, Ruze refined his design by re- 
quiring that a third, namely an axial point, also be cor- 
rected exactly. This he did by varying the refractive in- 
dex across the lens. 

A generalization of Ruze’s method was given recently 
by Sternberg’ for the design of two-point corrected solid 
dielectric lenses. Several three-dimensional circularly 
symmetric models were made by him from solid poly- 
styrene and were found to have good scanning properties 
over the range of source positions contained within the 
two extra-axial foci. By keeping to homogeneous dielec- 
tric material of constant refractive index, Sternberg 
could not insure a priori good correction for his lens at 
and close to the axis in a manner similar to that used by 
Ruze in the three-point corrected systems. Hence, the 
scanning properties of this lens, if it is to be used in 
scanning, are left somewhat to chance and its good be- 
havior outside the two focal points is somewhat fortui- 
tous. Bekefit reexamined Sternberg’s method of design; 


rection that can be achieved is governed by the number 
of optical surfaces available to the designer. A simple 
lens with two appropriately shaped surfaces can at best 
_be free from spherical aberration and ordinary coma. 
Astigmatism, which becomes prominent at the larger 
angles of displacement, cannot be eliminated from a 
single thin lens having an aperture stop coincident with 
its plane; this is so, whether the lens consists of a number 
of cemented components or whether it has spherical or 
_aspheric surfaces. Elimination of astigmatism requires 
“more complex optical systems such as an extremely 
thick lens,! a system of separated lenses, or a lens with 
an aperture stop displaced away from it.? Due to certain 
difficulties, little use has been made to date of these 
methods of correcting for astigmatism at microwave 
_frequencies;* most workers limit themselves to the sim- 
_ ple two-surface lens. Of the latter there are two distinct 
types which have found much application. There are the 
wide-angle aplanats (that is, lenses corrected for spheri- 
eal aberration and coma) based on the design of Martin? 


. * Manuscript received by the PGAP, November 5, 1956; revised 
~ manuscript received May 28, 1957. This work forms part of a project 
on microwave optics that is supported at McGill University, in the 
Raton Electronics Research Laboratory, by the USAF, Cambridge 
"Res. Ctr., Cambridge, Mass..on Contract AF 19(122)-81. 
| + Dept of Physics, McGill University, Montreal, P. Q., Can. 
4 t Mass. Inst. Tech., Cambridge, Mass., formerly at McGill Uni- 
_ versity, Montreal, P. Q., Can. By - 
1. C. Martin, “Wide aperature aplanatic single lenses,” Proc. 
: : ., vol. 56, pp. 104-113; 1944. 
- tet Hooking «Wave ee ot aberrations,” Clarendon Press, 
. land, pp. 132-136; 1950. 5 
d i. y e rita dn. “On the two-element separated doublet, 
paper no. 36, Symposium on Microwave Optics, McGill University, 
- Montreal, Can.; June, 1953. : 
4G. Bekefi, McGill University, | 
Final Report to the USAF Cambridge Res. Ctr., 
-19(122)-81; March, 1957. 


: 


5 F. G. Friedlander, “A dielectric lens aerial for wide-angle beam 
scanning,” J. JEE, p. 111A, vol. 93, pp. 658-662; 1946. 

6 J. Ruze, “Wide-angle metal plate optics,” Proc. IRE, vol. 38, 
pp. 53-59; January, 1950; also see “Correction,” vol. 39, p. 697; 
June, 1951. ; ; 

7R.L. Sternberg, “Elementary Methodsin the Numerical Design 
of Microwave Dielectric Lenses,” Lab. Electronics Rep., Boston, 
Mass.; 1954. Also see, J. Math. Phys., vol. 34, pp. 209-235; January, 
1956. 
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by permitting the presence of small errors in the two 
extra-axial object points and by examining the magni- 
tude of the various aberrations, he was able to choose 
that design which appeared to give good correction not 
only at the extra-axial foci but also at a third axial point. 
The method of design, in many respects similar to 
Sternberg’s, will be outlined below. 

On comparing the two types of design discussed above, 
there appears little doubt that an antenna engineer, 
interested in a two-dimensional lens giving a fan shaped 
beam and irradiated by a line source, should turn to the 
latter design in preference to the ordinary aplanat if 
wide angle scanning is desired. The two- or three-point 
corrected lens will give two good extra-axial images; in 
the intermediate region of the scan, errors are also not 
too large. If the design is applied to constructing a three- 
dimensional circular lens with a view to obtaining a 
pencil beam, the question arises: due to the presence of 
astigmatism, can much improvement be achieved 
through the use of this system in place of the aplanat? 
Astigmatism causes a very bad deterioration of the 
image;® it is virtually independent of the shape of the 
lens surfaces and it increases as the square of the source 
displacement and inversely as the focal length. Hence 
the two lens types, if of equal size and the same focal 
length, would be afflicted by about the same amount of 
astigmatism. Although the aplanat may at the same 
time suffer from higher orders of coma at fairly large 
source displacements, an error which would be absent in 
the two- or three-point corrected lens, effects due to these 
lens errors may be completely overshadowed by astig- 
matism. Whether or not this happens it was one of the 
main considerations which prompted our investigations. 
The two types of lenses were constructed and designed 
to be of equal size and focal length. Measurements of the 
respective radiation patterns, made in such a way as to 
insure a pencil beam of radiation, showed that the range 
of scanning angles of the two lenses was about the same 
and that astigmatism was the major factor which 
limited their performance. It was therefore concluded 
that in the design of lenses having focal numbers in the 
range of 1.5, nothing is gained by employing the 
laborious methods required in the design of the two- or 
three-point corrected systems in place of the much more 
easily designed aplanats. 

The procedures used in the design of the above types 
of lenses are outlined briefly below. For the sake of 
brevity, the Martin type of wide angle aplanat will be 
denoted as Lens A; the modified form of Sternberg’s 
design as Lens B. 


Lens A—WIDE-ANGLE APLANAT 


Let S; and S; be the two lens surfaces whose shape is 
to be determined [Fig. 1(a) ]. A general ray proceeds at 
an angle y from the paraxial focus Fo which is situated 


8M. P. Bachynski and G. Bekefi, McGill University, Eaton 
Electronics Res. Lab., Tech. Rep. No. 35, February, 1955; also 
see, “Aberrations in circularly symmetric microwave lenses,” IRE 
TRANS., vol. AP-4, pp. 412-421; July, 1956. 
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(b) 
Fig. 1. 


at some distance po from the apex of surface S;. The lens 
is made from a homogeneous material of refractive in- 
dex » and has a thickness ¢; 7 and r are the angles of 
incidence and refraction, respectively, at the first surface 
Si, and (x, y) are the coordinates of the point of inter- 
section of the ray with surface S:. Freedom from axial 
spherical aberration requires path equality between all 
rays proceeding from Fy, namely 


—a + pl + p = tu — 1). (1) 


Freedom from ordinary coma is obtained by satisfying 
the sine condition which states that 


(2) 
for an image at infinity. F, is the paraxial focal length. 
From equations above, from Snell’s law of refraction and 
from the lens geometry, one finds that the differential 
form of the equation for the first surface S; is given by5 
1 dp sin B | 
p dy pcosB—1 


y = F, sin y 


where 6 = (i—r); the angles 6 and y are found to be re-. 
lated through 


(Fp — p) sin ¥ 
t(u — 1) + p(cos y — 1) 


sin (y — 8) ; 
u — cos (y — 8) 


1957 


Cloutier and Bekefi: Scanning Characteristics of Microwave Aplanatic Lenses 


_ subtended by the lens at the focus. If the maximum value 


of y, say Yo, does not exceed about 25°, coefficients 
higher than ds are not needed. Difficulties arise if the 
lens is wider than this; higher order coefficients are then 
required, which are not available due to difficulties in 
their evaluation. Trial and error methods of estimating 


- such coefficients are discussed by Martin. However, in 


™~ 


such a case, it is preferable and more accurate to deter- 
mine the equation of the surface by direct numerical 
integration of the differential equation as was done by 


- Friedlander. 


A lens 50 cm in diameter having a focal number of 
about 1.5 was made from solid polystyrene plastic. The 
design parameters are given in the first column of Table 


J. All lengths are normalized to a unit of one meter. A 


check on the accuracy of the lens design was made in 


_ the following manner: the departure from axial stigma- 


tism (i.e., freedom from axial spherical aberration) was 


- found by calculating the path length difference between 
_ rays proceeding from the focus at various angles y and 


the axial ray. This path length difference in meters is 


_ given in the first row of Table II. The errors are seen to 


be extremely small. Departures from the exact fulfil- 
ment of the sine condition (known as the offense 
against the sine condition, or osc) was obtained from 
the relation?® 


9 Martin set F, equal to unity and took this value as his unit of 


length. 
iat A. E. Conrady, “Applied Optics and Optical Design,” Oxford 


_ University Press, New York, N. Y.; 1929 
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An elimination of 8 between the two equations results TABLE I 
in a differential equation for p in terms of the angle y 
and the lens constants yp, F,, t, and po. Martin! solved tae (aries Lens A (ete 
x é . ens Parameters engths normalized Wide- che 
the equation by a power series expansion to one meter as unit of length) a woe are ‘ 
ecte: 
p= poll + ary? t+ ayttay+-+-) (3) pSSRLY raaaieal) 
: Diameter 2a 
and evaluated the coefficients ae, as, etc., in general form mien focal length Fp earl: hieesas 
in terms of the lens constants.° efractive index of material p 1.58900 1.58900 
Fete the-ch f ee 3 Rea ae have been Lens thickness ¢ (as measured along the 
, the shape of the other surface is found by principal axis) ; 0.14432 | 0.14432 
numerical ray tracing. Position of extra-axial foci F\: 

In the above design procedure, numerical values must ¥o —- 050000 
be assigned eventually to the parameters u, F>, t, po. For ’ number ne He hner 
j Ps s 5y 2 0 . . 

a given refractive index and focal length the choice of a 0.52638 | 0.45799 
values to be assigned to pp and ¢ still remains open. This a 023960 ated 
. . . ° . 6 2 ome 
leaves room for some judicious guessing and several trial qo — 1.58482 
designs must be made to discover where the surfaces S; i ag ayieieee 
and Sz meet, before one arrives at a lenshaving the cor- Shape factor X= —(Ri+R:)/(Ri—R:) : 
rect aperture dimensions. The degree to which the final end Haan eee eal ve ca She ag 1.10585 
system is truly aplanatic depends on the angular size 
ABE E cul 
y in radians 0 0.05 0.10 0.15 0.20 Q.25 0.30 0.325 
Path length error (in units = - 1 
of 10-* meter) 0 —1 0 0 1 0 0 (equals —200 for Lens B) 
osc X 10° —235 
0 —59 +0.7 +11 +21 —29 —113 (equals 12700 for Lens B) 


ee i 
osc = rs +—), (4) 


Dp Po 


This expression is a dimensionless quantity; Fy is the 
“focal length” at an angle y and is defined as y/sin y 
where y is the height of the ray above the axis and 
parallel to it. Y is a quantity which determines the 
amount of spherical aberration in the system; it is given 
in terms of the distance between the paraxial focus and 
the point where some incident ray at height y will cross 
the axis after passing through the lens. If the lens is 
stigmatic or very nearly so, as in the present case, then 
the contribution from the term Y/pp is negligible. It will 
be seen from the second row of Table II that the osc is 
very small. In fact, both the errors tabulated lie well 
within the accuracy to which microwave lenses can be 
machined. 


Lens B—THE Two-PoIntT CORRECTED LENS 


As was mentioned previously, the Sternberg method 
of design which corrects exactly for two off-axis object 
points was modified somewhat in order to insure good 
correction also at an axial point. Consider a symmetrical 
lens [Fig. 1(b)] of size A142=2a@ again made from a ma- 
terial of refractive index uw. Ff: and F: are the two extra- 
axial foci situated at the coordinate points (xo, yo), 
(*o, —o) with respect to the origin 0. With the source at 
F,, the plane wave whose wave front is CiC2, is required 
to emerge at some angle y to the optic axis. Having 
chosen the position of Fi(xo, yo), the requirement of path 
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equality of the rays [F:A1C:]=[/iA2C2] determines the 
angle (W) which the emergent plane wave will make with 
the axis. The slopes of the surfaces at Ai and Az are 
found next. This is done by considering each apex as a 
small prism and by using Snell’s Law for the two rays 
C,AiF; and C.A2F, assumed to traverse the prisms in- 
finitesimally close to the apices A; and A». By invoking 
symmetry conditions of the two halves of the lens about 
the principal axis, a unique answer for the slope of the 
surfaces at A; and A, results. From this knowledge the 
marginal focus Fo is chosen as the center of the polar 
coordinates (p, y), (g, Y) of the two lens surfaces. The 
equation of these surfaces is expressed as a power series 
in y. To reduce the complexity of the calculations, it is 
assumed that three terms of the series suffice to define 
the surface, 


p = pod + azy? + asy*) 
g = go(1 + boy? + buy*). (5) 


The lens thickness, ¢, is then given by (qo—po). Six 
parameters, namely po, do, @2, @4 be, b4, must be determined. 
Four quantities which can be substituted into (5) are 
known already: when y has its maximum value Yo, 
p=q=FyA,. The knowledge of the slopes at the surfaces 
A, and A» determines the values of 


A complete solution of the unknown parameters re- 
quires the tracing of one more ray through the system. 
For reasons of simplicity a ray BoBF; passing through 
the center of the second surface was chosen. In place of 
satisfying exact path equality between this ray and the 
ray C,AiF\, a small path error A, was allowed to be 
present. It is defined by 


A> = [Ci AiF, — BoBF;,]. (6) 


It was then found convenient to study variations of this 
path error as a function of another path error (Aj) 
which is allowed to exist between the marginal and axial 
rays, namely 


A; = [PiAiFo — PBF)|. (7) 


If the latter error were zero, the system would be free 
from spherical aberration for the axial and marginal 
rays. However, the paraxial rays could still deviate from 
path equality and this would depend on the radii of 
curvature of the surfaces near the axis. 

For a given preassigned path error Aj, A» and all the 
lens parameters were evaluated. The amount of spheri- 
cal aberration and the offense against the sine condition 
were then computed. This was repeated for various 
values of A; and the design chosen resulted in the most 
reasonable combination of lens errors. A lens was de- 
signed in this manner, and the final parameters are 
given in the second column of Table I. All physical 
distances are normalized to twice the aperture size 
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Fig. 2—Aberrations in the “three-point” corrected lens B as function 
of small variations in the shape factor X. The arrow indicates the 
most promising design, details of which are given in Table I 
[osc is the offense against the sine condition; Y is the magnitude 
of spherical aberration A; and Ae defined by (6) and (7); all lengths 
normalized to the unit length of one meter]. 


(2a = A;A2) which is taken as the unit of distance of one 
meter. The variations of the lens errors as a function of 
the design are shown in Fig. 2. All lengths are expressed 
in units of one meter. It was found convenient to plot 
these variations vs the shape factor of the lens defined 
as 


Ri + Re 


gee eS 
Puoars 


where Ri; and R: are the paraxial radii of curvature of 
the two surfaces. The graphs denoted by A: and A; are 
the path length errors as defined by (6) and (7), re- 
spectively. Y is the spherical aberration expressed as 
the distance between the positions of the paraxial and 
marginal foci; it is taken to be positive when the mar- 
ginal focus lies between the lens and the paraxial focus 
and negative when it lies outside this region. It should 
be remarked that this time Y has been calculated only 
for the marginal ray in contrast to the previous case 
where it was computed for intermediate rays also. A 
reference to the last column of Table II will also show 
that the present lens has about 200 times as much 


spherical aberration as the Martin aplanat. The offense 
against the sine condition (for the marginal ray only) is — 


also depicted in Fig. 2. It is seen to vary little with 
change in the shape factor and it is nearly 55 times 
greater than in the previous lens (see Table II). How- 
ever, it must be stressed that despite the considerably 
greater aberrations this lens shows on and near the axis 
as compared with the Martin aplanat, the errors remain 
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LENS-B 


y 


79 .8ems- 
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Fig. 3—The scanning arc followed by a primary source irradiating lenses A and B, respectively. V and H are the planes containing the vertical 
and horizontal focal lines, respectively; C is the central plane and represents the arc followed by the source thus insuring a pencil beam of 
radiation; NV is the line followed by the source and giving a symmetrical image that is sensibly free of coma. 


very small when dealing with centimeter waves. 

While the offense against the sine condition varies 
little with the shape factor, the spherical aberration 
changes very rapidly. It will also be noted that when 
the extra-axial image becomes “free” of image errors 
(A:=0), as is the case in Sternberg’s design, the axial 
spherical aberration Y is fairly significant. This was 
found to be particularly so when lenses having a smaller 
f number were designed by this method. Hence it would 
appear that in order to insure a good axial focus, it is 
advantageous to allow some residual error A; to remain 
in the extra-axial focus. When lenses with very small f 
numbers are desired with greatly displaced extra-axial 
foci, the present method of design is too crude and more 
terms are required in the series expansion (5). The com- 
putational work, however, increases immensely and 
such refinements in design are discussed by Sternberg. 

A polystyrene lens was constructed according to the 
design data given in the second column of Table I. The 
lens is a shallow miniscus with its concave side facing 
the source. This can be compared with the previous 
lens, which is found to be doubly convex. 


RADIATION PATTERNS OF THE LENSES 


The radiation patterns of the lenses were measured 
on an outdoor site at a wavelength of 1.25 cm. The lenses 
were mounted in a large metal screen which acted as the 
aperture stop and prevented energy other than that 
traversing the lens from reaching the image space. The 
primary source was an open waveguide; this insured 
virtually constant illumination over the lens surface. 
Radiation patterns were then obtained for various 
angular displacements @ of the source from the principal 
axis. That path traced out by the primary source was 
chosen [Figs. 3(a), 3(b) ] which gave in each case a pencil 
beam of radiation. In order that this may be achieved, 
the source must be moved along a line lying midway 
between the astigmatic focal lines, namely along an arc 
in the central plane of the optical system. This is the arc 
C shown in Fig. 3(a) and 3(b). Departures from this arc 
will cause an elongation of the radiation pattern in a 


plane corresponding to the particular focal line which 
the source approaches. Corrections for coma can only be 
made for a conjugate pair of object and image distances. 
Hence, with the image distance fixed, there must exist 
an arc in a coma corrected system along which a dis- 
placement of the source will produce a symmetrical 
image sensibly free from coma. Such an are [denoted by 
N in Fig. 3(a) and 3(b)], was found experimentally. If 
the source is displaced on either side of this arc, the 
image becomes asymmetrical and these asymmetries 
are nearly mirror images of each other. The reason for 
this is, of course, that coma changes sign when the 
source crosses the arc NV. If these arcs (V and C) are 
widely removed from each other, one is faced with the 
problem as to which arc to choose for purposes of scan- 
ning. If path N is chosen, the radiation pattern will be 
free from coma but elongated in one plane; if path C is 
chosen, the main lobe is almost rotationally symmetrical 
but asymmetry of the side lobe levels appears in the 
plane of scanning due to the presence of coma. For the 
purposes of comparing the scanning characteristics of 
lenses A and B, it was found necessary to move the 
source along the arc C since this arc was identical in 
both systems, astigmatism being the same. Since in 
Lens B the arc WN lies further away from the central 
plane than in Lens A, asymmetries in the pattern of 
this lens are somewhat greater. 

In Fig. 4(a) and 4(b) are shown radiation patterns of 
the two lenses for several angular source displacements 
6. Both lenses behave equally well on axis and computa- 
tions show that the patterns approach very closely the 
ideal aberration-free image obtained from a perfect 
uniformly illuminated lens. As the angle of scan in- 
creases, slight asymmetries, mainly in Lens B, are ob- 
served; they are due to coma. The filling in of the 
minima, the broadening of the main lobe, and the in- 
crease in side lobe levels at the greater angles 0, is the 
result of astigmatism.® At an angle of 30° these effects 
are so bad that the pattern would not be usable in most 
applications. The circles shown in the graphs represent 
data computed from the diffraction theory of optical 
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Fig. 4—Radiation patterns obtained from lenses A and B, respectively, as function of the scan- 
ning angle 6; a is the angle of diffraction; the source was displaced along the central plane. 
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Fig. 5—Scanning characteristics of lenses A and B, respectively. The uppermost curve shows the decrease in gain as function of 6, The middle 
curve gives the variation of the first side lobe levels relative to the central maximum, and the bottom curve is that of the half-power 


beam widths. 


aberrations.!! The predictions from theory are seen to be 
extremely good except at the large angles of source dis- 
placement. 

The measured loss in gain relative to that when the 
source is on axis, the variations of the height of the first 
side lobes, and the change in the half-power beam width 
of the main lobe are shown graphically in Fig. 5(a) and 
5(b), as a function of the scanning angle @. These and 
previous results indicate that there is little to choose be- 
tween the actual performance of the lenses. As men- 
tioned before, astigmatism is the predominant factor 
which limits the degree of scanning that can be ob- 


11 G. G. Cloutier, “Scanning characteristics of some microwave 
aplanatic lenses,” M.S. thesis, McGill University; 1956. 


tained; in fact, it is so pronounced at the larger angles 
that it is impossible to tell whether the more elaborate 
design of Lens B results in a system which is better cor- 
rected for coma than is the simple aplanat. 
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The Long Distance Horizontal Radiation Pattern of a 
High-Frequency Antenna* 
RICHARD SILBERSTEIN+ 


Bureau of Standards, to compare the radiation 

pattern of a directive high-frequency antenna ob- 
served at a great distance with that determined locally. 
The array consisted of two three-element vertical yagis 
spaced a half wavelength apart with the center of the 
array about a half wavelength above ground. The hori- 
zontal beamwidth of the array was 50° and the operating 
frequency near 13.7 mc. 

In the first series of tests, transmissions were made 
from Sterling, Va., at various azimuths in accordance 
with a prearranged schedule with reception by a number 
of foreign and domestic laboratories. The antenna array 
was rotated through 340° in steps of 10° once each day 
for two weeks. At the distant stations records were made 
of the field strengths of these emissions and simultane- 
ous emissions on a vertical doublet. Data for each day 
were scaled and smoothed and the field strength values 
on each azimuth expressed as a ratio of amplitude for the 
array to amplitude for the doublet. Medians for the test 
period were then taken for each azimuth and normalized 
to the value for the axis of structural symmetry of the 
array. These data were then plotted and compared 
with the local radiation pattern. 

In the second series of tests the antenna was used for 
reception of a distant station and rotated once a minute 
for about an hour in each test. Medians were taken of 
the data on each azimuth and the values normalized. 

In both series of tests there was no appreciable distor- 
tion of the median pattern by the ionosphere as judged 


(pss were conducted in 1953, at the National 


* Manuscript received by the PGAP, June 23, 1956; revised 


manuscript received, May 22, 1957. 
+ Natl. Bureau of Standards, Boulder, Colo. 


by comparison with the pattern measured locally, even 
under poor propagation conditions. Fig. 1 illustrates the 
long-distance pattern at 3650 km compared with the 
local pattern. 
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Fig. 1—Sterling, Va., yagi radiation pattern for reception at 
San Diego, Calif., February, 1953. 


It is concluded that many instances of unexpected 
station interference attributed to differences between 
the distant and local patterns of an antenna may be 
more properly interpreted as instances of unusually 
high values of field strength, which is a time-distributed 
quantity. Unusually low values of field strength from 
the station it is desired to receive would also be a factor. 

Details of the experiment and the analysis are to be 
found in a report.! The report also gives some informa- 
tion on instrumentation problems under the severe re- 
quirement of minimizing interference in the reception 
tests. 

1R, Silberstein, “The Long Distance Horizontal Radiation Pat- 


tern of a High-Frequency Antenna,” Natl. Bureau of Standards 
Rep. No. 3538; July 6, 1955. 
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Surface Fields Produced by a Slot on a Cone* 
G. HELD} anp G. HASSERDJIANt 


field distribution on the surface of a cone excited 
by aslot, and the degree of interaction between the 
slot and the tip of the cone. 

In order to measure the field distribution on the cone 
surface, a small rod scatterer was supported perpendicu- 
lar to the cone surface. By measuring the “secondary” 
field radiated by the scatterer, placed at desired points, 
the field distribution was determined. The field intensity 
at a point can be measured by detecting a “secondary” 
field caused by a small scatterer placed at the field point.! 
The experimental setup for these measurements is 
shown in Fig. 1. 

As a source of excitation on the cone surface, two 
types of resonant slots were considered. First, a slot of 
length \/2 was studied which was cut along an element 
of the cone. Fig. 2 shows the measured distribution of 
the normal electric field along various elements of the 
cone for such a slot. Since the electric field distribution 
along the elements between the slot and the tip indi- 
cated no standing waves, it was concluded that a longi- 
tudinal slot did not interact with the tip. 

Second, a A/2 slot was cut circumferentially and the 
field distribution on the cone surface was measured. 
Fig. 3 shows the measured distribution along various 
elements of the cone for such a slot located 2.81 from 
the tip. The ¢=0° curve represents the distribution 
along the element passing through the center of the slot. 
These curves show standing waves which indicate inter- 
action between the tip of the cone and the slot. 

In order to present the experimental results in the 
form of tip to slot interaction, the cone surface was un- 
folded into a planar surface with the ¢=180°-element 
of the cone constituting the edge of this plane. It was 
assumed then, that the normal electric field on this sur- 
face could be represented by the expression 


E, — f(0')Gi(Rr) + AG»(kro) (1) 


[te paper presents an experimental study of the 


where the first term represents the radiation from the 
slot and the second term, the reflection from the tip. 
For a half-wavelength slot the coefficient of G; would be 


cos [1/2 cos 6’ 


sin 6’ 


) 


6’ being the angle that the radius would make with the 
slot. In (1) 7 is the distance from the slot and 7» the 


* Manuscript received by the PGAP, January 15, 1957. 

{+ University of Washington, Seattle, Wash. 

1R. Justice, “A Reflection Method of Measuring Near Field Dis- 
tribution,” Ohio State Univ., Antenna Lab. Rep. No. 486-28; 
September 30, 1954. 
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Fig, 2—Field distribution due to elemental slot excitation. 
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Fig. 3—Field distribution along elements of cone. 


distance from the tip. It was deduced further, that the 
coefficient of Gz could be defined as A =TGi(kd), with T 
representing the reflection coefficient of the tip (d= the 
distance between the slot and the tip). 

Theoretical evaluation of the fields of a cone excited 
by sources on the surface has yielded expressions con- 
sisting of the sum of an infinite seriess of Hankel func- 


1957 


tions.? However, such series are slowly convergent, thus 
quite difficult to evaluate, particularly on the surface of 
the cone. Assuming that one could represent the sum of 
such series by single functions, combinations of Hankel 
functions of fractional and integer order were tried for 
G, and G» to fit the experimental data. The best fit to the 
experimental distribution was obtained with G,= 
H;3?(kr) and G2= H*(kro). Fig. 4 shows a comparison be- 
tween an experimental curve along the ¢=0° and the 
assumed analytic expression. The choice for the coeffi- 
cient A was arrived at by a number of trials. With this 
value for A and the definition of the reflection coeffi- 
cient, was found to be equal to 0.725e-‘@"/, Using this 
same value of the reflection coefficient, computed values 
of the field distribution along different elements were 
compared with the data of several other slot-to-tip dis- 
tances. The agreement between the assumed expression 
and the experimental results was good. 

* L. L. Bailin and S. Silver, “Exterior electromagnetic boundary 


value problems for spheres and cones,” IRE Trans., vol. AP-4, pp. 
5-16; January, 1956. 


Cheng: On the Simulation of Fraunhofer Radiation Patterns in the Fresnel Region 
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Fig. 4—Field distribution along element —¢=0°. 


It should be noted that these measurements were 
made on a 28° cone. However, the outlined method can 
be applied to cones of any angle. 


On the Simulation of Fraunhofer Radiation Patterns 
in the Fresnel Region* 
DAVID K. CHENG} 


Summary—Physical limitations on the size of obstacle-free test 
sites give rise to the need of making radiation-pattern measurements 
on high-gain antennas at a reduced distance. The general practice is 
to defocus the primary source along the principal axis of the antenna 
reflector by a small distance so that Fraunhofer patterns may be 
simulated in the Fresnel region. This note summarizes and compares 
three different approaches with which the proper amount of defocus 
may be determined. 


specified in terms of the characteristics of its radia- 

tion pattern in the Fraunhofer region. This is be- 
cause of the fact that Fraunhofer radiation patterns do 
not change with the distance from the antenna as long 
as the far-zone approximations are satisfied. Radiation 
patterns in the Fresnel region tend to be very complex, 
and they change considerably with distance. 

While there is no clear-cut boundary between the 
Fraunhofer and the Fresnel regions, a common and ac- 
ceptable criterion is that 2D?/\ represents a safe far- 
zone distance, where D is the maximum dimension of 


A ss PERFORMANCE of an antenna is usually 


* Manuscript received by the PGAP, February 27, 1956; revised 
manuscript received February 4, 1957. This work was supported by 
the Rome Air Development Center, USAF, under Contract No. AF 


30(602)-1360. 
{ Elec. Eng. Dept., Syracuse Univ., Syracuse, N. Y. 


the antenna aperture and ) is the operating wavelength. 
At a distance of 2D?/X the maximum path-length dif- 
ference between the contribution from the edge of the 
aperture and that from the center corresponds to \/16 
or 7/8 radians. In practice, an unobstructed, open space 
with a dimension of 2D?/) is often not available for test- 
ing high-resolution antennas. For example, the 2D?/n 
distance for a 20-foot antenna at 3 cm would be about 
1.53 miles. Higher gain requirements would demand 
even larger test sites. Unfortunately, calculation of 
Fraunhofer radiation patterns from measurements in 
the Fresnel region is not of practical value at the present 
time both because of the difficult and laborious process 
of extrapolation and because of the inherent difficulties 
in making accurate amplitude and phase measure- 
ments.! The need for the technique of testing micro- 
wave antennas at reduced ranges is, therefore, both real 
and urgent. 

A commonly used practice is to displace the primary 
source of the antenna assembly slightly from the focal 
position in a direction away from the reflector. The 


1A. F. Kay, “Far Field Data at Close Distances,” Final Rep. for 
Contract No. AF 19(604)-1126, Tech. Res. Group, New York, N. Y.; 
October, 1954. 
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present note summarizes three different methods with 
which the required amount of defocus may be deter- 
mined.” The results are plotted and compared. 


THE GEOMETRICAL APPROACH 


It is a well-known fact that for an effective point 
source of excitation, the best radiation pattern will be 
obtained from a paraboloidal reflector at a field point in 
the far zone when the source is located at the focal 
point of the reflector. Geometrically, this may be ex- 
plained by equal path length from the source to all 
points in an aperture plane by virtue of the inherent 
property of a focused paraboloid. If the field point is 
far enough away from the reflector, the path lengths 
from the aperture points to the field point will be again 
approximately equal, resulting in an optimum additive 
effect. When the field point lies in the quasi-near zone 
(Fresnel region), the path-length differences from the 
points in an aperture plane to the field point must be 
compensated in some way in order that the measured 
radiation pattern may approach the true far-zone 
(Fraunhofer) pattern. This is done by slightly defocus- 
ing the source along the reflector axis in the direction 
away from the reflector. Since the amount of on-axis 
defocus is the only adjustable variable here, one cannot 
expect to achieve equal path length for all points in the 
aperture plane. For simplicity, the conventional ap- 
proach is to make the path length from the source to the 
field point by way of the apex of the paraboloid equal to 
that by way of the points on the edge of the reflector. 

Refer to Fig. 1, which represents a cross section of a 
symmetrical paraboloidal reflector with focal length 
OF =f and defocused point source at F’, the above re- 
quirement is equivalent to making 


F'O +00’ = F’'A+ AB (1) 


Call AA’=D (aperture diameter), O’P=BP=R (dis- 
tance of measurement), FF’ =e (defocus distance). Eq. 
(1) can be reduced to give 


joer ol 


When (f/R)?<1, it is accurate enough to write (2) as 


Avi(N} 


Normalizing all quantities with respect to the focal 
length and introducing new notations e’ =e/f, R’=R/f, 
and D’=D/f, one can rewrite (3) as 


ie} 


The normalized amount of defocus needed is seen to in- 
crease when R’ decreases and when D’ increases. As R’ 
approaches infinity, e’ correctly goes to zero. 


2D. K. Cheng, “Microwave aerial testing at reduced ranges,” 
Wireless Eng. vol. 33, pp. 234-237; October, 1956. 
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Fig. 1—Geometry of on-axis defocusing arrangement for 
paraboloidal reflector. 


THe APERTURE-PHASE APPROACH 


The defocusing problem can also be approached from 
a consideration of the phase distribution in an aperture 
plane of the reflector together with the diffraction in- 
tegral for the field at a point in space. When the point 
under consideration isin the quasi-near zone of a parab- 
oloidal reflector, the normalized diffraction integral 
which gives the field pattern in a horizontal plane can be 
approximated as? 


1 
I(u) = {. F(r)~i#2?r? 18 J (ur) dr. (5) 
0 


In (5), 7 is the radial dimension of the aperture plane 
normalized with respect to D/2; u=(aD/\) sin 0, @ be- 
ing the azimuth angle; k=27/\; and F(r) is the circu- 
larly symmetrical amplitude illumination function over 
the aperture. The explicit exponential term is the Fres- 
nel-region contribution; terms above the second order 
are neglected. When R is very large, (5) reduces to the 
far-zone pattern function 


To(u) =f Fe rdourar (6) 


When the primary source is displaced from the focus 
of a paraboloidal reflector along the reflector axis in the 
direction away from the reflector with a view to simu- 
lating far-zone patterns in the quasi-near zone, there 
will be a relative phase variation over the aperture. It 
has been found that this phase variation referred to the 
center point can be approximated satisfactorily by* 


9 


r 


LT apRRe rege 7) 
=r 1 
(5) + 
Eq. (7) is exact for r=0 (center) and r=1 (edge of aper- 


ture). For other values of 7, the | 6| given by (7) is 
slightly too large; the error decreases when the (f/D) 


6 = — 2e 


8S. Silver, “Microwave Antenna Theory and Design,” M.I.T. 
ae ope McGraw-Hill Book Co., Inc., New York, N. Y., vol. 12, 
ch. 6; : 
yf Doc kK. Cheng and S. T. Moseley, “On-axis defocus character- 
istics of the paraboloidal reflector,” IRE TRANs., vol. AP-3, pp. 
214-216; October, 1955. 
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ratio of the reflector increases. The diffraction integral 
now becomes 


1 
I(u) = f F (r) et (8—D*r?/8R ly J (ur) dr. (8) 

0 
In order to simulate Fraunhofer radiation patterns in 


the Fresnel region, the exponent under the integral sign 
in (8) should be made to vanish. This yields 


enona). 
érel+(G)] g 


which checks with (4) wien = 
able to write 


or 


R/f>1. If it is desir- 


R’ = nD?/N 


(10) 
with \’=)/f, n a numeric, then (9) reduces to 
Santas) baler. 
’ r nL\D' = 


Eq. (11) shows that for a given value of D’, (€/\) 
plotted vs 1 gives a hyperbola in linear scales, and a 
straight line in log-log scales.‘ It is noted that for n=2 
(R=2D?/h), appreciable defocus is still necessary. 


THE ELLIPSOIDAL-REFLECTOR APPROACH 


The purpose of defocusing the primary source in the 
case of a paraboloidal reflector is to simulate far-zone 
radiation patterns at points in the quasi-near zone. In 
terms of geometrical optics, it is quite easy to see that 
this could be achieved by means of an ellipsoidal re- 
flector. If the primary source is placed at one of the two 
foci of an ellipsoidal reflector, the reflected rays will con- 
verge at the other. 

The equation in the xz plane of a cross section of an 
ellipsoidal reflector with focal lengths fi and f2 is 


Tera s 
FSAI y/1- =], (12) 
2 fife 
Subject to the condition 
ae 2 
(/ fey se (13) 
fife — 2fife 
(13) can be approximated as 
pag Hits (14) 
Afi fe 
which is the equation for a parabola of focal length 
Werceet (15) 


fi + fe 
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or 


1 1 1 


= ee ee 

ie aren? 
Hence, for reflected rays to converge at R= fe, the pri- 
mary source should be placed at z=f;, and 


(16) 


e=f,— f= ue (17) 
aS, 
or, in normalized form, 
gee eae On ree aa 
RU BOR Ria (R58 


Eq. (18) should be compared with both (4) and (9). e’ 
can also be expressed in terms of the focal lengths as 


(19) 


which is extremely simple. 

The basis of the ellipsoidal-reflector approach lies in 
the fact that an ellipsoidal reflector of focal lengths f; 
and f2 approximates a paraboloidal reflector of focal 
length f as given by (15) or (16). An examination of 
(13) shows that it implies the condition (1/8) (x?/fife)?<1. 
Now the maximum value of x is D/2 S2f,. This reduces 
the condition to 


ty zt 
=e 2 


1 


(20) 


which is undoubtedly true in practice. An ellipsoidal re- 
flector with focal lengths f; and fz has its semimajor and 
semiminor axes equal to (f1+f2)/2 (arithmetical mean) 
and -/fifz (geometrical mean), respectively; it ap- 
proaches very closely a paraboloidal reflector when (20) 
is satisfied. As an example, with fe=R=50f1 the maxi- 
mum error is less than 0.09 per cent. 


COMPARISON OF DEFOCUSING METHODS 


Curves plotting e’ vs R’ based upon (4), (9), and (19) 
from the three different approaches discussed above are 
shown in Fig. 2. It is seen that except for small values 
of R’, the required e’ from the geometrical approach is 
nearly the same as that from the aperture-phase ap- 
proach, both of which increase with increasing D’. The 
required e’ from the ellipsoidal-reflector approach is the 
smallest of the three methods and is independent of D’. 

A review of the geometrical approach reveals that 
there is really no plausible justification in requiring 
equal path length from the source to the field point by 
way of the apex and by way of the points on the edge of 
the paraboloidal reflector only; the path lengths by way 
of the intermediate points on the reflector would then 
all be longer. Besides, there is no guarantee that the rays 
emanating from F’ will be reflected to pass through the 
point P except the ray along the principal axis. The ap- 
proximation (7) used in the aperture-phase approach is 
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Geometrical Approach 
— Aperture-Phase Approach 


e' = e/f (Defocus Distance in Multiples of Focal Length) 


“—————  Ellipsoidal-Reflector Approach 


10 20 50 100 200 500 


R' = R/f (Distance from Antenna Aperture in 
Multiples of Focal Length) 


Fig. 2—Comparison of defocusing methods. 
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exact for r=0 and r=1 only; for 0<r<1, 4, given 
by (7), is numerically too large resulting in an e’ which 
is also too large. It can be shown that the maximum 
error in 6 introduced by (7) is 


(ee (2 ee ) (21) 
pe iy Vig : 
where g=(D/4f)*. For a reflector with ¢g=0.35 or D’ 

= 2.36, the maximum error is about 5 per cent. 

Although the geometrical approach and the-aperture- 
phase approach yield approximately the same results, 
the aperture-phase approach makes it clear that this 
method would not be useful when R is too small because 
it would then be necessary to include terms higher than 
the second order in the exponent that appears in (5); the 
geometrical approach gives no indication of this resctric- 
tion. It is believed that e’ in (18) derived from the ellip- 
soidal-reflector approach gives the most nearly correct 
results because the approximation implied by (13) is 
very good; it does not restrict its correctness only to the 
edge of the reflector. 

It should be noted that in all three methods the re- 
quired amount of defocus is not a function of the oper- 
ating wavelength and that diffraction phenomena are 
neglected. 


Asymptotic Expansion of the Diffracted Wave for a 
Semi-Infinite Cone* 
L. B. FELSEN} 


the wave scattered by the tip of a semi-infinite cone 
has been obtained to first order in (1/kr), where k 
is the free-space wave number and 7 is the distance from 
the cone tip (see Fig. 1, opposite). It is shown in this 
note how one may obtain higher order terms in (1/kr). 


[i a previous paper,! an asymptotic expression for 


I. SCALAR PROBLEM (NEUMANN TyPE) 


In the domain (6+6’) <(20)—7), which excludes the 
rays reflected from the cone surface according to geo- 
metrical optics, the Neumann type Green’s function 
G'(r, r’) satisfying the inhomogeneous wave equation, 


SOR GAR te) it = nor antl) (1) 


* Manuscript received by the PGAP, April 26, 1957. This work 
was performed in connection with a subcontract with Univ. Michigan 
on Purchase Order No. 154700. 

ie Microwave Res. Inst., Polytech. Inst. of Brooklyn, Brooklyn, 


ot ae: Felsen, “Plane-wave scattering by small-angle cones,” 
(sy, TRawNs., vol. AP-5, pp. 121-129; January, 1957. See (24) and 


with 0G’/00=0 at 6=4, radiation condition at r—>« 
finiteness at r=0, can be represented exactly as? 
et*lr—r' | 


Gr) = + G,'(r, r’). (2a) 


der | r—r' | 
A time dependence exp (—iw#) is assumed throughout. 
The diffracted wave term G,’ arising from the scattering 
by the cone tip is given by 


ots ee 
VJrr'G,! = — — > em cos m(o — ¢’) 
Agr m=0 
’ f ax x Fiz (Rr) Az (hr’) Ges’ , (2b) 
aed a Kz ™(cos 6)K,—™(cos 6’) T'(ix + m + 3) 


- (—1)""! cosh raT(ix — m + 3) 
5 (d/d60)K;-"(—cos 60) 
(d/d0.)Kz-™(cos 0) 


2 Thid., see (8). 


(2c) 


oJ te» 
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ore p ef A 


Fig. 1—Cone geometry. 


5 etal cam Dl te mag Em = - m> 0; € = 3, 


(2d) 


: For plane-wave incidence, r’— « <, and we may employ 
_ the asymptotic formula’ 


Pe 
Hz“) (kr’) 7, y/ etkr’—a/4) tar]? kr’ > 0, (3a) 
wkr' 


_ If the wave is to be observed far from the cone tip so 
_ that kr is allowed to become very large, the complete 
_ asymptotic expansion for H;,"(kr) is appropriate in the 
range | kr|>>|x| : 


j Hes (br)~ 4/ 2 eitkr—m/4) pxn/2 s ack ' 
wkr St (—2ikr) 
s (4x?+- 1) (4x?+-3?) -- - [402+ (2n—1)?] 


27"! 


(3b) 


(3c) 
ooye 


n! 


(P3) (+342) «+ (@®+4+0(n=1)); 


(ix, 0) =1. (3d) 


The integrand in (2b) decays exponentially with in- 
creasing | «| in the domain (06+4’) <(26.—7), so that 
_ the error incurred by employing the expansion (3b) over 
the entire range in x is exponentially small as kr 
and therefore, may be ignored. One then obtains the fol- 
lowing asymptotic representation for the diffracted 
wave: 

i tk (rtr’) 2 


G,’ ~ —— 
Qr*k rr’ 


A,(9, 6; >, $') 


a 


n=0 


‘where the coefficients A, are given by 


7 > Em CoS m(o — ¢’) f dx x e?"(ix, n)Ges’. (4b) 
m=0 0 


3 W. Magnus and F. Oberhettinger, “Special Functions of Mathe- 
matical Physics,” Chelsea Pub. Co., New York, N. Y., p. 22; 1949. 
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Since the coefficients (ix, m) and (ix, m—1) are related 
via (3d) by 


a 
(ix, n) = eh [x + 3+ n(n — 1)] Gx, n — 1), 


1 e—N Mes Dies ena 


(5) 
and the K,™ functions satisfy the differential equation 
[Zm(8) — (4 + x°)]Kz"(cos 8) = 0, 


d m? 


Ln(@) = — sin 


sin 6 d@ do 


(6) 


? 
sin? 6 


one may construct the coefficient A, in (4b) from a 
knowledge of An_1 by means of the formula 


A, = 


1 
Ez = 


n 


Q2 
== iN) Ane 
sin? 6 d¢2 ae | : 


n= 1 RIP 
: 0 
: == ia) (y) = 
sin 6 00 00 


(7) 


Lo(9) = (7a) 


It is assumed that the interchange of the order of dif- 
ferentiation and summation (integration) necessary to 
deduce (7) from (4b) is permissible. Thus, all coefficients 
in the asymptotic formula in (4a) can be obtained from 
successive application of (7) if Ao is known. For small- 
angle cones (6) ~7), the coefficient A» in (4b) previously 
has been evaluated approximately to O[(#—6o)?] as! 


az — 00 {1+ cosé@ cos @’ 

“| 2 | hreseores cos 6’) 
2 sin @ sin 6’ 

(cos 6 + cos 6’)8 


Ao = 


Dah 6) | (8) 


II. VEcTOR PROBLEM (ELECTROMAGNETIC) 


The dyadic electric and magnetic Green’s functions 
z(r, r’) and y(r, r’), respectively, for a perfectly con- 
ducting semi-infinite cone can be expressed in terms of 
two scalar functions S,’(r, r’) and S;,’’(r, r’)5. In the 
domain (6+6’) <(20)—7), the total Green’s functions 
can be represented analogously to (2a) as the sum of the 
free-space Green’s function plus a diffracted wave term. 
The diffracted wave portions S;,’ and S;,’’ of Si’ and 
S,’’ contain a “residue contribution” which is 0(1) in 
kr and (kr’), plus an additional term to be denoted by 
Sis’ and S,,’’, respectively 36 


“ goa ss 
Sa (fet he verre me ye Em cos m(d — ¢’) 


T m=0 


2 x 
° dx Hz kr) Hix (Rr) Gos’, (9a) 
J Page : 


4 Felsen, op. cit., see (24). 
5 Ibid., see (32). 
6 Tbid., see (35). 
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£ cae yy 
Sip (at) = rr’ Pe > em cos m(p — ¢’) 


T m=0 


© x 
/ 
fiw e+ Hie (kr) Hiz™ (kr’)Ges, (9b) 


ri 
where Gy,’ is given in (2c). Ge, has the same form as Go,’ 
except that one replaces 


(d/d0o) Kz-”"(—cos 60) 
(d/d09)Kz-™(cos 60) 


K,-™(—cos 60) 


K,-™(cos 4) v9) 


As r'’— and subsequently kr>>1, the asymptotic ex- 
pansion of Sj,’ and S,’’ can be carried out as in Section 
I to yield 


ee —1 = Be 0, 0’; ) : 

Saif ies etk(r+r’) TD PE (10a) 
27k n=0 2ikr)” 
2 % ix, 2) 

Bn" = D> €m cos m(¢ — ¢’) i dx x e** \ Ges’, (10b) 
ed ae a 

and 

io —i = BO, 6; 6, 

es eik(rtr) oc AU (11a) 
Dak (2k 


B,’ =sameas B,,'’except that Go,’ isreplaced by Gg. (11b) 


The coefficients B,’’ and B,’ in (10b) and (11b) re- 
spectively, satisfy a recurrence relation identical with 
that for A, in (7): 


a : E ON ae 1)| mah 
Sy oe ae ete tLe 
BB," ae SA sin? 6 d¢? Bege 
pie A, es G2) 


In fact, the B, coefficients are the more basic since one 
may derive the A, coefficients from them through the 


relation [cf. (10b), (4b), and (6) ] 


a? | 
sap agi [Pa = Am = 12 19) 
sin 


Ex ae 


For the special case of an incident-plane wave polar- 
ized so that the electric vector is directed perpendicular — 
to the cone axis, it can be shown that the m=0 term in 
the representation of B,’ [see (11b) ] does not contribute 
to the scattered field. An approximate evaluation of the 
contributing portion of Bo’ and By’’ for small-angle 
cones has previously been shown to yield the following 
expressions to 0[(4—@:)?]:? 


6 6’ 
tan — tan —— 


E =| = fF ten see 
f cos (¢ — ¢’), a 
4 2 cos 6 + cos &’ ys, 


Tv 7 | leas 6 2 1 
By)!’ = =| | _ ——_————_._ — B,’. (14b) 
2 D cos6 + cos @’ 


The coefficients for ~>1 are obtained then, from suc- 
cessive application of (12). 

It is appropriate to point out at this time an error in 
algebraic sign which occurs in (45) of Felsen.t The 
first term on the right-hand side of (14b) corresponds 
to the above-cited reference and carries the correct 
algebraic sign. As a consequence of this error, (45), 
(48b), and (51) of Felsen must be corrected as follows: 

(45): multiply right-hand side by —1. 

(48b): multiply the ¢-independent term by —1. 

(51): replace 


(mf nt) (a 
Sins = €Os'0'COS* — sii? COS, sie 
2 2 7 2 2 


7 Tbid., see (44) and (46). 
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pany: January, p. 3 

Elliott, Robert, Is Founding Member of 
RANTEC Corp.: April, p. 175 

Feller, Robert, To Teach at University of 
South Carolina: January, p. 3 

Finke, H. A., Becomes Manager of Opera- 
tions, Polytechnic Research and De- 
velopment Corp.: January, p. 3 

Harrison, Bill, Joins RANTEC Corp.: 
April, p. 175 

King, D. D., Joins Electrical Communica- 
tions, Inc.: January, p. 3 

Krausz, Robert, Is Founding Member of 
RANTEC Corp.: April, p. 175 

Kurtz, Lou, Joins RANTEC Corp.: April, 

aS 

Rieti Don, Joins Systems Laboratories 
Corp.: April, p. 175 

Marsh, John, Elected Vice President of 
Systems Laboratories Corp.: April, p. 175 


Underberger, George, Joins Kearfott’s 
Microwave Group in Van Nuys, Calif.: 
April, p. 175 


Villard, O. G., Jr., Awarded 1957 Morris 
Liebmann Memorial Prize: January, p. 2 


AP TRANSACTIONS INDEX—5 


AIRBORNE ANTENNAS For Information Concerning 


Designers and producers of the finest 


in communication, navigation and ADVERTISING RATES 


countermeasures antennas, (Litera- 


Contact 
ture on request.) 


Mr. DELMER C. PoRTES 
DORNE AND MARGOLIN, INC. Jansky and Bailey, Inc. 
Westbury, Long Island, New York 1839 Wisconsin Ave., N.W. 

Washington 7, D.C. 


Exceptionally attractive openings for 


qualified microwave engineers. 


Telephone: Federal 3-4800 
f 


—_ eee eee ee eee 


PROCEEDINGS OF THE SYMPOSIUM ON 
ELECTROMAGNETIC WAVE THEORY 


Because of continued demand, consideration is being given to reprinting the Proceedings of the 
URSI Symposium on Electromagnetic Wave Theory, held at the University of Michigan in June, 
1955, and published in the IRE Transactions oN ANTENNAS AND PRopacartion, Vol. AP-4, No. 3, 
July, 1956. 


This 400-page volume contains 45 invited papers and five panel discussions, representing an 
unusually valuable and comprehensive cross-section of the latest theories of the world’s foremost 
authorities. Because of exceptionally heavy demand, stock copies were exhausted immediately 
after publication. 


To help assess the need for a reprint edition, your help is requested. If you would like to pur- 
chase a copy at a cost consistent with the size of the edition, please so indicate by writing to: 


Dr. K. M. Srecen 

Dept. of Electrical Engineering 
University of Michigan 

Ann Arbor, Michigan. 


The reprinting would be made possible through the cooperation and support of the University 
of Michigan Press and Engineering Research Institute. 


he solution to antenna problems begins when someone 

: “Let’s ask Kennedy!” 

A few of the many reasons why are shown on this page. 
These Kennedy antennas are setting new standards for all- 
weather reliability and versatility wherever they serve through- 
out the free world. 

Kennedy antennas come in many shapes, many’ sizes 
(the world’s largest scatter antenna is being built here). 
But whatever the type, and whatever the conditions under 
which it must serve, Kennedy can offer a design that fully 
measures up to specifications. And there are additional ad- 
vantages in Kennedy’s advanced construction techniques — 
like the extra ease in shipping and handling made possible 
by sectionalized aluminum construction, for example. 

Kennedy engineers are available to supervise installa- 
tions anywhere in the world. It’s a part of the complete, 
integrated service that is still another reason why more and 
more people in communications are “‘asking Kennedy” about 


antenna problems. 
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INSTITUTIONAL LISTINGS 


ed 


The IRE Professional Group on Antennas and Propagation is grateful tor 
the assistance given by the firms listed below, and invites application for 
Institutional Listing from other firms interested in the field of Antennas and 


Propagation. 
ANDREW CORPORATION, 363 E. 75th St., Chicago 19, Ill. ed 


Antennas, Antenna Systems, Transmission Lines, Development and Production. 


ANTLAB, 4950 North High St., Columbus 14, Ohio 
Antenna Pattern Range Systems—Recorders & Mounts. 


BLAINE ELECTRONETICS, INC., 14757 Keswick St., Van Nuys, Calif. 
Antennas, Paraboloids, Scale Models, Antenna Radiation Pattern Measurement Towers 


COLLINS RADIO COMPANY, Cedar Rapids, Iowa 
Antenna Design and Propagation Research Related for Airborne and Ground Communication Systems. 


rn 


DEVELOPMENTAL ENGINEERING CORP., 1001 Conn. Ave. N.W., Washington, D. C. and Leesburg, Ve 
Research, Development, Installation of Antennas and Antenna Equipment for Super Power Stations. t 


THE GABRIEL LABORATORIES, Div. of the Gabriel Co., 135 Crescent Road, Needham Heights 94, Mas: 
Research and Development of Antenna Equipment for Government and Industry. 


HUGHES AIRCRAFT COMPANY, Culver City, Calif. | 
Research, Development, Mfr.: Radar, Missiles, Antennas, Radomes, Tubes, Solid State Physics, Computers. 


I-T-E CIRCUIT BREAKER CO., Special Products Div., 601 E. Erie Ave., Philadelphia 34, Pa. 
Design, Development and Manufacture of Antennas, and Related Equipment. 


JANSKY & BAILEY, INC., 1339 Wisconsin Ave. N.W., Washington 7, D.C. 
Radio & Electronic Engineering; Antenna Research & Propagation Measurements; Systems Design & Evaluation 


MARYLAND ELECTRONIC MANUFACTURING CORPORATION, College Park, Md. 
Antenna and System Development and Production for Civil and Military Requirements. 


THE RAMO-WOOLDRIDGE CORPORATION, Los Angeles 45, Calif. 


TRANSCO PRODUCTS, INC., 12210 Nebraska Ave., Los Angeles 25, Calif. 
Res., Design, Dev., & Mfr. of Antenna Systems & Components for Missile, Aircraft & Ground Installations. 


WHEELER LABORATORIES, INC., 122 Cutter Mill Road, Great Neck, N. Y. 
Consulting Services, Research and Development, Microwave Antennas and Waveguide Components. 


The charge for an Institutional Listing is $25.00 per issue or $75.00 for four 
consecutive issues. Application may be made to the Technical Secretary, 
The Institute of Radio Engineers, 1 East 79th Street, New York alee ae. 


